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DIURNAL POLAR MOTION 


Paul McClure 


ABSTRACT 

An analytical theory is developed to describe diurnal polar motion in the 
Earth which arises as a forced response due to lunisolar torques and tidal de- 
formation. Doodson's expansion of the tide generating potential is used to 
represent the lunisolar torques. Both the magnitudes and the rates of change 
of perturbations in the Earth's inertia tensor are included in the dynamical 
equations for the polar motion so as to account for rotational and tidal 
deformation. 

It is found that in a deformable Earth with Love's number k = 0.29, the 
angular momentum vector departs by as much as 20 cm from the rotation axis 
rather than remaining within 1 or 2 cm as it would in a rigid Earth. This 20 cm 
separation is significant in the interpretation of sub-meter polar motion obser- 
vations because it necessitates an additional coordinate transformation in order 
to remove what would otherwise be a 20 cm error source in the conversion be- 
tween inertial and terrestrial reference systems. 
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DIURNAL POLAR MOTION 


1. INTRODUCTION 

Interaction of the lunar and solar tidal potentials with the Earth’s equatorial 
bulge generates torques on the Earth. Besides producing the well known phenom- 
ena of astronomical precession and nutation the lunisolar torques cause the rota- 
tion pole to travel within the Earth in a nearly diurnal epicycle which has a radius 
that varies from a few centimeters to a maximum of 62 cm. The diurnal motion 
of the rotation pole is superimposed upon longer period motions consisting mainly 
of the 14 month Chandler wobble and the annual and semiannual polar motion. 

Unlike the Chandler, annual and semiannual motions the diurnal motion of the 
pole has not been observed conclusively because of its comparatively small ampli- 
tude and high frequency. The astronomical methods for observing polar motion 
have uncertainties that are about equal to the amplitude of the diurnal polar motion 
and are able to produce pole positions only at 2 to 5 day intervals. Laser tracking 
of artificial satellites is now able to give pole positions at intervals of 6 hours 
[Smith et al., 1972] , thus providing an opportunity to observe diurnal motion of 
the rotation axis within the Earth. The satellite observations have noise levels 
of about 1 m and it is expected that observational uncertainties can be reduced 
to 10 cm in the future. 

In order to interpret polar motion observations with sub-meter noise levels 
it is necessary to model the diurnal motion of the pole. Woolard [1953] derived 
expressions for the diurnal polar motion in a rigid Earth. His results do not in- 
clude the effects of rotational and tidal deformation and the effect of the lunar and 
solar mean motion upon the coefficients in the diurnal polar motion terms is 
neglected. Melchior and Geeris [1968] use Doodson's [1922] expansion of the 
tide generating potential in order to obtain expressions for the lunisolar torques. 
The effect of mean motion of the disturbing bodies upon the diurnal polar motion 
is included as a second order correction. Their dynamical equations for the polar 
motion are for the case of a rigid Earth and do not include the effects of rotational 
and tidal deformation. 

The theory of diurnal polar motion presented here is for the ease of a deform- 
able Earth. The response of the Earth to deforming potentials is characterized 
by Love’s number k [Love, 1911]. Terms involving both the magnitudes and the 
rates of change of perturbations in the Earth's inertia tensor are included in the 
dynamical equations for the polar motion so as to account for the rotational ar»H 
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tidal deformation. Doodson’s expansion of the tide generating potential is used 
to represent the lunisolar torques, and mean motion of the tide generating bodies 
is included in the solution. The polar motion due to lunisolar torques is combined 
analytically with that due to rotational and tidal deformation in order to form a 
single set of coefficients for the diurnal polar motion. These coefficients along 
with the corresponding tesseral diurnal tidal arguments are arranged in tabular 
form so as to permit rapid computer evaluation of the diurnal polar motion at 
any instant of time and for a given Love number and set of astronomical constants. 

In addition to the results giving the motion of the rotation pole, solutions are 
obtained for the diurnal motion within the Earth of the angular momentum vector 
or principal axis of inertia. It is found that, in a deformable Earth, the angular 
momentum vector departs by as much as 20 cm from the rotation axis rather 
than remaining within 1 or 2 cm as it would if the Earth were rigid. The pre- 
cession-nutation theory represented in the Explanatory Supplement to the Ameri- 
can Ephemeris and Nautical Almanac [ 1961, p. 44] includes only the 1 to 2 cm 
rigid-Earth correction for the departure of the rotation axis from the angular 
momentum vector. An actual separation of 20 cm is significant in the interpre- 
tation of sub-meter polar motion observations because it necessitates an addi- 
tional coordinate transformation in order to remove what would otherwise be a 
20 cm error source in the conversion between inertial and terrestrial reference 
systems. Three alternative methods for making the additional coordinate trans- 
formation are discussed in detail. Each method makes use of tabulated coeffi- 
cients and arguments and can be readily programmed for use on an automatic 
computer. 

Woolard's theory of precession and nutation gives the best available repre- 
sentation for the direction of the Earth’s angular momentum vector in space. 

The diurnal polar motion theory presented here must be used in conjunction 
with such a precession-nutation model in order to give the orientation an 
observatory-fixed terrestrial reference frame. It is important that Woolard's 
results be understood in the context of the present development and, for this 
reason, a discussion of his method of solution is included. 
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2. POLAR MOTION DYNAMICS 

The rotational motion of a general mass distribution M is described by 


Liouville's equation 

L i = + e ijk "j »k 

(2.1) 

where 

H. = I tJ co. + h. 

(2.2) 


*ij = f (\\ S ij “ x i x j> dm 

J u 

(2.3) 


h i = f e iik x j \ dm 

(2.4) 


M 


Repeated indices indicate summation. The subscripts (i = 1, 2, 3) refer to a set 
of axes (x, y, z) having their origin at the center of mass and an angular velocity 
with components co. . L and H denote components of the net external torque 
and the angular momentum respectively. I . . is the inertia tensor. h. is the 
part of H. arising from motion relative to the x, y, z system. 

Although Liouville's equation is valid in any coordinate system, the axes 
shown in Figure (2.1) are especially useful for geophysical problems. The x, y, 
z system is attached to a set of terrestrial observatories in some prescribed 
manner. The x and y axes define the Earth's equator and the z axis is placed so 
as to remain nearly aligned with the rotation axis. 1 

The following perturbation scheme due to Munk and Macdonald [1960, p. 38] 
serves to simplify Liouville's equation. 


A + c n Cj 2 

C 13 


C 1 2 A + C 22 

C 23 

(2.5) 

C 1 3 C 23 

^ + C 33 


COj = IHj 0. 


(2.6) 

C0 2 = m 2 Q 


(2.7) 

co 3 - (1 + m 3 ) D 


(2.8) 
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Figure 2.1. Terrestrial Coordinate System. 



For the Earth, the quantities m., h./Cfl, and c. . /C are of order 10” 6 or smaller. 
Neglecting second order terms in the small quantities therefore produces first 
order equations for the polar motion that are exact to about one part in 10 6 . Sub- 
stitution of (2.5) through (2.8) into (2.1) gives 

Lj = Arrij Q + c J3 Q + hj 

+ m 2 (C - A) fl 2 - c 23 n 2 - h 2 fl (2.9) 

e « • 

L 2 = Am 2 fi + c 23 fl + h 2 

- ntj (c - A) a 2 + Cl 3 a 2 + n (2.10) 

L 3 = c 33 n + Cm 3 fi + h 3 (2.11) 

Equations (2.9) and (2.10) are written as a single complex equation in the form 

m = i < 2 1 (m - i/») (2.12) 
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where 


i =/n 

(2.13) 

m = nij + i m 2 

(2.14) 

cr ^ - (C - A) 12/ A 

(2.15) 

The polar motion excitation function 0 is 


i/r — iL c i c 

(C - A) Q 2 C - A (C - A) Q 


, h - i h 

(c ^ A) n ( C _ ^ n 2 

(2.16) 

where 


4> = + i 0 2 

(2.17) 

L — + i L 2 

(2.18) 

C = C 13 + 1 C 23 

(2.19) 

h = h x + i h 2 

(2.20) 

Equation (2.11) is rewritten as 


L, c h 

_ 3 33 3 

m3 _ co ” ~ co 

(2.21) 

The complex form of the polar motion equation makes solutions easier to 

visualize. For example, if there is no excitation then 0 = 
(2.12) is 

0 and the solution to 

i a t 

m 3 m Q e r 

(2.22) 


The motion is a counterclockwise circular path about the origin as shown in 
Figure (2.2). This "free" or "Eulerian" motion has a period of 10 months. Its 
amplitude and phase are determined by the initial condition 
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(2.23) 


m(0) = m 0 

where m 0 is a complex constant of integration. 


Y 



Figure 2.2. Eulerion Motion of the Pole. 


The first order expansion of the angular momentum vector is 

H = Afim + He + h (2.24) 

H 3 = Cfi(l + m 3 ) + c 33 n + h 3 (2.25) 

where 

H = Hj + iH 2 (2.26) 

The direction cosines of the rotation axis and the angular momentum vector are 


P| = ("V "V 


(2.27) 
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H 


1 


( 2 . 28 ) 


» 2 _ 

h| \cq’ cn 

in which second order terms are omitted. The complex numbers m and H/CQ 
represent equatorial projections of unit vectors directed along the rotation axis 
and the angular momentum vector respectively. 
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3. POLAR MOTION KINEMATICS 


It is necessary to define a set of coordinates to relate the terrestrial axes 
x, y, z to a set of inertial axes X, Y, Z. For this purpose Woolard [1953, p. 15] 
defines the Euler angles shown in Figure (3.1). The XY plane is that of the 
ecliptic at a prescribed epoch. The angle >p is the longitude of the equinox °P 1E 
and is measured in the ecliptic of epoch eastward from X. The obliquity of the 
ecliptic is denoted by 9 . The Earth's diurnal rotation is described by the angle 
4>, which is measured eastward from T 1E to the x axis. 



The rates of change of the Euler angles are related to the angular velocity 
components by Euler's kinematic equations, 
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(3.1) 


0sin0 = -&J 1 sin ^ - Wj cos ^ 

6 - ~ a Jj cos 0 + o> 2 sin <)i (3*2) 

4> - <^ 3 - 'j' cos 6 (3*3) 

Equations (3.1) and (3.2) are expressed in complex form as 

6 + i ^ sin 0 = •- Qm e l< ^ (3.4) 


A direction cosine vector u which is nearly aligned with the z axis can be 
described by the perturbations S<// and 3 <9 in the Euler angles as shown in Fig- 
ure (3.2). In order to relate the Euler angle perturbations to the direction 
cosines u. , it is convenient to introduce the node axes ^ , 77, £ shown in Figure 
(3.3). The £ axis points toward the equinox T and the 17 axis lies in the 
equator 90° to the east. The £ axis coincides with z. Neglecting second order 
terms in 30 and 3 6 , the node axis components of u are 

u* = - 30 sin 6 (3.5) 

u = 80 ( 3 . 6 ) 

V 

The Euler angle perturbations are defined in the sense 

30 = 0^ - 0 (3.7) 

U 

3 6 = 3-0 (3.8) 

U 

The x and y components of u are related to the node axis components by 

u ? + iu 7? = el<f> ( u x + U y^ ( 3 *9) 

Combining (3.9), (3.5), and (3.6) gives 

3 6 + iS0 sin 6 - - ie 1 ^ 1 (u^ + i ; u y ) (3.10) 
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4. LUNISOLAR TORQUES ON THE EARTH 


The torque exerted on the Earth by a point -mass disturbing body with geo- 
centric position r d and mass m d is 


3V 


L> — * £ i . . in. r j * i 
1 Ijk d dj 1 3 r 


(4.1) 


dk; 


The Earth's gravitational potential V is defined as the integral 


y _ f G dm 

■ JK-f 


(4.2) 


over all mass elements dm located at positions r within the Earth. V is ex- 
panded in terms of Legendre functions as 




n“2 


co n 


E E r] ^ 


cos m k, 

nm d 


n'2 m-1 


■+ S n m sin mA..) 


(4.3) 


The coordinates r d , 4> A , and A d of the disturbing body are shown in Figure (4.1). 
P m ( fx ) is the Legendre associated function of degree n and order m defined by 


P *(A0 = (1 -^ 2 ) m/2 — (P n (A^)) 
df-P- 


(4.4) 


where the P (/x) are Legendre's polynomials, 


1 d° (l ^ 2 - 1 ' 
P»=-T— ' 

n! dfjp 


(4.5) 
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z 



x d 

Figure 4.1. Disturbing Body Coordinates 


The torque components in the x, y, z system are 


t /BV\ . > m d / BV \ . , v 

L = - 1^—1 sin A. . + — — — — — sin®, cos A., 

cos 4> d \3\,/ d d 


3V \ i m<1 / BV \ . , . x 

— - cos A. + sin 0, sin A. 

*t>d) cos M 3 V 



(4.8) 


The second degree zonal term in V is larger by a factor of 10 3 than the higher 
degree zonal terms and the longitude dependent terms. Therefore, only the J 2 
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term is retained for the purpose of computing torques. This simplifies the ex- 
pressions for the equatorial torques L : and L 2 and results in zero torque about 
the z axis. The torques arising from the J 2 term are 


Gm £ 

= m. sin A. 

_ d i 



J 2 Pj (sin 0 d ) 


(4.9) 


L 


2 



m d cos A. d 



J 2 Pj (sin 4> d ) 


(4.10) 



(4.11) 


In order to use expressions (4.9) and (4.10) in the analysis of polar motion 
dynamics, the indicated functions of the disturbing body coordinates must be 
known as explicit functions of time. The time dependence of the lunar and solar 
coordinates is contained in Equation (A.22) which represents Doodson’s expan- 
sion of the tide generating potential. A formal spherical harmonic expansion of 
the tide generating potential is compared with (A.22) in order to determine the 
time dependence of the functions appearing in (4.9) and (4.10). 


The tide generating potential is 


0.^ V ftp 

r d [rj “ 


(cos y ) 


(4.12) 


n=2 


The local zenith angle y d of the disturbing body is given in terms of local lati- 
tude <f> and longitude A. by 


cos y d = sin <f> sin 4> d + cos 4> cos <£ d cos(A - A_ d ) 

Equation (4.13) is substituted into (4.12) to obtain 


(4.13) 


~ - n- , \ 

\ 1 \ < / r \ n 

-W„ P :< si "^P:( s in« d >cosm(\-\ d ) (4.14) 

" ^ ‘ 55 ’ \ 7 
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where 


W = - (n ~ m ) ! , (m = 1, 2, • ■ • , n) 
nm (n + m) ! v ' 


W „ = 1 

nO 


Linear independence of the spherical harmonic functions, 


(cos') 

r n P m (sin 0) -< Y 
" [s inj 




implies that the corresponding coefficients in (4.14) and (A.22) are equal. 


/l\ n+1 fcos't 

— \ W P m (sin0.)J V 

ir A nm nV d |sinj 




= (lY 1 T a ,/ cos l 

\ c d/ nn1Jd I s inj 


IT 


j- a), t - - (n - tn) 2 

Equation (4.18) is used to write (4.9) and (4.10) in the form 


Gm, r — i 

L 1 = Ve h L A 21 j COS < a, j t + 


um, ^ — i 

L 2 = 3 — m E a E J 2 ^ Sin(aJ j t +/? j ) 

C d i 

The complex form (2.18) of the torque is 

l = 22 a j e ’ i(V+v 


(4.15) 

(4.16) 

(4.17) 
Thus, 

(4.18) 

(4.19) 

(4.20) 

(4.21) 
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where 


Gm, 


S = ~ 3 — m E a lh A 21j 


(4.22) 


In evaluating the torques, the variations in J 2 due to rotational and tidal deforma- 
tion are neglected. From (5.16) with C 33 =0, 


J, = 


C - A 

m E a E 


(4.23) 


and (4.22) is written as 


Gm , 


A. = _ 3 -2. (c - A) A 21 . 


(4.24) 
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5. TIME VARIATIONS IN THE GEOPOTENTIAL AND IN THE 
INERTIA TENSOR DUE TO ROTATIONAL AND TIDAL 
DEFORMATION OF THE EARTH 

The Earth deforms as a result of tidal forces and the centrifugal force that 
arises from its spin about a shifting axis of rotation. Such deformations enter 
into the polar motion excitation function i (j of Equation (2.16) and into the right 
hand side of (2.21) by causing variations in the inertia tensor perturbations c , 
c J3 and c 33 . In order to write formulas for the inertia tensor perturbations, 
the parts of the Earth's gravitational potential due to rotational and tidal defor- 
mation are first found as functions of time. The inertia tensor perturbations 
are then related to the geopotential coefficients by Equations (D.13) through 
(D.17) and (D.26). 

The Earth's rotation axis moves relative to the terrestrial reference frame 
x, y, z of Figure (2.1). The resulting centrifugal force exerted on a mass ele- 
ment having the position vector r shown in Figure (5.1) is derivable from the 
rotational disturbing potential 


1 I — . — * L 2 

— \cj x r 

2 

The components of co are, from (2.6) through (2.8), 

nijil, m 2 fl, (1 + tn 3 ) Q 


(5.1) 


(5.2) 


in which Q is constant and the m. are of order 10" 6 . Neglecting second order 
terms in the m. , the rotational disturbing potential (5.1) is written as 


-jj- r 2 ^ 2 [1 - P 2 (sin </>)] 

O 


+ — r 2 fl 2 m 3 [1 - P 2 (sin 0)] 

3 

- — r 2 0 2 (nij cos A. + m 2 sin X.) P*(sin 4 > ) (5.3) 

3 
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z 



Figure 5.1. Generation of the 
Rotational Disturbing Potential 


The Earth responds by deforming so as to change its external potential by the 
amount 


V = 

V RD 



a E Q2 



P 2 (sin &) 




P 2 (sin <f>) 




cos \ + m 2 sin \) P*(sin <£) 


(5.4) 


The "secular" or constant part of the response takes place at a vanishingly small 
frequency and is therefore written in terms of the secular Love number k [Munk 
and Macdonald, 1960, p. 25] . The second and third terms on the right hand side 
of (5.4) represent response occuring at frequencies associated with the polar 
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motion and are given in terms of the "tidal -effective" Love number K [ Munk 
and Macdonald, 1960, p. 27] . 

The tide generating potential evaluated at r = a E is, from (A.22), 




P“ (sin <£) • 


The second degree zonal part of (5.5) contains one lunar and one solar term for 
which 

co. t + /S. =0 (5*6) 

J sec. * sec. 

The response to the secular part of U is given in terms of k rather than k. The 
external potential arising from response to the tide generating potential is 


Gm. /a y fa V 




n - 2 m=0 


a E V +1 K 


P: (sin <£) 



co. t + ft. + mA. + (n - m) ~ 

£ 


The external gravitational potential of the Earth is expanded as 



(5.7) 


(5.8) 
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The second degree zonal part of V is set equal to the combined zonal terms 
in (5.4) and (5.7) 




P 2 (sin 0) = 


\ k s a E (— 1 p 2 ( sin <£) 


2 /a ' 3 

<t i -i ra IB 


— k a| fl 2 m 3 [— ] Pj(sin <f>) 





k 




(sin 4>) • 



A 20j COS ("j t +^j) 


(5.9) 


sec. 

The inertia tensor perturbations corresponding to J 2 are found from Equations 
(D.13) through (D.17) and Equation (D.26), which take the form 


C- A 
■ + 

2C 33 " C 11 ' 

— = J 2 

(5.10) 

m E a E 

2m E 3 E 




c 22 ~ C 11 

= 0 

(5.11) 


C 1 2 

= 0 

(5.12) 


C X 3 

= 0 

(5.13) 


C 23 

= 0 

(5.14) 

C 1 1 

+ C 22 + C 33 

= 0 

(5.15) 
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Equations (5.15) and (5.10) are combined to give 


-^- + — = J 2 (5.16) 

m E a E ^ m E 3 E 

The first term on the right hand side of (5.9) is identified with the term 
(C-AjAr^a* in J 2 . This effectively defines the secular Love number in terms 
of C-A as 


k = ( 5 . 17 ) 

s 4 & 

The remaining terms on the right hand side of (5.9) are associated with the 
inertia perturbation c J3 . The part of c J3 due to rotational deformation is 


'33 


RD 


4 k ° 2m 3 4 

9 G 


(5.18) 


Equation (5.18) is written in terms of the secular Love number as 


' =— — (C - A) m 

33 rd 3 k v ' 3 

S 


The tidal contribution to c is 

33 


(5.19) 


c - ^_k d E A 

3 ^ K s A 20j 

C 3 


2 , m d a E 

H K 


3 C 3 

a }?) 


L 


A 20 j COS ("jt +/ 5 j) 


(5.20) 


Equations (5.11) through (5.15) give the remaining perturbations in the inertia 
tensor as 


''li 


= c 


RD 


22 


1 

C 


RD 


3 33 


RD 


(5.21) 
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( 5 . 22 ) 


1 

Cj < — ^55 — — C » » 

aa TD 3 00 TD 


The second degree tesseral part of V is set equal to a constant term char- 
acterized by the fixed coefficients C 2 i c and S 21c plus time varying terms from 
(5.4) and (5.7) due to rotational and tidal deformation. 




(5.23) 


In the case of a second degree tesseral potential, the inertia tensor relations are 


C 


21 


'13 


m E a E 


(5.24) 


S 


21 


'23 


m E 8 


2 

E 


The inertia tensor perturbations due to rotational deformation are 


(5.25) 


=13 =-P-(C-A)m. 

1 RD k 1 


(5.26) 
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(5.27) 


C 23 .= -(C-A)m 

/J RD k 1 

S 


In complex form, 


C RD = £ (C - A) m 


where 


Crd = Cl 3 rd + 1 C «rd 


The tidal perturbations in the inertia tensor are 


k m d a* 


13 


TD 


7^21 A 21 j sin (^j t +/3.) 


km d a I 


”23 


TD ^ 

d J 


22 A 2i j c ° s (v + ^) 


and these are written in complex form as 


km, at 


c - 

TD 


‘d 3 e 

c! 4 - 


A -i(«j t+^j-w/2) 

21 j 


where 


C TD “ C 1 3 + 1 C 23 

“td 2J td 


The constant inertia perturbations are 


c .3 0 '-"’ E a i c ,. t 

C 23„ = - m E a E S 21. 


In complex form, 


C 0 = " m E a i ( C 21„ + 1 S 21 


(5.28) 

(5.29) 

(5.30) 

(5.31) 

(5.32) 

(5.33) 

(5.34) 

(5.35) 

(5.36) 
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where 


/ 


/ 


c o = 


13 f 


+ 1 C 


23 - 


(5.37) 


The inertia tensor perturbations c c 22 and c do not enter explicitly 
into the excitation function (2.16) and therefore they have no direct effect on the 
polar motion. The lunisolar torques, however, depend ultimately upon every 
term in the geopotential. Before neglecting all but the J 2 term for the purpose 
of evaluating the lunisolar torques, it is necessary to know the magnitude of the 
rotational and tidal contributions to the C 22 and S 22 terms in the geopotential. 


The second degree sectorial term of V is set equal to a constant part plus 
a time varying part due to tidal deformation. 




• YL* 22 * 


P 2 (sin <f>) 


cos(a. + /3. ) cos 2\ 

1- sin(ajj -f /3. ) sin 2Al 


(5.38) 


In the case of a second degree sectorial potential. Equations (D.13) through 
(D.17) and (D.26) take the form 


C 22 - 


C 22 C 11 


4n V 


(5.39) 


S 


22 


^12 


2m E a E 


(5.40) 
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( 5 . 41 ) 


2 C 33 - C 11 ‘ C 22 = 0 


C 13 = 0 


C 23 = 0 


'11 + C 22 + C 33 ~ ® 


Relations (5*41) and (5.44) give 


C 3 3 ~ 0 


C 11 “ C 22 


From (5.39), 


C 

22 


22 


2m E a E 


The tidal perturbations in the inertia tensor are 


'22 


= - c 


2 km. 

d sl 


TD 


II 


TD 


J2 A 22j COS (^j t +^> 


(5.42) 

(5.43) 

(5.44) 

(5.45) 

(5.46) 


(5.47) 


(5.48) 


2km d a| 


"1 2 


TD c* 

d J 


22 * 22 > sin (oj. t + /?.) 


The constant inertia tensor perturbations are 


C 2 2 0 " C ll 0 - 2m E a E C 22 c 


C 12 Q “ ” 2m E a E S 22 c 


(5.49) 


(5.50) 

(5.51) 
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6. POLAR MOTION IN A RIGID EARTH WITH LUNISOLAR TORQUES 

The polar motion excitation function due to lunisolar torques is, from (2.16), 


0 L = — (6.1) 

(c - A> n 2 

where the complex torque L is given by (4.21). 

In a rigid Earth there is no rotational or tidal deformation to affect the 
inertia tensor, but constant products of inertia still enter into the excitation 
function (2.16). The excitation due to the constant products of inertia c Q of 
Equation (5.36) is 



C - A 


( 6 . 2 ) 


The differential equation (2.12) for the polar motion takes the form 


m = icr 


m 


c o 

C - A L _j 


^ A; ~i(co. t +/?. ) 
1 e J > 


r* (c - A> n 2 


(6.3) 


where A is given by (4.24). The general solution for the position of the rotation 
axis is 


m = m Q e 


“v + ^o_ + y th. 

C - A /_ -i AQ(o)^ + cr r ) 


-i("i t+^ ; ) 
e ’ 1 


(6.4) 


The complex constant of integration m 0 may be written in terms of an amplitude 
y 0 and a phase F 0 as 


m o = ?o eir ° 


(6.5) 


The axis of figure in a rigid Earth remains fixed relative to the x, y, z 
coordinate system. From Equation (E.9), the axis of figure is 


'Pf = 


C - A 


( 6 . 6 ) 
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The direction cosines of the angular momentum vector are found by substi- 
tuting the solution (6.4) for m into Equation (2.24). 


H 

Cfi 



icr 

e 


t 

r 



A 


+ 



iA. 
j 


-i(<y t+/3 ) 
e J J 


(6.7) 


The solutions for the rotation axis, the axis of figure and the angular 
momentum vector are shown in Figure (6.1). The axis of figure is fixed in the 
x, y plane. The constant m 0 represents the initial displacement of the Eulerian 
pole from the axis of figure. The Eulerian pole moves in a counterclockwise 
circular path about the axis of figure, completing one cycle in a period of 
27 t/ct = 10 months. Lunisolar excitation causes the rotation axis to move in a 
clockwise epicycle about the Eulerian pole position. Only one term of the sum- 
mation in (6.4) is represented by the circular epicycle shown in Figure (6.1). 
The frequencies a>. are grouped around the siderial frequency of 15.04107 de- 
grees per hour and lie in the range 

11.76554 deg. hr. -1 < co. < 17.69937 deg. hr.' 1 (6.8) 

The different terms alternately reinforce and cancel one another so that the 
rotation axis follows a path like the one shown in Figure (6.2). The amplitude 
of the nearly diurnal epicycle reaches about 62 cm when the principal terms 
are in phase. 


Ultimately it is necessary to locate the rotation axis, the axis of figure and 
the z axis of the "observatory-fixed" x, y, z system relative to inertially directed 
axes. The position of a particular axis in inertial space is specified by the Euler 
angles 6 and \p defined in Figure (3.1). The Euler angles for the angular 
momentum vector obtained by integrating Poisson's equations [Woolard, 1953, 
p. 34] are denoted by 6 and . The Euler angles Q f and for the rotation 
axis are given in complex form by 


ft 
vr r 


+ 


sin 6 - 6 


a . zp 


where, from the kinematical relationship (3.10), 


(6.9) 


Sd? r + i S \p T sin 6 = - ie 1 ^ 



( 6 . 10 ) 
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Figure 6.1. Polar Motion in o Rigid Earth 



Figure 6.2. Diurnal Polar Motion 
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The Euler angle perturbations for the axis of figure and the terrestrial z axis 
are given respectively by 


B6 { + i 8 sin 0 = 



( 6 . 11 ) 


8$ + i Si b sin 6 - i e 1 ^ 

z z cn 


( 6 . 12 ) 


Substitution of (6.4), (6.6) and (6.7) into (6.10), (6.11) and (6.12) gives 


8 6 1 + i8i/r r sin#--i 





i(cr t+4>) 

e r 


C - A j V" 1 A j 

\ C / “ An ("j. + °r) G , 


(6.13) 


50, 


■f i Si// f s i n 9 = l — 




i(^ t+<^> 
e r 



A. 

j 


C0(o;. +a r ) 


-i(<u t+/3 ~4>) 
e 1 J 


(6.14) 


Sfl • o / a ■ A i(cr t+<f>) 

ot> + l bvjj sin 0 = i — m e r 


+ i 


/ C 0 \ g i^ 
\ C " A / 


-L 


Cfi(^ +cr f ) 


“i(“i t+/3.“<£) 
e 1 1 


(6.15) 


The arguments oj. + t + /S. + and w._t + /?._ correspondingly to certain pairs 
of distinct indices j+ and j- are symmetric with respect to the Greenwich mean 
sidereal hour angle 0 in the sense that 

M 
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(6.16) 


M. ) + t +/3. + = (0 M +7T) + a. 

+ /3._ - (0 M + 77-) - a. (6.17) 

where a. denotes a linear combination of Doodson's standard variables defined 

j 

by Equation (B.13). The Euler angle perturbations (6.13) through (6.15) contain 
sums of the general form 



Ci). t +/S . ) 

J J ' 


(6.18) 


considered in Appendix F. The terms in (F.5) that contain sin(<£ M - 4 > ) involve 
products of the small difference (<p H ~ 4>) given by Equation (C.5) with factors the 
size of the Euler angle perturbations themselves. The sin(<£ M - <t>) terms are 
neglected and the resulting expressions are 



Si// r sin 6 = - 7 0 



cos (at + 4> + r o ) 



( 6 . 20 ) 


B6 f = - y Q — sin(c7 r t + 4> + r Q ) 


+ 



( 6 . 21 ) 


29 



Si// f sin 6 


80 - - 
Z 


+ 


Si sin 6 


7 ° 3 C ° S (°'r t + ^ + r o) 



' 0 ^-sin(cr r t +0 + r 0 ) 




( 6 . 22 ) 


(6.23) 


(6.24) 
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7. POLAR MOTION IN A DEFORMABLE EARTH WITH LUNISOLAR 
TORQUES AND TIDAL DEFORMATION 

The polar motion excitation functions due to lunisolar torques and to the 
constant part of the inertia tensor perturbations are 


and 


^L = 


iL 


(c - A> n 2 


'Pa - 


C - A 


(7.1) 


(7.2) 


respectively, where L is given by (4.21) and c 0 is given by (5.36). In a rigid 
Earth \p and i p Q are the only sources of excitation and their effects are con- 
sidered in Section 6. 

In a non-rigid Earth rotational deformation gives rise to the products of 
inertia c RD given by (5.28) and tidal deformation causes the products of inertia 
c^ given by (5.32). The excitation functions due to rotational and tidal defor- 
mation are, from (2.16), 


^RD “ 


"RD 


C - A 


A Crd 

IlC-A 


(7.3) 


/ _ C TD _A ^TD 

'^ TD = C - A C - A 

The rotation pole m satisfies (2.12), which takes the form 
m = i <r r (m - % - 0 L - - '/' TD ) 

Substituting for c^ from (5.28) into (7.3) gives 

, k i k ♦ 

rRD k flk 

S S 


(7.4) 


(7.5) 


(7.6) 
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Since ^ contains terms in m and m, it is necessary to modify the form of (7.5). 
Equation (7.5) becomes 


where 


m = 1 a Q m - 


i cr 

r 



( h + 'K + ^TD > 



(7.7) 


(7.8) 


Separate solutions of (7.7) are first obtained for each of the excitations 0 O , 
and i . These solutions are then combined with the free motion to form 
the complete solution. The solution due to 0 O is 


fixed products 
of inertia 



(7-9) 


The polar motion due to lunisolar torques is 


m 


Lunisolar 

torques 



i A. 

j 



A H (w. + cr 0 ) 


-»(<“« t +/3. ) 
e } 1 


(7.10) 


The coefficients A. are given by (4.24), 


Gm. 

*i'- 3 -r <c < 7 - u > 

c ; 

and the A 21j are related to Doodson's coefficients by Equation (A.23). The 
solution due to tidal deformation is 


32 



(7.12) 


tidal 

de f o rma t i on 


2 ] 


iCT r \i 


+/3. ) 


1 I 1 + ir v) ( “i ♦ <r o> 


where 


A,. =k 


m , 


,5 / _d 


n.i 


^ A„, . 


c 3/ C-A ,.fi 21 j 

d> 


(7.13) 


The frequency m is defined as 


n . ~ Q - a). 

J j 


(7.14) 


and is related to Doodson T s arguments by Equation (B.14). The free motion is 


m free = m 0 e 




(7.15) 


where m 0 is a complex constant of integration which is written in terms of its 
amplitude y Q and phase f o as 


m 0 = y 0 e 


iT„ 


(7.16) 


It is useful to make some observations about the separate solutions before 
combining them to form the complete solution. If the Love number k is zero. 
Equations (7.9), (7.10) and (7.15) reduce to the corresponding terms in (6.4) for 
a rigid Earth, and the polar motion (7.12) due to tidal deformation vanishes. 
Each coefficient A Xj is proportional to the angular rate n by which the fre- 
quency of the jth tidal component differs from the Earth's rotation rate. When 
M i =n > th e jth tidal component stands still in space and has no effect on the 
polar motion. The period of the free motion is lengthened by rotational defor- 
mation from 10 months to 277 / 0 ^ = 14 months. 


Tidal deformation produces polar motion having components with the same 
frequencies as those for the motion excited by lunisolar torques. The ratio of 
the jth coefficient in (7.12) to the corresponding coefficient in (7.10) is 


j • tidal deformation 
J. lunisolar torques 


" j i k 
n k 

s 


(7.17) 
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Equation (7.17) shows that the motion due to tidal deformation is smaller than 
that due to lunisolar torques and directed in the opposite sense. 

The solutions (7.9), (7.10), (7.12) and (7.15) are combined to give 


iCT t 

m = m 0 e + 



' AQ 2 


A j RA e 


i(“j t +/S. ) 


where 



The axis of figure is from Equation (E.9), 


(7.18) 


(7.19) 



c 

C - A 


(7.20) 


where c is the total of all perturbations in the inertia tensor. 

C = C 0 + C RD + C TD (7 ’ 21) 

The motion of the axis of figure due solely to rotational deformation is 


, c rd k 

\b c - — 1 = — m 

f RD C - A k 


(7.22) 


Substituting for m from (7.18) into (7.22) gives 


C - A 


+ 'Pr 


RD 


f - A ( + - m Q JL A. M e 

k . A-» k * An 2 JRA 

k 1 




(7.23) 
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The motion of the axis of figure due to tidal deformation alone is 



c 


TD 


C - A 


(7.24) 


Substituting for c TO from (5.32) into (7.23) and making use of (7.11) and (5.17) 
gives 


Pf 


TD 



. k A j -i("i 

i — ' e 1 

k , fl 2 (C - A) 




(7.25) 


The appearance of C-A in the denominator of (7.25) causes each term in to 
be much larger than the corresponding term in</> { . The jth terms of (7.25) and 
(7.23) have the ratio 150 


+ cT {7/i 

2i _k\ n 

n k ./ 

i 

The tidal effect is therefore oppositely directed from the rotational effect and 
about 300 times larger. In terms of displacement at the Earth’s surface, the 
amplitudes are of order 


'Pi 

j t tidal deformation 

Pf 

j T rotational deformation 


C - A, 


F l + 



| ~ 60 m 

tidal deformation 


(7.27) 


and 


Pf. 

it rotational deformation 


~ 20 cm 


(7.28) 
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The combined motion of the axis of figure is, from (7.23) and (7.25), 


= 




m o e 



i k 

Af la k 

r 9 


A. 


j AF 


+/ 3 . ) 

e J J 


(7.29) 


where 



(7.30) 


When k = 0, the solution (7.30) reduces to that given by (6.6) for a rigid Earth. 


The direction cosines (2.28) of the angular momentum vector are found by 
substituting the solution (7.18) for m into Equation (2.24). 


H 

cn 



[C - A 



m Q e 




+ 



i 

cn 2 



-i(o>. t + /3 ) 

A-.,, e J i 

j AM 


(7.31) 
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where 



The factor (l-k/k s ) = 2/3 in the summation term of (7.31) causes the diurnal 
terms for the angular momentum vector to be only 2/3 the size of the corre- 
sponding terms in the motion of the rotation axis. The result is that the angular 
momentum vector and the rotation axis are separated by as much as 1/3 x 62 cm, 
or 21 cm, instead of remaining within 1 or 2 cm of each other as they do in the 
case of a rigid Earth. 

The solutions for m, 0 f , and H/CH are shown in Figure (7.1). The fixed 
pole ¥ is defined by 


¥ = 



(7.33) 


An interesting consequence of Equation (2.24) is 


(7 ' 34 

so that the angular momentum vector, the axis of rotation and the axis of figure 
all lie in the same plane, for either a rigid or a deformable Earth. 

In the case of a deformable Earth, the Euler angle perturbations defined by 
Equations (6.10) through (6.12) are 
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Figure 7.1. Polar Motion in a Deformable Earth 
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-i (a). t+/ 3 .~(£) 

e J 1 


(7.36) 


S^ 2 + iS'/' z sin 0 




i(cr t+ 0 ) 

m Q e 




A 

1 AM 


-i(o>. t+£ -<£) 

e J J 


(7.37) 


The Euler angle perturbations contain sums of the form 


L 


A.e 

) 




(7.38) 


considered in Appendix F. The terms in Equation (F.5) that contain sin(0 M - 0) 
involve products of the small difference ( 0 M - 0) given by (C.5) with factors the 
size of the Euler angle perturbations themselves. The sin(0 - 0) terms are 
therefore of second order. When they are neglected the resulting expressions 
for the Euler angle perturbations are 
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C - A 
C 


s$ 

r 



7 0 sin(cr 0 t + 4 > + 


r o) 







S'/' sin 0 = - 


/C - A 1 



r 0 CO S (cr Q t + </> + 


r o) 



80 


f 


A 

C 



sin(cr o t + 4 > + r 0 ) 


+ 



k /A 

+ k (c 



+ A. 

3 AM 


C - A 


+ JifA 



cos or.. 

j 


(7.39) 


(7.40) 


(7.41) 
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y 0 cos (° 0 t + ^ + r 0 ) 


St p { sin 








bo 



1C - A' 



>o 


s i n ( a 0 t + 4> + r 0 ) 



St/» z sin 5 - 




y 0 cos(cr o t + 0 + r 0 ) 



The arguments a. are defined by Equation (B.13). 


(7.42) 


(7.43) 


(7.44) 
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8. COMPARISON OF THE PRESENT THEORY WITH WOOLARD'S 
DEVELOPMENT OF THE EULERIAN AND LUNISOLAR COR- 
RECTION TERMS 

Doodson's [1922] expansion of the tide generating potential is based upon 
Brown's [1905] lunar theory and Newcomb's [1898] theory of the sun. Woolard 
[1953] also uses the theories of Brown and Newcomb in his development of the 
luni solar precession and nutation. Since the solutions obtained in Sections 6 and 
7 are derived using Doodson's expansion, they may be expected to agree closely 
with Woolard's equations for the Eulerian and lunisolar correction terms. Wool- 
ard's paper is devoted mainly to a very elegant analytical integration of Poisson's 
equations in order to obtain the Euler angles and ^ of the angular momentum 
vector. No attempt is made here to comment on this part of his procedure. Of 
primary interest are the corrections S# z and S</» z that must be added to # H and 
p in order to obtain the position of the observatory-fixed z axis. 

Woolard assumes a rigid Earth, so his results must be compared with the 
solutions for a rigid Earth in Section 6, or equivalently, with those of Section 7 
after setting the Love number k equal to zero. In order to make a detailed com- 
parison, Woolard's development of the Euler angle perturbations and the polar 
motion coordinates is reproduced using the present notation. References to 
equation numbers in Woolard's paper are of the form (1), (2), etc. 

In a rigid Earth, the axis of figure remains fixed. Woolard takes advantage 
of this at the outset by aligning the z axis of the terrestrial coordinate system 
with the axis of figure. In Section 6, the axis of figure is given a fixed displace- 
ment c 0 /(C-A) relative to the z axis, and a corresponding term is carried along 
in the restatement of Woolard's equations for purposes of comparison. 


Woolard begins by differentiating his form (3) of Euler's kinematic equations. 
In the present notation, the kinematic equations are 


9 + ip sin# - - fim e 1 ^ 


and their derivative is 


— (9 ) + i — (di sin 9) - - Qme 1 ^ - Otnic^e 1 ^ 
dt v ' dt ' 


The time derivative m is replaced using the dynamical equation, 


m - icr 


m - 


iL 


c " A n 2 cc - A)J 


( 8 . 1 ) 


( 8 . 2 ) 


(8.3) 
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which is equivalent to Woolard's {6). Substitution of {8.3) into (8.2) gives 


— (9 ) , i 
dt ^ ; dt 


s in 8) = - + <fi) - e 1 ^ — + ie 1 ^ 

r A 


C() fi 2 


(8.4) 


The derivative <p in (8.4) is eliminated using Euler's third kinematic relation, 

= cj> + t p cos 9 (8.5) 

Since there are no torques about the z axis, the third dynamical equation reduces 
to 


a ^ = ft = constant (8.6) 

The polar motion m on the right hand side of (8.4) is eliminated using (8.1). The 
resulting expression is written as 


9 + i < p s in 0 = 


iL 

iA 

d 

eft" 

' Cft 

dt 


( 0 ) + 


— (di s in 9) 
dt K _ 


+ —(9 + i p sin 9) \p cos 9 - 


n 

e 1 


It is convenient to introduce the complex number * defined by 

Y = # + ii/< s in # 


In terms of y, Equation (8.7) is 


CO 


(8.7) 


( 8 . 8 ) 


(8.9) 


Equation (8.9) is equivalent to Woolard's (19). 

Poisson's equations for the motion of the angular momentum vector are 
stated as Woolard's Equation (30). In (8.9) the first and third terms on the left 
are neglected in comparison with the second to obtain Poisson's equations in 
the form 
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( 8 . 10 ) 


iL 

cn 


$ 


where the subscript H means that the corresponding Euler angles 9 and 0 
refer to the angular momentum vector. H H 

The terms "free" or "Eulerian" motion refer to the motion that would occur 
in the absence of the lunisolar torque L. Woolard's derivation of the Eulerian 
terms is discussed on pages 130 and 131 of his paper. The torque L is set equal 
to zero and it is assumed that 


4> = n 


( 8 . 11 ) 


in (8.9). Then 


X = 

which is the same as Woolard's 


• CQ • c o n2 id> 
1 T~V= i — ~e^ 

A A 


(52). The solution to (8.12) is 


( 8 . 12 ) 


C f) 

y = v 0 e { < cn/A) 1 (8.13) 

C **■ A 

in which y Q is a complex constant of integration and the assumption (8.11) is 
used for an approximate integration of the c Q term. Equation (8.13) may be 
written in terms of the polar motion m by making use of Euler's kinematic 
equation (8.1) in the form 


y = - (8.14) 

The polar motion solution corresponding to (8.13) is 


m = in e 1 [(C^/A) t-<£t0(o)] . _ c p. 

b C - A 


where 


% 


= - — e _i ^(0) 


n 


(8.15) 


(8.16) 
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The assumption (8.11) implies that the Euler angle 0 varies linearly with time. 
Thus 


0 = fit + 0(0) (8.17) 

Equation (8.17) is substituted into Equation (8.14) to obtain 


(8.18) 


which agrees exactly with the Eulerian terms in Equation (6.4). 

It is of interest that Equation (8.18) is in complete agreement with the 
Eulerian terms of (6.4) even though the assumption (8.11) is not used at all in 
deriving (6.4) while it is used twice in deriving (8.18) from (8.9). To see what 
is happening, Woolard's simplified differential equation (8.12) is written in polar 
motion form as 


m = i (cr f + 0 - 0) m - i 


C 0 H 


(8.19) 


Equation (8.19) differs from the Eulerian part of (6.3) by the inclusion of the 
extraneous term 


i(0 -.0) m (8.20) 

If the term iy (0-0) in (8.9) is retained at the outset, then the polar motion 
form is 


• . . c o n 

m = icr m - l 

r A 


( 8 . 21 ) 


which is identical with the Eulerian part of (6.3). Woolard T s initial neglect of 
the term iy( 0-0) in simplifying (8.9) is effectively cancelled by his assump- 
tion (8.17) of a linearly varying Euler angle 0 . 

In order to obtain the Euler angle perturbations, (8.10) is substituted for the 
torque term in (8.9) and the latter equation is written as 

= + < 8 - 22 > 
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From the definition (8.8), 


X-y H = se z + i8ijj z sin8 (8.23) 

where S# z and h\p t are the perturbations that must be added to the Euler angles 
of the angular momentum vector in order to obtain those of the z axis. The in- 
tegrals of the y and y terms on the right hand side of (8.22) are added to the 
integral of the Eulerian solution (8.13) to give 


S# + i£i// sin 6 - K + i — — e »(cf2/A)t 

1 1 c n 


+ i - C ° 

C - A 


(8 ' 24) 

In performing this integration, sine 0 is taken to be constant on the left hand side 
and the assumption (8.11) is used for the c 0 term. Equation (8.24) is equivalent 
to Woolard's Equations (53). Woolard's solution y of Poisson's equations is 
then substituted for y in Equation (8.24) in order to obtain his expressions (54) 
for S0, and S \p z . 

Woolard uses his solution for the Euler angle perturbations S 6 and S0 
in order to obtain the diurnal terms in the polar motion. From (2.24) the jth 
diurnal term m. in the polar motion is related to the corresponding angular 
momentum term by 


H. = Aflm. (8.25) 

) i 

in the case of a rigid Earth with c Q = 0. Substituting (8.25) into (6.12) gives 


m i =-i- 7-e -1 '^(S0 + i Si p sin 0) (8.26) 

j A z j z j 7 

which is the complex form of Woolard's (69). Woolard's Equations (70) for the 
main diurnal polar motion terms are obtained by substituting his solution (54) 
for the Euler angle perturbations into Equation (8.26). 
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The arguments in Woolard's diurnal polar motion expressions contain the 
Euler angle 0. The arguments ^ t +/3 j in the polar motion solution (6.4) con- 
tain the Greenwich mean sidereal time <* instead of 0 since, from (B.12), 

M 

"j t + /3. = (<p u + 7 T) + a. (8.27) 

Woolard's polar motion equations contain 0 in place of 0 M because the terms 
involving sin(0 M -0) are neglected in his Equation (54) for the Euler angle per- 
turbations upon which the polar motion is based. The same approximation in- 
volving the sin(^ M - 0) terms is made in deriving the Euler angle perturbations 
of Section 6, but there the polar motion is found first rather than being based 
upon the Euler angle perturbations as it is in Woolard's paper, and so the diffi- 
culty does not arise. The diurnal part of the Euler angle perturbations B8 z and 
50^ is of the form 

+ i S \p 2 sin 

i j 

Substituting (8,27) for the argument gives 


e )-L 


A. e 




(8.28) 



L K + “*■. 


s in 8} 





A. e 


(8.29) 


The diurnal terms in Woolard's Equations (54) are equivalent to (8.29) with 
(0 - <Pj neglected. Substituting (8.29) with the (0 - 0) term included into (8.26) 
gives 


(8,30) 

i j 

which contains 4^ in the arguments as it should. All terms involving the dif- 
ference (0 M -0) multiplied by the polar motion components or the Euler angle 
perturbations are of second order and are therefore negligible in numerical 
computations. However, it is more exact as well as more straightforward to 
use 0 . in the polar motion arguments instead of the Euler angle 0. 

The integral term in Equation (8.24) gives rise to the secular terms in 
Woolard's Equation (54) for 80 and in his polar motion equation (70). This 
integral term is effectively accounted for already by Woolard's use of the cor- 
rect Eulerian solution (8.18) in place of (8.13) and therefore should not be 
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included in (8.24) at all. Fortunately, the small size of the secular term makes 
it negligible in numerical computations. 

Of greater importance for numerical computations is the error inherent in 
Woolard's procedure for combining the complementary and particular solutions 
of the basic Equation (8.9). The nature of this error is readily determined by 
obtaining the particular integral of (8.9) with the help of expansion (4.21) of the 
lunisolar torque. Equation (8.9) is written as 


X - i(^ r + 4 >) Y = - ~ 

j 


' i ( ^ j t +/3j-<£) 


where cr^ is given by (2.15). The particular integral of (8.31) is 


(8.31) 


x 


-L 


iA. 

j 


A( w j + <r T ) 




(8.32) 


Woolard's procedure is to neglect the Y term in (8.31) and solve for Y to obtain 


v — * i A. 

= -2 

‘ ^ A f A\ i „ 




A(0 + o'') 


(8.33) 


The effect is to replace each tidal frequency go. in the coefficients of (8.33) by 
the frequency <j>, thus neglecting the motion of the sun and moon. The frequency 
4> is 


0 = Q- 0cos6 l = [1 + 0(10 n (8.34) 

so that solutions equivalent to Woolard's can be obtained to within 1 part in 10 7 
by substituting 


^ = n ( 8 . 35 ) 

into the coefficients of the rigid-Earth formulas of Section 6. Equations (6.4) 
and (6.7) with approximation (8.35) are used to obtain the coefficients given in 
Table (8.1) for the diurnal motion of the rotation axis and the angular momentum 
vector. The six largest terms in the diurnal polar motion given by Woolard's 
Equations (70) correspond to the terms listed in Table (8.2). The amplitudes in 
Table (8.2) agree to within the number of significant figures retained by Woolard, 
thus providing a check on the rigid-Earth polar motion theory of Section 6. 
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Table 8.1 


Coefficients for the Diurnal Motion of the Rotation Axis, 
the Axis of Figure and the Angular Momentum Vector in 
a Rigid Earth with Wollard’s Approximation; ^ . = & . 


CHEF F IC I ENTS IN SECONDS OP ARC 
SINE FOR X— COMPONENTS 
COSINE FOR Y-COMPONENTS 

INDEX TIDAL A RG . DIST. COEFFICIENTS OF ROTATION AXIS OF ANGULAR 



CODE NUMBER 

800 Y 

PHIM 

L 

LP 

F 

D 

OM 

AXIS 

FIGURE 

MOMENTUM 

I 

105,955 

M 

1 

-4 

0 

-2 

0 

-2 

-0.00000 18 

0.0 

-0.0000018 

2 

107*755 

M 

1 

-2 

0 

-2 

-2 

-2 

-0.0000075 

0.0 

-0.0000075 

3 

109.555 

M 

1 

0 

0 

-2 

-4 

-2 

-0.00000 46 

0.0 

-0.0000046 

4 

115.845 

M 

1 

-3 

0 

-2 

0 

-1 

-0.0000034 

Q.O 

-0.0000034 

5 

115.855 

M 

1 

-3 

0 

-2 

0 

-2 

-0,0000177 

0.0 

-0.0000177 

6 

117*645 

M 

1 

-1 

0 

-2 

-2 

-1 

-0 . 0000087 

0.0 

-0.0000087 

7 

117.655 

M 

1 

-1 

0 

-2 

-2 

-2 

-0 .00004 56 

0.0 

-0.0000455 

8 

118. 654 

M 

1 

-i 

1 

-2 

-2 

-2 

-0.0000 03 4 

0.0 

-0.0000034 

9 

119.445 

M 

1 

1 

0 

-2 

-4 

-1 

-0.00000 16 

0.0 

-0.0000016 

10 

119.455 

M 

1 

1 

0 

-*2 

-4 

-2 

-0.0000089 

o.o 

-0.0000068 

11 

124*756 

M 

1 

-2 

-1 

-2 

0 

-2 

0.0000021 

0.0 

0.0000021 

12 

125*745 

M 

1 

-2 

0 

-2 

0 

-1 

-0.0000295 

0.0 

-0.0000294 

13 

125.755 

M 

1 

-2 

0 

-2 

0 

-2 

-0.0001567 

0.0 

-0.0001562 

14 

12 6*55 f 

M 

1 

0 

-l 

—2 

-2 

“2 

0.0000026 

0.0 

0.0000026 

15 

126.655 

n 

1 

-1 

0 

- 2 

-1 

-2 

0.00000 18 

0.0 

0.0000018 

1^ 

126.754 

M 

1 

-2 

1 

-2 

0 

-2 

-0*0000025 

0.0 

-0.0000025 

17 

127.545 

M 

1 

0 

0 

-2 

-2 

-1 

-0.0000358 

0.0 

-0.0000357 

18 

127.555 

M 

1 

0 

0 

-2 

-2 

-2 

-0*0001892 

0.0 

-0.0001886 

19 

1 28. 544 

M 

1 

0 

1 

-2 

-2 

-1 

-0.0000023 

0.0 

-0.0000023 

20 

128.550 

M 

1 

0 

5 

—6 

2 

-6 

-0*0000130 

0,0 

-0.0000129 

21 

129.355 

M 

1 

2 

0 

-2 

-4 

-2 

-0.0000057 

0.0 

-0,0000057 

22 

133.855 

M 

1 

-3 

0 

-2 

2 

-2 

0.0000038 

0,0 

0.0000038 

23 

134,656 

M 

1 

-1 

-1 

-2 

0 

-2 

0.0000100 

0.0 

0.0000100 

24 

135.435 

M 

1 

1 

0 

—4 

0 

-2 

0.0000046 

0.0 

0*0000046 

25 

135. 635 

M 

1 

-1 

0 

-2 

0 

0 

0.0000069 

0,0 

0.0000069 

26 

135. 645 

M 

1 

-1 

0 

-2 

0 

-1 

-0.0002232 

0,0 

-0.0002224 

27 

135. 655 

M 

1 

-1 

0 

-2 

0 

-2 

-0.0011841 

0.0 

-0.0011802 

28 

135.855 

M 

1 

-3 

0 

0 

0 

0 

0.0000031 

0.0 

0.0000031 

29 

136.456 

M 

1 

1 

-1 

-2 

-2 

-2 

0.0000021 

0.0 

0.0000021 

30 

136.555 

M 

1 

0 

0 

-2 

-1 

-2 

0.0000064 

0.0 

0.0000064 

31 

136, 644 

M 

1 

-1 

1 

-2 

0 

-1 

-0.00000 18 

0.0 

-0.0000018 

32 

13 6. 654 

M 

1 

-1 

1 

“2 

0 

-2 

-0.0000112 

0.0 

-0.0000111 

33 

137.445 

M 

1 

1 

0 

“2 

-2 

-1 

-0.0000423 

0.0 

-0.0000422 

34 

137.455 

M 

1 

1 

0 

-2 

-2 

-2 

-0.0002250 

0.0 

-0.0002242 

35 

137*655 

M 

1 

-1 

0 

0 

-2 

0 

0.0000128 

0.0 

0.0000128 

36 

137.665 

M 

1 

-I 

0 

0 

-2 

-1 

-0.0000039 

0.0 

-0.0000039 

37 

138.444 

H 

1 

1 

1 

-2 

-2 

-1 

-0.00000 18 

0.0 

-0.0000018 

38 

138.454 

M 

1 

1 

1 

-2 

-2 

-2 

-0.0000 105 

0.0 

-0.0000105 

39 

139.455 

M 

1 

1 

0 

0 

-4 

0 

0.0000023 

0.0 

0.0000023 

40 

143.535 

M 

1 

0 

0 

—4 

2 

-2 

0.0000028 

0*0 

0.0000028 

41 

143.745 

M 

1 

-2 

0 

-2 

2 

-1 

0.0000033 

0,0 

0.0000033 

42 

143.755 

M 

1 

-2 

0 

-2 

2 

-2 

0.0000185 

0.0 

0.0000185 

43 

144.546 

M 

l 

0 

-1 

-2 

0 

-1 

0.0000025 

0.0 

0,0000025 

44 

144,55 6 

M 

1 

0 

-1 

-2 

0 

-2. 

0.0000 213 

0.0 

0.0000213 

45 

145.535 

M 

1 

0 

0 

-2 

0 

0 

0.0000358 

0.0 

0.0000357 

46 

145.545 

M 

1 

0 

0 

-2 

0 

-1 

-0.0011659 

0.0 

-0.0011621 

47 

145.555 

H 

1 

0 

0 

-2 

0 

-2 

-0.0061846 

0.0 

-0.0061644 

48 

145.755 

M 

1 

-2 

0 

Cl 

0 

0 

0.0000399 

0*0 

0.0000397 

49 

145.765 

M 

1 

-2 

0 

0 

0 

-1 

0.00000 66 

0.0 

0.0000065 

50 

146.544 

M 

1 

0 

1 

-2 

0 

-1 

-0.0000020 

0.0 

-0.0000020 

51 

146. 554 

M 

1 

0 

1 

-2 

0 

-2 

-0.0000 1 89 

0.0 

-0 • 0000 188 

52 

147.355 

M 

1 

2 

0 

-2 

-2 

-2 

0.0000034 

0.0 

0.0000034 

53 

147.545 

M 

1 

0 

0 

0 

- 2 

1 

-0.0000023 

0.0 

-0.0000023 

54 

147.555 

M 

1 

0 

0 

0 

-2 

0 

0.0000 80 6 

0.0 

0.0000803 

55 

147.565 

M 

l 

0 

0 

0 

-2 

-1 

-0.00001 76 

0.0 

-0.0000175 

56 

148.554 

M 

1 

0 

1 

0 

-2 

0 

0.0000054 

0.0 

0.0000054 

57 

152.656 

M 

1 

-1 

-1 

-2 

2 

-2 

0.0000023 

0.0 

0.0000023 

58 

153. 645 

M 

1 

-1 

0 

-2 

2 

-1 

0.000010 3 

0.0 

0.0000103 

59 

153.655 

M 

1 

-1 

0 

- 2 

2 

-2 

0.0000456 

0.0 

0.0000455 

<0 

1 54 . 65 6 

M 

1 

-1 

-1 

0 

0 

0 

-0.0000025 

0.0 

-0.0000025 
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Table 8 .1— (continued) 


COEFFICIENTS IN SECONDS OF ARC 
SINE FOR X -COMPONENTS 
COSINE FOR Y-COMPONENTS 


INDEX 

TIDAL ARG * 
CODE NUMBER 

DIST. 

BODY 

PHIM 

COEFFICIENTS 
L LP F 

OF 

D 

on 

ROTATION 

AXIS 

AXIS OF 
FIGURE 

ANGULAR 

MOMENTUM 

61 

155*435 

M 

1 

1 

0 

-2 

0 

0 

-0.0000028 

0.0 

-0.0000028 

tz 

155-445 

M 

1 

1 

0 

“2 

0 

“1 

0.0000323 

0.0 

0.0000322 

63 

155*455 

M 

1 

1 

0 

-2 

0 

-2 

0.0001748 

0.0 

0.0001742 

64 

155. 645 

M 

1 

-I 

0 

0 

0 

1 

-0.0000139 

0.0 

-0.0000139 

65 

155*655 

M 

1 

-1 

0 

0 

0 

0 

0.0004864 

0.0 

0*0004648 

66 

155. 665 

M 

1 

-1 

0 

0 

0 

“1 

0.0000 97 5 

0.0 

0.0000972 

67 

155*675 

M 

1 

“1 

0 

0 

0 

-2 

-0.0000028 

0.0 

-0*0000028 

68 

156.555 

M 

1 

0 

0 

0 

-1 

0 

-0.0000026 

0.0 

-0.0000026 

69 

156* 654 

M 

1 

-1 

1 

0 

0 

0 

0.00000 30 

0.0 

0*0000029 

70 

157*445 

M 

1 

l 

0 

0 

-2 

1 

-0.0000026 

0.0 

-0,0000026 

71 

157*455 

M 

1 

1 

0 

0 

-2 

0 

0.0000929 

0.0 

0*0000926 

72 

157.465 

M 

1 

1 

0 

0 

“2 

-l 

0 .0000 20 3 

0*0 

0.0000203 

73 

158*454 

M 

1 

1 

1 

0 

-2 

0 

0.0000039 

0*0 

0.0000039 

74 

1 61*557 

S 

1 

0 

-2 

-2 

2 

“2 

-0.0000069 

0.0 

-0.0000069 

75 

162*556 

s 

l 

0 

-1. 

-2 

2 

-2 

-0.0001692 

0*0 

-0.0001686 

76 

163*535 

M 

1 

O’ 

0 

-2 

2 

0 

-0.0000023 

0.0 

-0*0000023 

77 

163.545 

M 

1 

0 

0 

-2 

2 

-1 

0.0000327 

0-0 

0,0000325 

78 

163.555 

M 

1 

0 

0 

“2 

2 

-2 

-0.0000049 

0,0 

-0*0000049 

79 

163.555 

S 

1 

0 

0 

-2 

2 

-2 

-0.0028864 

0.0 

-0.0028769 

80 

163.557 

S 

1 

0 

-2 

0 

0 

0 

0.0000018 

0,0 

0*0000018 

81 

163.755 

M 

1 

-2 

0 

0 

2 

0 

0-0000043 

0.0 

0*0000043 

82 

164.554 

S 

I 

0 

1 

-2 

2 

-2 

0,0000241 

0.0 

0.0000241 

83 

164.556 

S 

1 

0 

“1 

0 

0 

0 

0.0000696 

0-0 

0*0000693 

84 

165.545 

M 

1 

0 

0 

0 

0 

1 

-0.0001723 

0*0 

-0.0001717 

85 

165.555 

S 

1 

0 

0 

0 

0 

0 

0.0027652 

0.0 

0.0027561 

8 6 

1 65.555 

M 

1 

0 

0 

0 

0 

0 

0.0059457 

0*0 

0.0059262 

87 

1 65 . 5 65 

M 

1 

0 

0 

0 

0 

-1 

0.0011785 

0.0 

0.0011747 

88 

165.575 

M 

1 

0 

0 

0 

0 

-2 

-0.0000253 

0*0 

-0.00002 52 

89 

1 66. 554 

S 

1 

0 

1 

0 

0 

0 

0.0000696 

0.0 

0.0000693 

90 

1 67*355 

M 

1 

2 

0 

0 

“2 

0 

0.0000043 

0.0 

0.0000043 

91 

167.553 

S 

1 

0 

2 

0 

0 

0 

0-00000 18 

0.0 

0.0000018 

92 

1 67.555 

s 

1 

0 

0 

2 

-2‘ 

2 

0.0001243 

0.0 

0.0001239 

93 

167*5 65 

M 

1 

0 

0 

2 

-2 

1 

-0.0000048 

0.0 

-0.0000047 

94 

1 67.575 

M 

1 

0 

0 

2 

“2 

0 

-0.0000023 

0-0 

-0.0000023 

95 

168.554 

S 

1 

0 

1 

2 

— 2 

2 

0.0000073 

0.0 

0.0000072 

96 

1 72 . 65 6 

M 

1 

-1 

-1 

0 

2 

0 

0.0000039 

0.0 

0.0000039 

97 

173.445 

M 

1 

1 

0 

-2 

2 

-1 

0*0000028 

0.0 

0.0000028 

98 

173.645 

M 

1 

-1 

0 

0 

2 

1 

-0.0000030 

0.0 

-0.0000029 

99 

173.655 

M 

1 

-1 

0 

0 

2 

0 

0*0000929 

0.0 

0.0000926 

100 

173- 665 

M 

1 

-L 

0 

0 

2 

-1 

0.0000184 

0.0 

0,0000183 

101 

174.456 

M 

1 

1 

-1 

0 

0 

0 

0*00000 30 

0.0 

0*0000029 

102 

174.555 

M 

1 

0 

0 

0 

1 

0 

-0*0000026 

0.0 

-0.0000026 

103 

175.445 

M 

1 

1 

0 

0 

0 

1 

-0.0000143 

0.0 

-0.0000142 

104 

175.455 

M 

1 

1 

0 

0 

0 

0 

0.0004864 

0*0 

0,0004848 

10 5 

175.4 65 

M 

1 

1 

0 

0 

0 

-1 

0.0000963 

. 0.0 

0.0000960 

106 

175.475 

M 

1 

1 

0 

0 

0 

-2 

-0*0000021 

0.0 

-0.0000021 

107 

175.655 

M 

1 

“1 

0 

2 

0 

2 

-0.0000075 

0.0 

-0.0000075 

108 

175.665 

M 

l 

“1 

0 

2 

0 

1 

-0.0000048 

0.0 

-0.0000047 

109 

175.675 

M 

1 

-1 

0 

2 

0 

0 

-0*0000028 

0*0 

-0.0000028 

110 

176,454 

M 

1 

1 

1 

0 

0 

0 

-0.0000025 

0.0 

-0.0000025 

111 

177.455 

M 

1 

1 

0 

2 

-2 

2 

-0.00000 20 

0.0 

-0.0000020 

112 

182.55 6 

M 

1 

0 

-1 

0 

2 

0 

0*0000053 

0.0 

0.0000052 

113. 

183.545 

M 

1 

0 

0 

0 

2 

1 

0*0000026 

0.0 

0.0000026 

114 

183.555 

M 

1 

0 

0 

0 

2 

0 

0.0000 807 

0.0 

0.0000805 

115 

183.565 

M 

1 

0 

0 

0 

2 

“1 

0.0000158 

0.0 

0,0000157 

116 

185.355 

M 

1 

2 

0 

0 

0 

0 

0.0000394 

0.0 

0.0000393 

117 

185.3 65 

M 

1 

2 

0 

0 

0 

“1 

0.0000079 

0.0 

0.0000079 

118 

185.555 

M 

1 

0 

0 

2 

0 

2 

0 .0002663 

0.0 

0*0002655 

119 

185.565 

M 

1 

0 

0 

2 

0 

1 

0.0001705 

0.0 

0.0001699 

120 

185.575 

M 

.1 

0 

0 

2 

0 

0 

0.0000358 

o_*.o. 

.0,0.0003.17, 

121 

185.585 

M 

1 

0 

0 

2 

0 

“1 

0.0000023 

0.0 

0.0000023 

122 

191. 655 

M 

1 

“1 

0 

0 

4 

0 

0.0000025 

0.0 

0.0000025 

123 

193.455 

M 

1 

1 

0 

0 

2 

0 

0.0000128 

0.0 

0.0000128 

124 

193.465 

M 

1 

1 

0 

0 

2 

“1 

0.0000025 

0.0 

0.0000025 

125 

193. 655 

M 

1 

-1 

0 

2 

2 

2 

0.00000 97 

0.0 

0.0000096 

126 

193* 665 

M 

1 

-1 

0 

2 

2 

1 

0.0000062 

0.0 

0-0000062 

127 

195.255 

M 

1 

3 

0 

0 

0 

0 

0.0000031 

0.0 

0.0000031 

128 

195.455 

M 

1 

1 

0 

2 

0 

2 

0.0000510 

0.0 

0.0000509 

129 

195.4 65 

M 

1 

1 

0 

2 

0 

1 

0.0000327 

0.0 

0.0000325 
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Table 8.1- (continued) 


COEFFICIENTS IN SECONDS OF ARC 
SINE FOR X -COMPONENTS . 
COSINE FOR Y -COMPONENTS 


INDEX 

TIDAL ARG. 
CODE NUMBER 

DIST. 

BODY 

PHIM 

COEFFICIENTS 
L LP F 

OF 

D 

m 

ROTATION 
..AXJS . 

AXIS OF 
FIGURE 

ANGULAR 

MOMENTUM.. 

130 

195*475 

M 

1 

1 

0 

2 

0 

0 

n. 0000069 

0.0 

0.0000069 

131 

1X3*555 

M 

1 

0 

0 

2 

2 

2 

0*0000002 

0.0 

0. 0000 OR 2 

132 

1X3*565 

M 

1 

0 

0 

2 

2 

1 

0.0000053 

0*0 

0.0000052 

133 

1X5.355 

M 

1 

2 

0 

2 

0 

2 

0.0000067 

0*0 

0.0000067 

13* 

1X5*365 

M 

1 

2 

0 

2 

0 

1 

0,000004* 

0.0 

0.0000044 

135 

IE3.455 

M 

1 

1 

0 

2 

2 

2 

0.0000020 

0.0 

0.0000020 


Explanation of symbols 

PHIM la the Greenwich mean aide rial time; ^ 

L, LP, F, D and OM are ‘Brown's ftindamental 
arguments; i' f F, D and D 


Constants 


k * 0 

K, = -71 r 552430 x 10 3 Julian century 1 
K e = -31'484150 x 10 3 Julian century" 1 
(C-A)/C * 3.272930 X 10" 3 
0 = 380^056 day" 1 
<C»/V = 2.343852 x 10* 

J* = 1.082645 x 10- 3 
(m^/nig) - 3,334320 x 10 « 
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The effect of Woolard's approximation (8.35) on numerical computations of the 
diurnal variation of latitude for Goddard Space Flight Center is shown in Figure 
(8.1). The diurnal terms in Equation (6.4) are substituted into Equation (9.14) in 
order to compute the latitude variation. The amplitudes differ by at most OV0008 
which is 3.7% of the diurnal polar motion amplitude or 2.5 cm at the Earth's 
surface. 


Table 8.2 

Comparison of Woolard's Diurnal Polar Motion Amplitudes 
with Those from Table 8.1 


Tidal Argument 
Code Number 

Amplitude from 
Woolard's Equations 
(70) 

Amplitude from 
the Column Headed 
"Rotation Axis "in 
Table 8.1 

135.655 

-0'.'0012 

-070011841 

145.545 

-0.0012 

-0.0011659 

145.555 

-0.0062 

-0.0061846 

163.555 

-0.0029 

-0.0028913 

165.555 

0.0087 

0.0087109 

165.565 

0.0012 

0.0011785 


The Euler angle perturbations in Table 8.3 are obtained from Equations (6.19) 
through (6.24) with a > = Q so as to make the results equivalent to Woolard's Equa- 

tion (54). The amplitudes of the largest terms in hd z and S< p z sin# from Table 
8.3 are compared with Woolard's Equation (54) in Table 8.4. The amplitudes 
agree to as many significant figures as Woolard retains, thus providing a check 
on the Euler angle perturbation theory of Section 6. 




LATITUDE (SECONDS OF ARC) 




Table 8-3 


Coefficients for Perturbations in the Euler Angles of the Rotation Axis, the Axis of Figure 
and the Angular Momentum Vector in a Rigid Earth with Woolard's Approximation; o>. = 

COEFFICIENTS FOR EULER ANGLE PERTURBATIONS 
IN SECONOS OF ARC 

COSINES FOR DELTA THETA 

SINES FOR DELTA PSt * S I N i THETA ) 

ROTATION AXIS OF TERRESTRIAL 

TIDAL ARGUMENT COEFFICIENTS DF AXIS FIGURE Z AXIS 

CODE NUMBERS L LP F D DM COS SIN COS SIM COS SIM 


10 5 • 955 

• 

4 

0 

2 

0 

2 

-0.00000001 

-0.0 000 00 o.l 

O.OQOOOISO . 

0.00000 IRQ. . 

.0. 00000180. 

O.Q 0000 IRQ 

10 7.^755 

* 

2 

0 

2 

2 

2 

-0.00000002 

-0.00000002 

0.00000752 

0,00000752 

0.00000752 

0.00000752 

109*555 


0 

0 

2 

4 

2 

-0, 00000002 

-0.0 000 000 2 

0.00000458 

0.00000458 

0.00000458 

O.OOQQ045R 

115.845 


3 

0 

2 

0 

1 

-0.00000001 

-0,00000001 

0.00000343 

0.00000343 

0.00000343 

0,00000343 

115*855 


3 

0 

2 

0 

2 

-0, 00000006 

-0.0000 000 6 

0.00001766 

0.00001766 

0.00001766 

0.00001766 

117- *45 


1 

0 

2 

2 

1 

-0.00000003 

-0.00000003 

0.00000867 

0.00000867 

0.00000867 

0.00000867 

117*655 


1 

0 

2 

2 

2 

-0*00000015 

-0.00000015 

0.00004547 

0.00004547 

0.00004547 

.0.00004547 

118. *54 


1 

-1 

2 

2 

2 

-0.00000001 

-0.00000001 

0. 00000343 

0.00000343 

0,00000343 

0.00000343 

119.445 


-1 

0 

2 

4 

1 

-0, 00000001 

-0.00000001 

0.00000164 

0.00000164 

0.00000164 

0.00000164 

119.455 


-1 

0 

2 

4 

2 

-0*00000003 

-0.00000003 

0 » 000008 B3 

0.00000883 

0.00000883 

0.00000883 

124,75* 


2 

1 

2 

0 

2 

0.00000001 

0.0000 0001 

-0.00000213 

-0.00000213 

-0.00000213 

-0.00000213 

125*745 

1X5.365 

2 

0 

2 

0 

1 

-0*00000008 

-0.00000011 

0.00002502 

0.00003386 

0.00002502 

0.00003386 

125-755 

1X5.355 

2 

0 

2 

0 

2 

-0.00000049 

-0.00000053 

0.00014949 

0.00016290 

0,00014949 

0,00016291 

12 f ♦ 5 5 6 


0 

l 

2 

2 

2 

0.00000001 

0.00000001 

-0.00000262 

-0.00000262 

-0.00000262 

-0.00000262 

126* 655 


1 

0 

2 

1 

2 

0*00000001 

0.00000001 

-0.00000180 

-0.00000180 

-0.00000180 

-0.00000180 

12 *. 754 


2 

-1 

2 

0 

2 

-0*00000001 

-0.00000001 

0.00000245 

0.00000245 

0.00000245 

0.00000245 

127.545 

1X3.565 

0 

0 

2 

2 

1 

-0.00000010 

-0.00000013 

0.00003042 

0.00004089 

0.00003042 

0,00004089 

127.555 

1X3*555 

0 

0 

2 

2 

2 

-0*00000059 

-0.00000065 

0.00018040 

0.00019676 

0.00018041 

0.00019676 

128.544 


0 

-1 

2 

2 

1 

-0* 00000001 

-0.00000001 

0.00000229 

0. 00000229 

0.00000229 

0.0000022 9 

128*550 


0 

-5 

6 

-2 

6 

-0.00000004 

-0.00000004 

0,00001292 

0.00001292 

0.00001292 

0.00001292 

129.355 


-2 

0 

2 

4 

2 

-0.00000002 

-0.00000002 

0.00000572 

0.00000572 

0.00000572 

0.00000572 

1 33.855 


3 

0 

2 

-2 

2 

0.00000001 

0.00000001 

-0.00000376 

-0.00000376 

-0,00000376 

-0.00000376 

134.656 


1 

1 

2 

0 

2 

0,00000003 

0.00000003 

-0.00000998 

-0.00000998 

-0.00000998 

-0 .00000998 

135.435 


-1 

0 

4 

0 

2 

0.00000002 

0.00000002 

-0,00000458 

-0.00000458 

-0.00000458 

-0.00000458 

135. *35 

195.475 

1 

0 

2 

0 

0 

0. 00000005 

0.0 

-0.00001374 

0.0 

-0.00001374 

0 .0 

135. *45 

195.4 65 

1 

0 

2 

0 

1 

-0.00000062 

-0.00000084 

0.00018989 

0.00025499 

0.00018989 

0.00025499 

135. 655 

195.455 

1 

0 

2 

0 

2 

-0. 00000371 

-0.00000404 

0.00112936 

0.00123109 

0.001L2938 

0.00123111 

135.855 

195.255 

3 

0 

0 

0 

0 

0.00000002 

0.0 

-0.00000622 

0.0 

-0.00000622 

0.0 

136.45 6 



1 

2 

2 

2 

o.ooooonoi 

0.00000001 

-0.00000213 

-0.00000213 

-0.00000213 

-0 .0000021 3 

1 3**555 


0 

0 

2 

1 

2 

0.00000002 

0.00000002 

-0.00000638 

-0.0000063B 

-0.00000636 

-0.00000638 

136, 644 


1 

-1 

2 

0 

1 

-0. OOOOOOO 1 

-0.00000001 

0.00000180 

0,00000180 

0.00000100 

0.00000180 

13*. *54 


1 

-1 

2 

0 

2 

-0.00000004 

-0.00000004 

0.00001112 

0.00001112 

0.00001112 

0.00001112 

137.445 

193. 665 

-1 

0 

2 

2 

1 

-0.00000012 

-0.00000016 

0.00003598 

0.00004841 

0.00003598 

0.00004841 

137,455 

193.655 

-1 

0 

2 

2 

2 

-0.00000070 

-0.00000077 

0,00021459 

. 0.00023389 

0.00021459 

0 • 00023389 

137,655 

193.455 

1 

0 

0 

2 

0 

0.00000008 

0.0 

-0.00002551 

0.0 

-0.00002552 

0.0 

1 37. * 65 


1 

0 

0 

2 

1 

-0.00000001 

-0.00000001 

O.OOOOCI393 

0.00000393 

0.00000393 

0.00000393 

138 .444 


-1 

-1 

2 

2 

1 

-0. 0000000 1 

-0.0000 000.1 

o.oooooieo 

0.00000180 

O.P0000180 

0,00.000.180 

138-454 


-1 

-1 

2 

2 

2 

-0.00000003 

-0.00000003 

0.00001047 

0.00001047 

0.00001047 

0.00001047 

139.455 

191 . 655 

-1 

0 

0 

4 

0 

0. 00000002 

-0.00000000 

-0.00000474 

0.00000016 

-0.00000474 

0.00000016 

143.535 


0 

0 

4 

-2 

2 

0.00000001 

0.00000001 

-0.00000278 

-0.00000278 

-0.00000278 

-0.00000278 



Table 8-3— (continued) 


COEFFICIENTS FflR RULER ANGLE PFR TUR RAT I DNS 
IN SECONDS OF ARC 

COSINES FOR DELTA THETA 

SINES FOR. flELTA PSI * S I N < THFTA ) 

RnTATlON AXIS OF TERRESTRIAL 

TIDAL ARGUMENT COEFFICIENTS OF AXIS FIGURE l AXIS 


CODE NUMBERS 

L 

LP 

F 

D 

OM 

CDS 

SIM 

CDS 

SIM 

cos 

SIN' 

143*745 

• 

2 

0 

2 

-2 

1 

O.OOOOOOOl 

0.00000001 

-0.00000327 

-0.00000327 

-0,00000327 

-0.00000327 

143.755 

• 

2 

0 

2 

-2 

2 

0.00000006 

0.00000006 

-0.00001848 

-0.00001848 

-0, 00001848 

-0,00001848 

144,546 

• 

0 

1 

2 

0 

1 

0.00000001 

O.OOOOOOOl 

-0.00000245 

-0.00000245 

-0.00000245 

-0.00000245 

*144 . 55 f 

* 

0 

1 

2 

0 

2 

0.00000007 

0.00000007 

-0,00002126 

-0.00002126 

-0. OO0O2 12 6 

-0 ,0000 2126 

145.535 

185.575 

0 

0 

2 

0 

0 

0.00000023 

0.0 

-0.00007131 

0.0 

-0.00007131 

0 • 0 

145.545 

185.565 

0 

0 

2 

0 

1 

-0.00000326 

-0.00000437 

0.00099214 

0,00133201 

0.00099215 

0,00133203 

145.555 

185.555 

0 

0 

2 

0 

2 

-0.00001937 

-0.00002111 

0.00589886 

0.00642977 

0.00589893 

0.00642984 

145.755 

185.355 

2 

0 

0 

0 

0 

0. 00000026 

0,00000000 

-0.00007900 

-0,00000049 

-0.0000790 0 

-n ,oooo(J049 

145.7 65 


2 

0 

0 

0 

1 

0.00000002 

0.00000002 

-0,00000654 

-0.00000654 

-0.00000654 

-0.00000654 

146. 544 

• 

0 

-1 

2 

0 

1 

-0. 00000001 

-0.0000 0001 

0,00000196 

0. 00000196 

0.00000196 

0 ,00000196 

146.554 

♦ 

0 

-1 

2 

0 

2 

-0.00000006 

-0.00000006 

0.00001881 

0.00001881 

0 « 00001 8R 1 

0 *00001 i: 81 

147.355 

, 

-2 

0 

2 

2 

2 

0. 00000001 

0.00000001 

-0,00000343 

-0.00000343 

-0.00000343 

-0 .000(10 343 

147.545 

183.565 

0 

0 

0 

2 

-1 

0.00000004 

-0.00000006 

-0.00001341 

0*00001799 

-0.00001341 

O.OOOOL799 

147.555 

183.555 

0 

0 

61 

2 

0 

0.00000053 

-0.00000000 

-0 • 000 16078 

0.00000016 

-0.0001607* 

0.0 0000(116 

147.565 

163.545 

0 

0 

0 

2 

1 

-0.00000005 

-0.00000007 

0.000014R3 

0.00002012 

0.00001488 

0.00002012 

148.554 

182.556 

0 

-l 

0 

2 

0 

0.0000000 3 

0.00000000 

-0,00001063 

-0.00000016 

-0.00001U63 

-0 ,00000016 

152.656 

♦ 

1 

1 

2 

-2 

2 

0.00000001 

0.00000001 

-0.00000229 

-0.00000229 

-0.00000229 

-0.00000229 

153. 645 

• 

1 

0 

2 

-2 

1 

0.00000003 

0.00000003 

-0.00001030 

-0. 00001030 

-0.00001030 

-0.00001030 

153.655 

177.455 

1 

0 

2 

-2 

2 

0.00000014 

0.00000016 

-0.00004351 

-0,00004743 

-0.00004351 

-0.00004743 

154.65 6 

176.454 

1 

1 

0 

0 

0 

-0.00000002 

0.0 

0.00000491 

0.0 

0. 0000049L 

n.o 

155.435 

175. 675 

-1 

0 

2 

0 

0 

-0.00000002 

0.0 

0.00000556 

0.0 

0.Q0U0O556 

0.0 

155.445 

175. 665 

-1 

0 

2 

0 

1 

0.00000009 

0 .000000 1 2 

-0.00002748 

-0.00003696 

-0.00002748 

-0.00003694 

155.455 

175. 655 

-1 

0 

2 

0 

2 

0.00000055 

0.00000060 

-0.00016666 

-0.00018171 

-0.00016667 

-0.0001*171 

155. 645 

175.465 

1 

0 

0 

0 

-1 

0.00000027 

-0.00000036 

-0.00008211 

0.00010991 

-0. 000.082 H 

0.00010991 

155. 655 

175.455 

1 

0 

0 

0 

0 

0.00000318 

0.0 

-0.00096957 

o.n 

-0.00096958 

0.0 

155.665 

175.445 

1 

0 

0 

0 

1 

0.00000027 

0.0 000003 7 

-0.00008292 

-0. 000 11138 

-0.00008292 

-0.00011138 

155.675 

. 

1 

0 

0 

0 

2 

-0.00000001 

-0.00000001 

0.00000278 

0.00000278 

0.00000278 

0.00000274 

156.555 

174.555 

0 

0 

0 

1 

0 

-o.ooooono? 

0.0 

0.00000523 

0.0 

0.00000523 

n.o 

156. 654 

174.456 

1 

“1 

0 

0 

0 

0.00000002 

0.0 

-0.00000589 

o.o 

-0.00 000589 

0.0 

157.445 

173.665 

“1 

0 

0 

2 

-1 

0.00000005 

-0.00000007 

-0.00001570 

0. 00002094 

-0.00001570 

0.C0OO2U94 

157.455 

173. 655 

-1 

0 

0 

2 

0 

0.00000061 

0.0 

-0.000185 15 

0.0 

-0.00018515 

0.0 

157,465 

173. 645 

-1 

0 

0 

2 

1 

0. 0000000 A 

0.00000008 

-0.00001734 

-0.00002323 

-0.00001734 

-0.00002 32 3 

158.454 

172. 656 

-1 

-1 

0 

2 

0 

0.00000003 

0.0 

-0.00000785 

0.0 

-0.00000785 

0.0 

161.557 

. 

0 

2 

2 

-2 

2 

-0.00000002 

-0.0000000 2 

0,00000687 

0, 00000687 

0. 00000687 

n.o 0000 6 87 

1 62.55 6 

1 68.554 

0 

1 

2 

-2 

2 

-0.00000053 

-0.00000058 

0.00016140 

0.00017588 

0,00016140 

0.00017588 

163.535 

167*575 

0 

0 

2 

-2 

0 

-0.00000002 

0.0 

0.00000458 

0.0 

0.00000458 

0.0 . 

L 63.545 

167.5 65 

0 

0 

2 

-2 

1 

0.00000009 

0.00000012 

-0.00002780 

-0.00003729 

-0.00002781 

-0.00003729 

163.555 

. 

0 

0 

2 

-2 

2 

-0.00000002 

-0.00000002 

0.00000491 

0.00000491 

0.00000491 

0.00000491 

1 63.555 

1 67.555 

0 

0 

2 

-2 

2 

-0.00000904 

-0.00000985 . 

0.00275300 

0.00300079 

0.00275303 

0.00300082 

163.55 7 

167.553 

0 

2 

0 

0 

0 

0.00000001 

0.0 

-0.00000362 

0.0 

-0.00000362 

0.0 



Table 8-3— (continued) 


COEFFICIENTS FOR EULER ANGLE PERTURBATIONS 
IN SECONDS DF ARC 

COSINES FOR DELTA THETA 

SINES FDR DELTA PSI * SI M< THETA) 


TI0AL 

CODE 

ARGUMENT 

NUMBERS 

COEFFICIENTS 
L LP F D 

OF 

OM 

ROTATION 

AXIS 

COS SIN 

AXIS DF 
FIGURE 

COS SIN 

TERRESTRIAL 

Z'AXIS 

COS SIM 

163.755 

167.355 

2 

0 

0 

-2 

O 

0.00000003 

0.0 

-0.00000850 

0.0 

-0.00000851 

o.p 

1 64. 554 

• 

0 

-1 

2 

-2 

2 

0.00000008 

0.00000008 

-0.00002407 

-o' . 00002407 

-0.00002407 

-0. 00002407 

164.556 

1 66.554 

0 

1 

0 

O 

0 

0. 00000046 

0.0- 

-0.00013R68 

0.0 

-0.00013869 

0.0 

1 65 • 545 

1 65.5 65 

0 

0 

0 

0 

-1 

0*00000329 

-0 . 00000442 

-0.00100293 

0.00134640 

-0.00100294 

0.00134642 

1 65 .555 

* 

D 

0 

0 

0 

u 

0.00000905 

0.00000905 

-0.00275609 

-0.00275609 

-0.00275613 

-0.00275613 

1 65.555 

• 

O 

0 

0 

0 

0 

0.00001946 

0.00001946 

-0.00592617 

-0,00592617 

-0.00592624 

-0.00592624 

• 

1 65.57 5 

0 

0 

0 

0 

-2 

-0 • 0.000000 B 

0.00000008 

0.000Q25.1.9 . 

-Q. 00002519 

0.00002519 

-0.0Q002319 

» 

173.445 

1 

0 

-2 

2 

-1 

0. 00000001 

-0.00000001 

-0.00000278 

0.00000278 

-0.00000278 

0.00000278 

• 

175.475 

1 

0 

0 

0 

-2 

-0. 00000001 

0.00000001 

0.00000213 

-0.00000213 

0.00000213 

-0.00000213 

♦ 

105.3 65 

2 

0 

0 

0 

-1 

0.00000003 N 

-0.00000003 

-0.00000785 

0.00000785 

-0.00000785 

0.00000785 

« 

185.585 

0 

0 

2 

0 

-1 

0.00000001 

-0.00000001 

-0.00000229 

0.00000229 

-0.00000229 

0.00000229 

t 

193.4 65 

1 

0 

O 

2 

-1 

o.oooaoooi 

-0.00000001 

-0.00000245 

0.00000245 

-0,00000245 

0.00000245 

♦ 

163.455 

1 

0 

2 

2 

2 

0.00000001 

-0.00Q00001 

-0.0000Q196 

0.00000196 

-0.00000196 

0.00000196 


Cn 

CD 


Explanation of aymbolg 


I* LP, F. D and OM are Brown's fundamental 
arguments; l t l’ , F. D and 0 


Constants 

k = 0 

= -71'552430 x 10 s Julian century -1 
- -3' '48 4150 x 10 3 Julian century -1 
<C-A)/C = 3.272930 x 10" 3 
n = 3603S56 day -1 
<VJ\) = 2.343862 x 10 4 
J a = 1.082845 x 10" 3 
= 3.334320 x 10* 



Table 8.4 

Comparison of Woolard's Euler Angle Perturbations 
with those from Table 8.3 


Code Numbers of 

Terrestrial z 

Axis Euler Angle 

Euler Angle Perturbations 

Symmetric Tidal 

Perturbations from Table 8.3 

from Woolard 

’s Equation (54) 

Arguments 

z 

h\p z s in d 

hd 

Z 

h\p z s in 6 

117.655 


0700004547 

0V00004547 


0700005 

125.755 

1x5.355 

0.00014949 

0.00016291 

0.00015 

0.00016 

127.545 

1x3.565 

0.00003042 

0.00004089 


0.00004 

127.555 

1X3.555 

0.00018041 

0.00019676 

0.00018 

0.00020 

135.645 

195.465 

0.00018989 

0.00025499 

0.00019 

0.00025 

135.655 

195.455 

0.00112938 

0.00123111 

0.00113 

0.00123 

137.445 

193.665 

0.00003598 

0.00004841 


0.00005 

137.455 

193.655 

0.00021459 

0.00023389 

0.00021 

0.00023 

145.545 

185.565 

0.00099215 

0.00133203 

0.00133 

0.00099 

145.555 

185.555 

0.00589893 

0.00642984 

0.00590 

0.00643 

147.555 

183.555 

-0.00016078 

0.00000016 

-0.00016 


153.655 

177.455 

-0.00004351 

-0.00004743 


-0.00005 

155.455 

175.655 

-0.00016667 

-0.00018171 

-0.00017 

-0.00018 

155.645 

175.465 

-0.00008211 

0.00010991 


-0.00011 

155.655 

175.455 

-0.00096958 

0.0 

-0.00097 


155.665 

175.445 

-0.00008292 

-0.00011138 


-0.00011 

157.455 

173.655 

-0.00018515 

0.0 

-0.00018 


162.556 

168.554 

0.00016140 

0.00017588 

0.00016 

0.00018 

163.555 

167.555 

0.00275303 

0.00300082 

0.00275 

0.00300 

164.556 

166.554 

-0.00013869 

0.0 

-0.00014 


165.545 

165.565 

-0.00100294 

0.00134642 

-0.00100 

-0.00135 

165.555 


-0.00868237 

0.0 

-0.00868 
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9. COMPUTATIONAL FORMULAS AND APPLICATIONS 

For numerical computations the diurnal terms in Equations (7.18), (7.29), 
(7.31) and (7.39) to (7.44) for m, > H/C D, and the Euler angle perturbations 
are written in terms of the common multipliers of the lunar and solar terms 
that arise in the theory of precession and nutation [Woolard, 1953, pp. 124-125] . 
The common multiplier of the solar terms is 


Ke = -3 




(9.1) 


= - 3484 " 15 (Jul i an century) -1 
and the common multiplier of the lunar terms is 


K 


9 




7552 "4295 (Jul ian century) -1 


(9.2) 


From Equation (4.24), 


which is written in 



= - 3 


Gm d 


terms of the constants 



K and K* as 

0 


(9.3) 



(9.4) 


In order to make the value of the secular Love number k , given by (5.17), con- 
sistent with the values (9.1) and (9.2) of K and K , k is written as 

© s 



The computational form of Equation (7.18) for the polar motion is 


(9.5) 
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VC- 

m = e + v 


C-A / 

n\ 



(9.6) 


where (A s /AH 2 ) and k s are computed from (9.4) and (9.5) respectively. The 
direction cosines of the axis of figure are computed using Equation (7.29) in the 
form 




The computational form of Equation (7.31) for the angular momentum vector is 
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The computational forms of Equations (7.39) to (7.44) for the Euler angle per- 
turbations are obtained similarly in terms of the dimensionless ratio (Aj/AO 2 ). 

Equations (9.6), (9.7) and (9.8) are used to compute the coefficients listed in 
Table 9.1 for the diurnal motion of m, , and H/C ft in a rigid Earth. Table 9.2 
gives coefficients for evaluating the diurnal motion of m, 4> { , and H/CQ in a de- 
formable Earth with k = 0.29. The computational forms of Equations (7.39) to 
(7.44) are used to obtain the coefficients given in Tables 9.3 and 9.4 for the Euler 
angle perturbations. Table 9.3 is for a rigid Earth and Table 9.4 is for a de- 
formable Earth with k = 0.29. 


The effect of polar motion on latitude and time is shown in Figure 9.1. The 
true equinox of date T t corresponds to the ascending node of the true equator of 
date on the mean ecliptic of date. The true sidereal system has its x axis directed 
toward T' t , its y axis 90° eastward frorn'T T in the true equator of date and its z 
axis along the rotation axis‘s . The transformation from the true sidereal system 
into the terrestrial or "observatory-fixed" system x, y, z is given by 


Serrestrial ^2 ^ — m i ) R l( m 2 ) R 3 (GASTI) x t 


rue siderial 


(9.9) 


where R. (a.) denotes the rotation of a coordinate system about its j axis through 
the angle a, and GAST1 is the Greenwich apparent sidereal time corrected for 
polar motion. Let x T , y T , and z T correspond to the system obtained by rotating 
the true sidereal system through GAST1 about a > . Then 


^terrestrial = R 2 (" m l) R 1 ( m 2 ) *T (9 - 10) 

Latitude and longitude are denoted respectively by <t> T and A T in the x T , y T , z T sys- 
tem and by $ a A in the terrestrial system. Equation (9.10) is used to obtain 

cos $ T cos A T - cos $ cos A - sin <I> (9*11) 

cos $ T sin Aj. = cos $ sin A - m 2 sin $ (9*12) 

sin c|> = sin $ + cos ^>(m 1 cos A + m 2 sin A) (9.13) 

in which second and higher order terms in m x and m 2 are neglected. From 
(9.13), the first order expression for the latitude variation is 


= m 1 cos A + m 2 sin A (9.14) 

Equations (9.11) and (9.12) are combined to give 

A T -A=tan<I>(m 1 sin A - m 2 cos A) (9.15) 
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Table 9-1 


Coefficients for the Diurnal Motion of the Rotation Axis, the Axis of Figure 
and the Angular Momentum Vector in a Rigid Earth 


COF FF I C I F NT S IN SFCQMDS np ARC 
SIMP FDR x-Cn^pnMFMTS 
COS INF F flR Y-CnMPHNFNTS 


INDEX 

TIDAL ARC* 

cnne number 

OIST. 

BODY 

ph r m 

CHEFF [CI ENTS 
LLP F 

DF 

n 

fix 

RDTAT ION 
AXIS 

AXIS OF 
F I CHRP 

ANOilLAR 

MnMFNTUM 

I 

105.955 

M 

1 

— 4 

0 

-2 

0 

-2 

-0.0000023 

o.o 

-0.0000023 

2 

107.755 

M 

1 

-2 

0 

-2 

-2 

-2 

-0.00000 06 

0.0 

-0.0000096 

3 

109.555 

M 

1 

U 

0 

-2 

-4 

-2 

-0*000005 ft 

0.0 

-0.0000058 

4 

115.845 

M 

1 

-3 

0 

-2 

0 

-I 

-O. 00000 42 

O.U 

-0.0000042 

5 

115.855 

M 

1 

-3 

0 

-2 

0 

-2 

-0.0000 21 6 

0.0 

-0.00002 16 

A 

117.645 

M 

L 

-1 

0 

-2 

-2 

-1 

— 0 ,0000106 

o.o 

-0*0000106 

7 

117.655 

M 

1 

-1 

0 

-2 

-2 

-2 

—0*00005*4 

o.o 

-0.0000662 

8 

lift. 654 

M 

1 

-1 

1 

-2 

-2 

-2 

-o.o oo no *2 

0 , 0 

-0.0000042 

9 

119.445 

M 

1 

1 

0 

-2. 

-4 

-1 

-0.0000020 

0.0 

-0.0000020 

Ll> 

119.466 

M 

1 

1 

0 

-2 

-4 

-2 

-0.0000107 

0.0 

-0.0000107 

LI 

124.756 

M 

1 

-2 

-1 

-2 

u 

-2 

0.0000025 

0.0 

0.0000026 

L2 

125.745 

M 

1 


0 

-2 

0 

-1 

-0.0000345 

0.0 

-0.0000344 

13 

125.755 

M 

1 

-2 

0 

-2 

0 

-2 

-0.0001 R3 3 

0.0 

-O.OOCHR27 

1* 

126.556 

M 

1 

0 

-1 

-2 

-2 

-2 

o.ononoM 

0.0 

0.0 000031. 

15 

126.655 

M 

1 

-1 

0 

-2 

-1 

-2 

0.0000021 

0 . 0 

(1 .0000021 

16 

126.754 

M 

1 

-2 

L 

-2 

a 

-2 

-0*0000029 

o . n 

-0,0000020 

17 

127.545 

H 

1 

0 

0 

-2 

-2 

-1 

-0.00 00 41. fa 

0.0 

-0.0000415 

1« 

127.555 

M 

1 

U 

0 

-2 

-2 

-2 

-0 . 000 2200 

0.0 

-0.0002193 

19 

1 2 H. 544 

9 

1 

0 

1 

-2 

-2 

-1 

-0.0000027 

0 . 0 

-0.0000027 

20 

12«. 550 

M 

1 

0 

5 

— 6 

2 

-6 

-0.000O16U 

0.0 

-0 . 0000 1 50 

2 L 

129.355 

N 

1 

2 

0 

-2 

-4 

-2 

-0,0000066 

0,0 

-0.0000066 

22 

] 33 ♦ H 64 

M 

1 

-3 

0 

-2 

2 

-2 

0 .0000043 

0.0 

0.000004? 

23 

1 34.656 

M 

1 

- L 

'1 

-2 

(J 

-2 

n.mmn n 3 

n . o 

0.0000112 

24 

135.435 

M 

1 

1 

u 

— 4 

0 

-2 

0.00000*2 

0.0 

0.O0000M 

25 

135.635 

M 

1 

-L 

u 

-2 

0 

u 

0.0000077 

o.o 

0.0000077 

26 

135.645 

n 

1 

-l 

u 

-2 

0 

-1 

-0.0002505 

0.0 

-0.0002496 

27 

135.655 

M 

1 

-1 

0 

-2 

0 

-2 

-0.0013287 

0.0 

-0.0013244 

2ft 

135.855 

M 

1 

-3 

0 

0 

0 

U 

0.0000035 

0.0 

0.0000035 

29 

136.456 

M 

1 

1 

-l 

-2 

-2 

-2 

0.0000024 

0.0 

0.0000024 

30 

136.555 

M 

L 

U 

(j 

-2 

-l 

-2 

0.0000072 

0.0 

0.000007 1 

31 

136.644 

M 

1 

-1 

l 

-2 

0 

-1 

-0*0000020 

0*0 

-0.0000020 

32 

136.654 

M 

L 

-1 

1 

-2 

0 

-2 

-0 .0000125 

0.0 

-0.0000124 

33 

137.445 

M 

1 

1 

0 

-2 

-2 

-1 

-0.0000473 

0.0 

-0.0000471 

34 

137.465 

M 

3. 

1 

0 

-2 

-2 

-2 

-o.ono25 n 

0.0 

-0.0002503 

35 

137.655 

M 

1 

-1 

0 

0 

-2 

0 

n.ooooi43 

0.0 

0.0000142 

36 

137.665 

M 

1 

-1 

0 

0 

-2 

-1 

-0 .0000044 

0.0 

-0.0000044 

37 

138.444 

M 

1 

1 

1 

-2 

-2 

-L 

-0.0000020 

0.0 

-0. ncioon?o 

38 

138.454 

M 

1 

L 

1 

- 2 

-2 

-2 

-0.0000117 

0.0 

-0.0000116 

39 

1 39.455 

M 

1 

1 

0 

0 

-4 

0 

0.0000025 

0.0 

0.0000026 

40 

143.535 

M 

t 

0 

0 

— 4 

2 

-2 

0.00000 30 

0.0 

0.00000 30 

41 

143.745 

M 

1 

-2 

0 

-2 

2 

-1 

0.0000036 

0.0 

0.0000035 

42 

143.755 

M 

1 

-2 

0 

-2 

2 

-2 

0.0000201 

0.0 

0.0000200 

43 . 

1 44.546 

M 

1 

0 

-1 

-2 

0 

-1 

0*0000027 

0.0 

0.0000027 

44 

144.556 

M 

1 

0 

-1 

-? 

0 

-2 

0.0000231 

0.0 

0.0000230 

45 

L45.535 

M 

1 

0 

0 

-2 

0 

0 

0.0000386 

0*0 

0.0000 3$ 5 

46 

145. 545 

M 

L 

0 

0 

-2 

0 

-1 

-0.0012576 

0.0 

-0.0012535 

47 

145.555 

M 

1 

0 

0 

-2 

0 

-2 

-0.00 66 700 

0.0 

-0.0066481 

4ft 

145.755 

M 

1 

-2 

0 

0 

0 

0 

0.0000430 

o.o 

0.0000 42 ft 

49 

145.765 

M 

1 

-2 

0 

0 

0 

-1 

0.0000071 

o.o 

0.0000070 

50 

146.544 

M 

1 

0 

1 

-2 

0 

-1 

-0*00000 21 

0.0 

-0.0000021 

51 

146,554 

M 

1 

0 

1 

-2 

0 

-2 

-0.0000 20 3 

(J . 0 

-0.0000202 

52 

147.355 

P-1 

1 

2 

0 

-2 

-2 

-2 

0 .00000 "*7 

n.o 

0,0000037 

63 

147 . 5*5 

M 

1 

U 

0 

0 

-2 

1 

-0.0000025 

0.0 

-0.0000026 

54 

147.555 

M 

1 

0 

0 

0 

-2 

0 

0.0000864 

0.0 

(!. 00008 6 L 

55 

147.565 

M 

l 

0 

0 

• 0 

-2 

-1 

-0.0000 18 8 

0.0 

-0.O00018P 

56 

148.554 

M 

1 

0 

1 

0 

-2 

0 

0,0000068 

o.o 

0.0000068 

57 

1 52.656 

M 

1 

-1 

-1 

-2 

2 

-2 

0.0000024 

0.0 

0. <30000 2 4 

5H 

153.645 

M 

1 

-1 

0 

-2 

2 

-1 

0.000010ft 

0 . n 

0.00001 OR 

59 

153.665 

M 

1 

-L 

0 

-2 

2 

-2 

0.0000476 

0.0 

0.0000474 

60 

154. 656 

M 

1 

-1 

-1 

0 

0 

0 

-n .00000 26 

0.0 

-0.0000026 

61 

155.435 

M 

l 

1 

0 

-2 

0 

0 

-0,0000029 

0.0 

-0.0000029 

62 

155.445 

M 

1 

1 

0 

-2 

0 

-1 

0.0000336 

o . o 

0.0000336 

63 

155.455 

M 

1 

1 

0 

-? 

0 

-2 

0.0001814 

o.n 

0.000180ft 
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Table 9-1- (continued) 


CHFFFICIFNTS SFCO^OS OF ARC 
SINE FOR X-CnMPOMFNTS 
CHS INE FDR Y-CHMpnWFMTS 


JLND.EX... 

TIDAL ARG.« . 

COO E number 

„ D 1 S T. . 
BODY 

PHIM 

CQE.FF LCIENTS. 
L LP F 

. OF 
D 

OM 

ROTATION 

AXIS 

AXIS OF 

F IGURF 

ANGULAR 

momentum 

64 

155.645 

M 

1 

-1 

0 

0 

0 

1 

-0. 00001*5 

o.n 

-0.0000144 

65 

155.655 

M 

1 

-l 

0 

0 

0 

0 

0.00050*6 

0.0 

0.0005029 

66 

155.665 

M 

1 

-1 

0 

0 

0 

-1 

0.0001011 

0.0 

0 .0001 OOP 

67 

155.675 

M 

1 

-1 

0 

0 

0 

-2 

-0. 00 00029 

0.0 

-0.0000029 

6ft 

156.555 

M 

1 

0 

0 

0 

-1 

0 

-0.0000027 

(1.0 

-0.0000027 

69 

156.654 

M 

1 

-1 

1 

0 

0 

0 

0.0000031 

0 . 0 

o.nnoon 3 n 

70 

L57.445 

M 

l 

1 

0 

0 

-2 

1 

-n .0000027 

0.0 

-0.0000027 

71 

157.455 

M 

1 

1 

0 

0 

-2 

0 

0 . 0000 95 9 

n.o 

0.0000956 

72 

157.465 

M 

1 

1 

0 

n 

-2 

-1 

0 • 00002 3 0 

0.0 

0.0000209 

73 

158.454 

M 

1 

1 

1 

0 

-2 

0 

0.0000041 

0.0 

0.0000040 

74 

161 . 557 

5 

1 

0 

-2 

-2 

2 

-2 

-0.0000070 

0.0 

-0.0000069 

75 

162.556 

S 

1 

0 

-1 

-2 

2 

-2 

-0.0001 706 

0.0 

-0.0001700 

76 

163.535 

M 

L 

0 

0 

-2 

2 

0 

-0 .0000023 

0.0 

-0.00000?3 

77 

163.545 

M 

1 

0 

0 

-2 

2 

-1 

0.0000 328 

0.0 

0.0000327 

7 B 

163.555 

M 

1 

0 

0 

-2 

2 

-2 

-0 .0000049 

0.0 

-0.0000049 

79 

163.555 

S 

1 

0 

0 

-2 

2 

- 2 

-0.0029022 

0.0 

-0.0028927 

80 

163.557 

$ 

1 

0 

-2 

0 

0 

0 

0.0000018 

0 . 0 

o.nooooi p 

81 

163.755 

M 

1 

-2 

0 

0 

2 

0 

0.0000043 

c 

o 

0 . 00000*3 

82 

164. 554 

S 

1 

0 

1 

-2 

2 

-2 

0.00002^2 

0.0 

0.00(11124] 

85 

.164.556 

S 

1, 

0 


0 

0 

0 

0.0000 698 

0.0 

0.0000693 

84 

165, 545 

M 

1 

0 

0 

0 

0 

1 

-0 • 000 1723 

0.0 

-0.000 1'7 IP 

85 

165.555 

5 

1 

0 

0 

0 

0 

0 

0.0027652 

0.0 

0.0027561 

86 

165.555 

M 

1 

0 

0 

0 

0 

0 

0 . 00 5 94 57 

0.0 

b.00592 62 

87 

165.565 

M 

1 

0 

0 

0 

0 

-1 

0.0011 78 4 

0.0 

0,0011743 

88 

165. 575 

M 

1 

0 

0 

0 

0 

-2 

-0.0000253 

0.0 

-0.0000252 

89 

166.554 

S 

1 

0 

1 

0 

0. 

0 

0.0000694 

0.0 

0.0000692 

90 

167.355 

M 

1 

2 

0 

' 0 

-2 

0 

0.0000042 

0.0 

0.0000042 

9 L 

167.553 

S 

1 

0 

2 

0 

0 

0 

0.000001 8 

0.0 

o.nouooiP 

92 

167.555 

S 

1 

0 

0 

2 

-2 

2 

0.0001236 

0.0 

0.0001232 

93 

167.565 

M 

1 

0 

o- 

2 

-2 

l 

-0.0000047 

0.0 

“0.0000047 

94 

167.575 

M 

1 

0 

0 

2 

-2 

0 

-0,0000023 

0.0 

-0. on 00 02 3 

95 

168.554 

5 

1 

0 

1 

2 

-2 

2 

0.0000072 

0.0 

0.0000072 

96 

172.656 

H 

1 

-1 

-1 

0 

2 

0 

o. 0000038 

0.0 

n.0000038 

97 

173.445 

M 

1 

1 

0 

-2 

2 

-1 

0.0000027 

0.0 

0.0000027 

98 

173.645 

M 

1 

-L 

0 

0 

2 

1 

-0.000 no 29 

0 . 0 

-0. 0000029 

99 

173.655 

M 

1 

-l 

0 

0 

2 

0 

0.0000 901 

o.o 

o.oooopqp 

LOO 

173.665 

M 

1 

-L 

0 

0 

2 

-1 

0.0000178 

0.0 

0 • 0000178 

1QL 

174.456 

M 

1 

1 

-1 

0 

0 

0 

0.0000029 

0.0 

U.OOOU02P 

102 

174.555 

M 

1 

0 

0 

0 

1 

0 

-0. 000 00 25 

0.0 

-0.0000025 

103 

175.445 

M 

1 

1 

0 

0 

0 

l 

-0.0000138 

0.0 

-0.0000137 

104 

L75.455 

M 

1 

i 

0 

0 

0 

0 

0.0004694 

0.0 

0.0004679 

105 

175.465 

M 

l 

1 

0 

0 

0 

-L 

0.00 00 930 

0.0 

0.0000927 

106 

175.475 

M 

i 

1 

0 

0 

0 

- 2 

-0.0000021 

0.0 

-0.0000021 

107 

175.655 

M 

1 

-1 

0 

2 

0 

2 

-0.0000073 

0.0 

-O.OUOUC73 

1D8 

175.665 

M 

L 

-l 

0 

2 

0 

1 

-0 .0000046 

(1 . 0 

-0.0000046 

109 

L75.675 

M 

1 

-1 

0 

2 

0 

0 

-0.0000027 

0.0 

-0.0000027 

110 

176.454 

M 

l 

L 

1 

0 

0 

0 

-0. 00000 24 

0.0 

-0.0000024 

L 11 

177.455 

H 

1 

1 

0 

2 

“2 

2 

-0.0000019 

0.0 

-0*0000019 

1 L 2 

182. 556 

M 

1 

0 

-1 

0 

2 

0 

0.0000049 

0.0 

0.0000049 

113 

183.545 

M 

1 

0 

0 

n 

2 

1 

0.0000025 

0.0 

0.0000025 

114 

183.555 

M 

1 

0 

0 

0 

2 

0 

0.0000756 

0.0 

0.0000764 

115 

183.565 

M 

1 

0 

0 

0 

2 

-1 

0.0000148 

0.0 

0.0000147 

116 

1B5.355 

M 

1 

2 

0 

0 

0 

0 

0.0000367 

0.0 

0.0000366 

117 

185.365 

M 

1 

2 

0 

0 

0 

-1 

0.0000073 

0.0 

0.0000073 

118 

185.555 

M 

1 

0 

0 

2 

0 

2 

0.0002483 

0.0 

0.0002475 

119 

185.565 

M 

1 

0 

0 

2 

0 

1 

0.0001589 

0.0 

0.0001584 

120 

185.575 

M 

1 

0 

0 

2 

0 

0 

0.0000333 

0.0 

0.000033? 

121 

185.585 

M 

.1 

0 

0 

2 

0 

-1 

0.0000021 

0.0 

0.0000021 

122 

191.655 

M 

1 

-1 

0 

0 

4 

0 

0.0000022 

0 . o 

0.0000022 

123 

193.455 

M 

1 

1 

0 

0 

2 

0 

0.0000 116 

0.0 

0.0000116 

124 

193.465 

M 

1 

1 

0 

0 

2 

-1 

0.0000022 

0,0 

0.000002? 

125 

193.655 

M 

1 

-1 

0 

2 

2 

2 

0.0000088 

0.0 

0.0000087 

126 

193.665 

M 

1 

-1 

0 

2 

2 

l 

0.0000056 

0.0 

0.0000066 

127 

195.255 

M 

1 

3 

0 

0 

0 

0 

0.0000028 

0.0 

0.0000023 

128 

195.455 

M 

1 

l 

0 

2 

0 

2 

n. 0000460 

0.0 

0,0000459 

129 

195.465 

M 

1 

1 

u 

2 

0 

1 

0.0000294 

0.0 

0.0000293 

130 

195.475 

M 

1 

1 

0 

2 

0 

0 

0.0000062 

o.o 

0 . 0000062 

1 31 

1X3.555 

M 

L 

0 

0 

2 

2 

2 

0.0000072 

0.(1 

0.0000072 

132 

LX3.565 

M 

1 

0 

0 

2 

2 

1 

0.0000046 

0.0 

0.0000046 
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Table 9-l-(continued) 

C DF FF I C I f NT S I N SFCPMOS nF ARC 
S1NF Fnp X-CnMpPMFNTS 
CDS INF FOR Y -COMPONENTS 


INDEX. 

TIDAL ARG. 

C 0D F NUMBER 

DIS7. 

RPDY 

PHIM 

COEFFICIENTS 
L LP f 

HF. 

n 

nn 

ROTATION 

AXIS 

AXIS OF 
F 1FUPF 

ANGULAR 

momentum 

133 

1X5-355 

M 

L 

2 

U 

2 

0 

2 

0 .0000059 

0.0 

0.0000059 

134 

LX5.3A5 

M 

1 

2 

u , 

2 

0 

1 

0*0000039 

0.0 

0.0000 039 

135 

1E3-455 

M 

1 

L 

0 

2 

2 

2 

0.00 0001 7 

0.0 

0.0000017 


Explanation of symbols 

PHIM is the Greenwich mean sitierial time; ^ 

L, LP, F, D and OM are Brown’s fundamental 
arguments; l, V 9 F, D and ^ 

Constants 

k " 0 

K, - -71*552430 x IQ 5 Julian century- 1 
K m = -31*484150 x 10 J Julian century" 1 
<C-A)/C = 3.272930 x 10" 1 
n = 380!9856 day” 1 
(C^/a,) = 2.343852 x 10 < 

J 2 = 1.082645 x 10-* 

(Oi 9 /a±) - 3.334320 x 10 # 
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Table 9-2 


Coefficients for the Diurnal Motion of the Rotation Axis, the Axis of Figure 
and the Angular Momentum Vector in a Deformable Earth 

COfrFP ICIFIMTS TM SECONDS DF ARC 
s r mf for x-cnMpnMFMTs 
COSINE FDR y-compo mfnts 


INDEX TIDAL ARC. DIST. COEFFICIENTS OF 

CODE NUMBER BODY PHIM L IP F I) 


1 

10 5*955 

M 

1 

-4 

0 

-2 

0 

2 

107*755 

M 

1 

-2 

0 

-2 

-2 

3 

109.555 

M 

1 

0 

0 

-2 

-4 

4 

115.845 


1 

-3 

0 

. - 2 

0 

5 

115.855 

M 

1 

-3 

Li 

-2 

0 

6 

117.645 

M 

1 

-1 

0 

-2 

-2 

7 

117.655 

M 

1 

-1 

0 

~2 

-2 

8 

1 1 8 * 6 54 

M 

1 

-1 

1 

-2 

-2 

9 

119.445 

M 

1 

1 

0 

-2 

— 4 

10 

119.455 

M 

1 

1 

0 

-?■ 

-4 

11 

124.756 

M 

1 

-2 

-1 

-2 

o 

12 

125.745 

M 

l 

-2 

0 

-2 

0 

13 

125.755 

M 

1 

-2 

0 

-2 

u 

14 

126.556 

M 

1 

0 

-1 

-2 

-2 

15 

126.655 

M 

1 

-1 

0 

'2 

-1 

IB 

126.754 

M 

1 

-2 

1 

-2 

0 

17 

127.545 

M 

1 

0 

u 

-2 

-2 

IB 

127.555 

M 

1 

0 

0 

-2 

-2 

L 9 

128.544 

M 

1 

0 

1 

-2 

-2 

20 

128*550 

M 

1 

(J 

5 

-6 

2 

21 

129.355 

M 

1 

2 

0 

-2 

-4 

22 

133.855 

M 

1 

-3 

0 

-2 

2 

23 

134.656 

M 

1 

-1 

-1 

-2 

o 

24 

135*436 

M 

1 

1 

11 


0 

25 

135.635 

M 

1 

-1 

0 

'2 

0 

26 

135.645 

M 

1 

-1 

0 

-2 

0 

27 

135.655 

M 

i 

-1 

0 

‘ -2 

0 

2 R 

135.855 

M 

1 

-3 

0 

0 

0 

29 

136.456 

M 

1 

1 

-1 

“2 

-2 

30 

136.555 

M 

1 

0 

0 

- 2 

-1 

31 

136. 64 4 

M 

1 

-1 

1 

-2 

0 

32 

136.654 

M 

1 

-1 

1 

- ? 

o 

33 

1 37.445 

M 

1 

1 

0 

-? 

-2 

34 

137.455 

M 

1 

1 

0 

-2 

-2 

3 5 

137.655 

M 

1 

-1 

0 

0 

-*2 

36 

137.665 

M 

1 

-1 

0 

r> 

-2 

37 

138.444 

M 

1 

1 

1 

-2 

-2 

38 

138.454 

M 

1 

1 

1 

-2 

-2 

39 

1 39.455 

M 

1 

1 

0 

0 

-4 

4 fl 

143.535 

M 

1 

n 

0 

-4 

2 

41 

143.745 

M 

1 

-2 

0 

-2 

2 

42 

143,755 

M 

1 

-2 

0 

-2 

2 

_43 

. 14 . 4,546 

M 

1 

0 

-1 

-2 

u 

44 

144.556 

M 

1 

0 

-1 

-2 

0 

45 

145,535 

M 

1 

0 

u 

-2 

o 

46 

145.545 

M 

1 

0 

0 

-2 

0 

0 

47 

145.565 

M 

1 

0 

0 

-2 

48 

145.755 

M 

1 

-2 

0 

0 

o 

49 

145.765 

M 

1 

-2 

0 

0 

o 

50 

146.544 

M 

1 

0 

1 

-2 

o 

51 

146.554 

M 

1 

0 

1 

-2 

o 

52 

147.355 

M 

1 

2 

0 

-2 

-2 

53 

147.545 

M 

1 

0 

0 

0 

-2 

54 

147.555 

M 

1 

0 

0 

0 

-2 

55 

147.565 

M 

1 

0 

0 

0 

-2 

56 

148.554 

M 

1 

0 

1 

0 

-2 

57 

152-656 

M 

1 

-1 

-1 

-2 

2 

5« 

153.645 

M 

1 

-1 

0 

-2 

2 

59 

153.655 

M 

1 

-1 

0 

-2 

2 

60 

154, 656 

M 

1 

-1 

-1 

0 

0 

61 _ 

155.435 

M 

1 

l 

0 

-2 

0 

62 

155.445 

M 

1 

1 

0 

-2 

0 

63 

155.455 

M 

1 

1 

0 

'2 

0 


ON 

ROTATION 

AXIS 

AXIS OF 
F I Oil R F 

angular 
muwfntmm , 

-2 

-ft .0000022 

0.0001699 

-0. O0000 16 

-2 

-0.00000 QQ 

0.0007 105 

-0.0000066 

-2 

-0.0000054 

0.0004325 

-0 . 0000040 

~ 1 

-n .ononoM) 

0.0003244 

— 0 . 0000029 

-2 

- 0. 00 00 20 4 

0.00 166 R4 

-0.0000149 

'1 

-0*00001 DO 

0 . 0008 188 

-o. 0000073 

-2 

-o. 0000 524 

0.0042946 

-0.0(10038 L 

-2 

-0.0000039 

0 • 0003244 

-0.no 00020 

- 1 

-0*0000019 

0.0001545 

-0.0000014 

-2 

-n. ooon iru 

0*0008342 

-o . 0000074 

-2 

0.0000024 

-0.0(102008 

0.0000017 

- 1 

-n .00003*0 

0 . 002 7 8 10 

-0 . 0000238 

-2 

-o, 0001 75 L 

0*0147545 

-0.0(i01262 

-2 

0.0000029 

-0.0002472 

O.oo no o?l 

-2 

0.0000020 

-0.0001699 

U. 0000014 

-2 

-0.0000027 

0.0002317 

-n. 0000020 

-l 

-n.oooo3°« 

0.0033681 

-0.0000286 

-2 

-n.0002105 

0. 017 K 1.39 

-0. non 15 1 s 

-1 

v -0 . 0000 02 6 

0,0002163 

-o.noouoi* 

— 6 

v -0,0000 1^4 

0.0012206 

-o.nnooj n3 

-2 

-0.0000064 

0.0005408 

-0.0000046 

-2 

0 .000 Du A 1 

-0 • 00(13 55.4 

0.0000070 

-2 

0.0000 109 

-0.0009425 

0 . 0000078 

-2 

n .ooooo so 

—0 . 0004326 

0. 0000036 

0 

o.oo on 07 s 

-0. 0006490 

0 . 0000053 

-1 

-0.0002420 

0.02J. 014 0 

-0*0001725 

-2 

-0.0012840 

0.11 14976 

-0.00091 49 

0 

n. ooo no *4 

-0,0002936 

0.0000074 

-2 

0.0000023 

-0.0002009 

0.0000016 

-2 

0,00000 69 

-0 • 0006 02 6 

0 . 0000049 

-i 

-0.0000020 

0.0001700 

-0,0000014 

-2 

-0.00001 21 

0 . 0010507 

-0 * 0(10008 6 

- 1 

-0.0000457 

0 . 00 39R6-5 

-n. 0000326 

-2 

-0 .000 24 30 

0.021 1842 

-0.0001729 

0 

0.0000138 

-0.0012052 

0.0000098 

-1 

-0*0000043 

0.0003708 

-0.0000030 

-1 

-n.ooonu) 9 

0.0001700 

-0.0000014 

-2 

-o. 00001 1 3 

0.0009RR9 

-0.0000080 

0 

0*0000025 

-0.0002163 

0. nooooiR 

-2 

0 . OOOOO *0 

-0.0002627 

0.00(10021 

-1 

0.0000035 

-0.000309 1 

0 .000002 5 

-2 

0.0000196 

-0.0017462 

0.0000] 38 

'1 

0.0000026 

-0.0002316 

U. OOOOO 18 

-2 

0.0000225 

-0. 0020 OR 9 

0.0000189 

0 

0.0000377 

-0.0033687 

0.0000266 

-1 

-0.0012292 

0.1097934 

-0.00086 39 

-2 

-0.0065199 

0.5824076 

-0.0045924 

0 

0.00004 20 

-0.0037551 

0 . 0000296 

-1 

0.0000069 

-0.0006 18 1 

0.000 0049 

-1 

-0.0000021 

0.0001R54 

-0.0000015 

-2 

-0.0000199 

0.0017771 

-0.0000140 

-2 

0.00000*6 

-0.0003245 

0. 0000075 

1 

-0. 00 00024 

0.0002163 

’-0.00000 1 7 

0 

n. 0000846 

-0.0075R75 

0.0000595 

-1 

-0.0000 184 

0 . 00 L 65 3 5 

-0. 0000 130 

0 

0.0000057 

-0.0005 TOO 

0. 0000040 

-2 

0.0000024 

-0.0002164 

0.0000017 

-1 

0.0000106 

-0.0009736 

0.0000074 

-2 

0.0000470 

-0.0042963 

0 . 000032 R 

0 

-0.0000025 

0. 0002 3 18 

-0. 0000018 

0 

-0.0000029 

0.0002627 

-G.OC00020 

-1 

0.0 IJOO 3 3? 

-0.0030445 

0.0000231 

-2 

0.0001794 

-C. 0164589 

0.0001249 
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Table 9-2— (continued) 


COEFFICIENTS IN SFCONOS OF APC 
SlWr FOR y-cnwpnwFNTS 
COSINF FOR Y-CnMPOMFNTS 


lnh£x 

TTn&l ARO. 

OI£I~ 

BODY 

PHI 

-COEFF IC I.EN7S 
M L LP F 

OF 

D 

OM 

ROTATION 

AXIS 

AXIS OF 
F IC.URF 

ANGULAR 

j*n»FMT><H 


CODE NUMBER 

64 

155.645 

M 

1 

-1 

0 

0 

0 

1 

-0.0000143 

0.0013136 

-0.0000100 

65 

1 55.655 

M 

L 

-1 

0 

0 

U 

0 

0.0004990 

-0.0458069 

0,0003474 

66 

155.665 

M 

1 

-1 

0 

0 

0 

-1 

o.oooiono 

-0.0091799 

0.0000696 

67 

155.675 

M 

1 

-1 

0 

0 

0 

-2 

-0.0000029 

U. 0002627 

-0.0000020 

68 

156.555 

M 

L 

0 

0 

0 

-1 

0 

-0.0000027 

0.0002473 

-o. oo no oi q 

69 

156.654 

M 

1 

-1 

1 

0 

0 

0 

0.0000030 

-0.0002782 

0.0 000021 

70 

157.445 

M 

1 

1 

0 

n 

-2 

1 

-0.0000027 

0.0002473 

-o.oonnoip 

71 

157.455 

M 

1 

L 

u 

0 

-2 

0 

0,00 00 949 

-0.0087473 

0.0000660 

72 

157.465 

M 

1 

1 

0 

0 

-2 

-1 

0.0000206 

-0.0019164 

0.0000145 

73 

158.454 

M 

1 

1 

1 

0 

-2 

0 

0.0000040 

-0.0003709 

0.0000028 

74 

161 , 557 

S 

1 

0 

-2 

-2 

2 

-2 

-0 .0000069 

0.00064*8 

-0.0000048 

75 

162.556 

s 

1 

0 

-1 

-2 

2 

-2 

-0.0001702 

0.0159356 

-0.0001174 

76 

163.535 

M 

1 

0 

0 

-2 

2 

0 

-0.0000023 

0,0002 164 

-0.0000016 

77 

163.545 

M 

1 

0 

0 

- 2 

2. 

-1 

0.0000 328 

-0.0030757 . 

.0.0000226 

78 

163.555 

M 

1 

0 

0 

-2 

2 

-2 

-0.0000049 

0.0004637 

-0.0000034 

79 

163.555 

s 

1 

0 

0 

-2 

2 

-2 

-0.0028973 

0.2718555 

-0.0019980 

80 

163.557 

s 

1 

0 

-2 

0 

0 

0 

0.0000018 

-0.0001711 

0.0000013 

81 

163.755 

M 

1 

-2 

0 

n 

2 

u 

0.0000043 

-0.0004018 

0.0000030 

82 

164.554 

s 

1 

0 

1 

-2 

2 

-2 

0.0000242 

-0.0022745 

0.0000167 

83 

164.556 

s 

1 

0 


. 0 

.-0 

0 

0.0000697 

-0.0O65526 

0.0000480 

84 

165.545 

M 

1 

0 

0 

0 

0 

1 

-0.0001723 

0.0162285 

-0.0001186 

85 

165.555 

s 

1 

0 

0 

0 

0 

0 

0.0027662 

-0.2604434 

0.0019037 

86 

165.555 

M 

1 

u 

0 

0 

0 

0 

0.0059437 

-0,5600078 

0.0040934 

87 

165.565 

N't 

1 

0 

u 

0 

0 

-1 

0.00 1) 784 

-0.11 10032 

0,0008113 

88 

165.575 

M 

1 

u 

0 

* 0 

0 

-2 

— 0 .0000263 

0.0023802 

-0.0000174 

89 

166.554 

s 

1 

0 

1 

0 

0 

0 

0.0000694 

-0.0065627 

0.0000478 

90 

167.355 

M 

1 

2 

0 

0 

-2 

0 

0.00000 43 

-0.0004019 

0.0000029 

91 

167.553 

s 

1 

0 

2. 

0 

0 

0 

0.000001 8 

-0.0001711 

0.0000012 

92 

167.555 

s 

1 

0 

0 

2 

-2 

2 

0.0001238 

-0.0117079 

0.0000851 

93 

167.565 

M 

1 

0 

0 

2 

-2 

1 

-0.0000047 

0.0004482 

-0.0000033 

94 

167.575 

M 

1 

0 

0 

2 

-2 

0 

-0.00000 23 

0.0002 164 

-0.0000016 

95 

168.554 

s 

1 

0 

1 

2 

-2 

2 

0.0000072 

-0.0006845 

0.0000050 

96 

172.686 

M 

1 

-1 

-1 

0 

2 

U 

0.0000039 

-0.0003710 

0.0000026 

97 

173*445 

M 

1 

1 

0 

-2 

2 

-1 

0.0000027 

-0.0002628 

0. 0000019 

98 

173.645 

H 

1 

-1 

0 

0 

2 

1 

-n .ooooo?9 

0. 00027P2 

-o, 0000020 

99 

173.655 

M 

1 

-1 

0 

0 

2 

0 

0.0000909 

-0.0087485 

0.0000620 

100 

173.665 

M 

1 

-1 

0 

0 

2 

-1 

0.00001 PQ 

-0 . 0017312 

0.0000123 

101 

174.456 

M 

1 

1 

-1 

0 

(J 

0 

0.0000029 

-0.0002782 

0.0000020 

10 2 

174.555 

M 

1 

0 

u 

n 

1 

0 

-0. 0000026 

0.0002473 

-0.0000017 

103 

175.445 

M 

1 

1 

0 

0 

0 

1 

-0.0000 139 

0.0013448 

-0.0000095 

104 

175.455 

M 

1 

L 

0 

0 

0 

0 

0.0004747 

-0.0468 145 

0.0003237 

105 

175.465 

M 

1 

1 

0 

0 

u 

-1 

o.oo on 940 

-0.0090732 

0.0000640 

10 6 

175.475 

M 

1 

L 

0 

0 

0 

-2 

-0.00000 21 

0.0002009 

-0.0000014 

107 

175.655 

M 

1 

-L 

0 

2 

0 

2 

-0.0000074 

0*00071 1U 

-0.0000050 

108 

175.665 

M 

1 

-1 

0 

2 

0 

1 

-0.0000046 

0.0004483 

-0.0000037 

109 

L75.675 

M 

1 

-1 

0 

2 

0 

0 

-0 • 000002 7 

0.0002628 

-0.0000019 

110 

176.484 

M 

1 

1 

1 

0 

0 

0 

-0 .00000 24 

0.0002319 

-0.0000014 

111 

177.455 

M 

1 

1 

0 

2 

-?. 

2 

-0.0000(119 

0.0001855 

-0.0000013 

112 

182.556 

tt 

1 

0 

-1 

o 

2 

0 

0.00000 60 

-0.0004947 

0.0000034 

113 

183.545 

M 

1 

u 

0 

0 

2 

1 

O.0OOO025 

-0.0002473 

0.0000017 

114 

183.555 

M 

1 

0 

0 

0 

2 

0 

0,0000772 

-O. 0076053 

0.0000521 

115 

183.565 

M 

1 

0 

0 

0 

2 

-1 

0.0000 161 

-0,0014840 

0.0000102 

116 

185.355 

M 

1 

2 

0 

0 

0 

0 

0.0000376 

-0.0037099 

0.0000253 

117 

185.365 

rt 

1 

2 

0 

n 

0 

-1 

0.0000075 

-0.0007420 

0.0000051 

118 

185.555 

M 

1 

0 

0 

2 

0 

2 

0.0002539 

-0.0250885 

0.0001709 

119 

1 85.565 

M 

1 

0 

0 

2 . 

0 

1 

0.0001625 

-n .0160610 

0.0001094 

120 

185.575 

M 

1 

0 

0 

2 

0 

0 

0.0000341 

-0.0033699 

0.0000230 

121 

185.585 

M 

1 

0 

0 

2 

0 

-1 

0.0000022 

-0.0002164 

0.0000015 

122 

191. 655 

M 

1 

-1 

0 

0 

4 

0 

0.0000023 

-0.0002319 

0.00000)5 

. 123 

193.455 

M 

1 

1 

0 

0 

2 

0 

0.0000 120 

-0.0012058 

0. QUO 0080 

124 

L 93 * 465 

M 

1 

1 

u 

0 

2 

-1 

0 .00000 23 

-0.0002319 

0,0000015 

125 

193,655 

M 

1 

-1 

0 

2 

2 

2 

0.0000001 

-0.0009121 

0.0000060 

126 

193.665 

M 

1 

-1 

0 

2 

2 

1 

0.0000068 

-0.0005874 

0.0000039 

127 

195.255 

M 

1 

3 

0 

0 

,0 

0 

0.0000029 

-0.0002937 

0.0000019 

128 

195.455 

M 

1 

1 

0 

2 

0 

2 

0.0000476 

-0.0048078 

0.0000317 

129 

195.465 

M 

1 

1 

0 

2 

t> 

1 

0.0000 304 

-0.0030764 

0.0000203 

130 

195.475 

M 

1 

L 

u 

2 

0 

0 

0 . 00000 64 

K). 0006493 

n. 0000043 

131 

1X3.555 

M 

1 

0 

0 

2 

2 

2 

0.0000075 

-0.0007730 

o. non oo. 50 
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Table 9- 2- (continued) 

COFFF ICIFNTS IN SFCQNDS GF ARC 
S1NF FGR X-CnvipDMFMTS 
CDS INF FGP Y-CnNpriNFNTS 

JLMDEX . -T lhAL ARG, DIST. CH6 FF 1C I ENTS flF ROTA T 1 ON AXIS DF ANGULAR 

CODE NUMBER BODY PHIM L LP F 0 OM AXIS FIGURF ^(J^F NT JIM 


132 

1X3.565 

M 

1 

U 

0 

2 

2 

1 

n.ooooo*ft 

-0 .0004347 

0.0000032 

133 

1X5. 355 

M 

1 

2 

0 

2 

0 

2 

0.0000061 

-0.0006339 

0.00000*0 

134 

1X5.365 

M 

1 

2 

0 

2 

0 

1 

0.00000*0 

-0.0004174 

0. 000002 7 

135 

163.455 

M 

1 

I 

0 

2 

2 

2 

0 . 000001 ft 

— 0 . 000 lH 5 5 

0.0000012 


Explanation o 1 symbols 

PHIM is the Greenwich mean side rial time; 

L, LP, F, D and OM are Brown's fundamental 
arguments; ^',F, D and n 

Constants 

k = 0 

K, = -71‘552430 X 10 3 Julian century- 1 
K a = -3T484150 x 10 3 Julian century' 1 
<C-A)/C - 3.272930 x 10“ 3 
II = 36029856 day^ 1 
(Q^) = 2*343852 x 10 A 
J 2 = 1.082645 X 10- A 
(m^/m,) - 3*334320 x 10° 
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Table 9-3 


Coefficients for Perturbations in the Euler Angles of the Rotation Axis, the Axis 
of Figure and the Angular Momentum Vector in a Rigid Earth 

COEFFICIENTS FOR F-ML fr ^ Awr; LF P F(5 Tllft *A T I Hi-S 
1 to i)FCn,\in$ OF ARC 

C.oSTmfs FU» nptTft thfta 

SIWFS FO& nFLTA ftST * SIMIThFTaJ 

ft HT AT I f)M AXIS HF TFKkFbTKTAL 

TIQAL ARGUMENT COEFFICIENTS OF AXIS FIGtlftF t. AXIS 


CODE NUMBERS 

l 

LP 

F 

0 

UM 

CDS 

SIN 

cos 

STM 

LOS 

STu 

' 1U5.955 


4 

0 

2 

u 

2 

-o.o ooo on oi 

-n .ooo uf >ooi 

0. 00 0002 30 

0. 00000230 

0.OOUOU23O 

o . u n nuo23ii 

107.755 


?. 

0 

2 

2 

2 

-o, nonuono^ 

-0 .ODOOOUO 3 

0 . OOU 0095 6 

0,00000955 

11.0001)09 54 

o.onouiwSS 

109.555 


0 

0 

2 

4 

2 

-o.oonuono? 

-n.onounno? 

o.nnuons78 

0. 00000678 

0.00000478 

[' . 1*0)100676 

115. R45 


3 

0 

2 

U 

1 

-0. 00000 00] 

-0.01 >00000 1 

0. 00000419 

O.OOOO041O 

0.00000419 

0 .000004 1 4 

LL5.H55 


3 

0 

2 

0 

2 

-0.0 Of H I OM 07 

-O.OOOOOOU7 

0. noun 2 157 

0. 000021 67 

0 . nOi>U2 147 

n 4 V 

11 7 . BAS 


1 

f) 

2 

2 

1 

-o. nuoon no 3 

-0. 00000003 

O.OOUO 1052 

0.00001052 

0,001)01042 

0,0(11)01 i > 4 7 

117. 655 


L 

t) 

2 

2 

2 

-o.o ooo on 1R 

-o.ooo on oi ft 

0.000055 19 

0.00005519 

0 . 0000551 s 

0.OOOOS5 i'J 

118.654 


1 

-1 

2 

2 

2 

-n. nungnoo i 

-o.oooonoo l 

0 . 000004 L 6 

0.0000041 6 

0. O0UU041 6 

n .000004 1 t 

1 19. AA5 


-1 

0 

2 

4 

l 

-owning on m - 

-o.onnuonni 

0.00000197 

o.o nnoo ),9 7 

0.00000197 

0. 00000197 

119. A55 


-1 

0 

2 

4 

2 

-0. 0U 000 00 3 

-o .or 'uo oo 0 3 

O.OOUO 1066 

o. oooo in 66 

0 , OOOU 1 IJ 66 

0 .00001066 

1 2 A. 756 


2 

1 

2 

0 

2 

0.0 000 OfK.n 

o. oo non ooi 

-n.onoM(i?49 

-0.00000240 

-n. 00000249 

-0. OUOOU/49 

125. 7A5 

1X5.365 

2 

0 

2 

0 

1 

-o. ooooo mo 

-o.oonuooi 3 

O.OOUO 3 05 8 

0. 00003829 

0.000030 

0 .0000 3 62 \> 

125*755 

LX5.355 

2 

0 

2 

U 

2 

-0.00000058 

-o.on nun 062 

0. OOI) 176 8 1 

o. non 1.8863 

0 , nouj. 7 68 1 

0,000 1 8'ih3 

126.556 

4 

0 

1 

2 

2 

2 

o. no o oo no i 

o.ooou one l 

-0 , n0u0f!306 

-0, 0000030 5 

-n , OOO MU 3 05 

-n ,00000305 

126.655 

4 

L 

0 

2 

1 

2 

o.ooou on oi 

0 . oo non oo l 

-0.00000? 10 

—<>♦00000 2 1.0 

-O . 000(10 2 lo 

-0*1)0000210 

126. 75A 


2 

-1 

2 

0 

2 

-n. non on no j 

-0 .0000000 1 

0. 00000266 

0.00000 28 6 

0. 000002 8 6- 

0 ,oo(joo^h r 

127.545 

1 X3.566 

0 

0 

2 

2 

1 

-o.ooou on i? 

-0. on (UK) 01 6 

o.onon^fftft 

0. 00006.404 

n. 0000 3 668 

0 . 0000460^ 

127.555 

1X3.555 

0 

0 

2 

2 

2 

-o. onooon70 

-0 .00000074 

n. 0007121 3 

0.00022648 

0.00(121213 

0 ,uf!02? 1.4 K 

128.544 


0 

-1 

2 

2 

1 

-o.ooooonoi 

-o.onuonooi 

0.nnunn?65 

o.nonoo? 65 

n. 000007 64 

n .o0tio.)7 r"> 

128.550 


0 

-5 

6 

-2 

6 

-o. no ooo oo 5 

-n .uonuonos 

n . 001)01498 

0.00001498 

O. 0 (.11.10 1.496 

0.00001^48 

I 29.355 


-2 

0 

. 2 

4 

2 

-o.ooooooo? 

-o.ooo no on? 

0.00000662 

o.o noon 66? 

0.00000 662 

0. O0o(lf'ro62 

133.855 


3 

f) 

2 

-2 

2 

o. ouoonooi 

o .onuoooo 1 

-O.OOUO 04 24 

-0. 0000042 4 

-0.00000424 

-n . U 0 00047.4 

134.656 


1 

1 

2 

0 

2 

0.00000004 

0.00000004 

-n. ononii 23 

-o.n non u?3 

-O.OOOU 112 3 

-0.000^1173 

135.436 


“I 

1) 

4 

0 

2 

o, no ooo no? 

0.000(1000 2 

-0.00000614 

-0.0000051 4 

— 0 . OOoOU 4 1 4 

-0 ,(100(10 *>14 

) 35.635 

195.^75 

1 

0 

2 

0 

0 

O.OOOU 00 05 

o.oonooooo 

-o.noooi 390 

-0.0000015? 

-u. 0000 1390 

-0,00000 l 42 

135.645 

195.465 

1 

0 

2 

0 

1 

-0 • 0000007? 

- 0.0 not) no 9 2 

0.00022030 

0.00027899 

0.00027030 

0. 00027400 

135.655 

195.455 

1 

0 

2 

0 

2 

-0.0 OOO 04 20 

-o.no non 4 so 

o. no 12 7 R 60 

0.00 1.37 026 

U. 001 27851 

O.0Ol37u?r> 

135.855 

195.255 

3 

0 . 

0 

0 

0 

o. nuoon no 2 

o.unouoouo 

-0. 000(10629 

-0. 00000068 

-0.00000629 

- 0 ,00000068 

136.456 


-1 

1 

2 

2 

2 

o.ooou oooj 

o.oooooooi 

-n.noonn?3« 

-0.00000238 

-0.00000238 

-0 . 00000238 

136.556 


0 

(1 

2 

L 

2 

o. on ooo oo? 

0.00000002 

-0.00000714 

-0. 00000714 

-0. 0(10007 L 4 

-0,00000714 

136.644 


1 

-1 

2 

0 

l 

-0. 00000001 

-O.OOOOOOOI 

0.00000201 

0 • 00000?! u 

0. 00000201 

0 . 00000201 

136.654 


1 

-1 

2 

0 

2 

-0. 00000 004 

-0.00000004 

0.00001244 

0.00001244 

0. UUU01 244 

0. 00001244 

1 37.445 

193.665 

-1 

0 

2 

2 

1 

-0.00000014 

-0.00000017 

0.00o04l48 

0.00004273 

O.O0OU414K 

n . 00006273 

137.455 

193.655 

-1 

0 

2 

2 

2 

-0. OOOU0O79 

-0 .oouo on 8 5 

0. 00U?4 153 

0.0002590} 

0* 000241 63 

0.00024401 

1 37.655 

193.455 

1 

0 

0 

2 

0 

0.00000008 

n. oo ooo ooi 

-0. 00002679 

-0 ■ 000002 67 

-0.00002579 

-0,000002 nl 

137.665 


1 

0 

0 

2 

1 

-o. ooo oo no i 

-0 .0000000 1 

■0. 000 0043 8 

0.00000438 

0 . 00000438 

0.00000438 

L38.444 


-1 

-1 

2 

2 

1 

-o.ooooonoi 

-0.00000001 

0.00000200 

O.OOOU 0200 

0.00000200 

0.00000200 

138.454 


-1 

-1 

2 

2 

2 

-0. 00000004 

-0 .0000000 4 

0. 0000 1165 

0.00001 1 65 

0,00001165 

0 ,01)00 L l 64 

1 39.455 

191.655 

-1 

0 

0 

4 

0 

o.o ooo on 02 

o.ooo uo ooo 

-0.OOO0O477 

-f). OOO OOO 31 

-0.00000477 

-O.UOOOOU31 

143.535 


0 

{1 

4 

-2 

2 

0. 01)000001 

o.oooooooi 

-n. ooo 00302 

-0. OOO0030? 

-0 . OOU 00 3 02 

-0,00000.402 



Table 9- 3— (continued) 


COEFF 1C IFNTS FOR FlfL ER ANGLE P ER TURRATI OlM$ 
im seconds np arc 


C.f IS 1 MF S FOR Op L T A TMFTA 

SINFS FOR DELTA PSI * ST N( THETA ) 


Oi 

00 


tidal 

CODE 

ARGUMENT 

NUMBERS 

L 

COEFF TCI ENTS 
L P F IJ 

OF 

OM 

ROTATION 

AXIS 

COS SIM 

AXIS OF 
FIGURE 

CHS SI* 1 

'143.745 


2 

0 

2 

-2 

1 

o. ooooooo l 

0.0000000 1 

-0.00000355 

-0.00000398 

143,755 

. 

2 

0 

2 

-2 

2 

O.OOOUOOU7 

0.00000007 

-0. 00002004 

-n.00002004 

1 44 . 546 

. 

0 

1 

2 

0 

1 

o. ooooooo i 

0.0000000 1 

—0 ■ 00000265 

-D. 00O00265 

144.556 

. 

0 

1 

2 

0 

2 

0.00000008 

o.ou ooorn ’8 

-o. nnon23oo 

-0 . 00002300 

L45. 535 

185.575 

0 

0 

2 

(J 

0 

0,00000024 

o. ooooooo? 

-0.0nu07l69 

-0. 000005? 4 

145,545 

185.565 

0 

0 

2 

0 

1 

-0. 00000360 

-0.00000464 

0.00109508 

0.00141185 

145.555 

185.555 

0 

0 

2 

0 

2 

-0. 00002102 

-0 .0000 2264 

0.006^0060 

0.00680950 

145.755 

1 85.355 

2 

0 

0 

U 

0 

0.0O0OD026 

0.00000002 

-0.00007945 

-o.noooo 62? 

145.765 

. 

2 

0 

0 

0 

1 

0.0000000? 

0.(1 000000 2 

-0.00000705 

-0.00000705 

146.544 

. 

0 

-1 

? 

0 

1 

-o.oooooooi ' 

-0.00000001 

0.00000211 

o.onnoo?ii 

146.554 

• 

0 

-1 

2 

0 

2 

-0. 00000007 

— 0 .(J 000 000 7 

0. ooo 0202 3 

0. 00002023 

147.355 

. 

-2 

0 

2 

2 

2 

0.00000001 

0. on ooo ooi 

-0.00000369 

-0.00000369 

147,545 

183.565 

0 

0 

0 

2 

-1 

0. 0000000 4 

-0.0 000 00 06 

-0. 00001 225 

0.000017 16 

147.555 

183.555 

(i 

0 

0 

2 

0 

0.00000053 

0. 00000004 

-0.00016150 

-0 • 0000 107 1 

147.565 

L83.545 

0 

0 

0 

2 

1 

-0. 00000 00 5 

-0.0 000 000 7 

0.00001631 

0.000021 ?] 

148.554 

182.556 

0 

-1 

0 

2 

0 

0.00000004 

0.00000000 

-0. 00001069 

-0.00000085 

152.656 

• 

1 

1 

2 

-2 

2 

o. nonooooi 

0.0 noon oo i 

-0.00000240 

-0. 00000240 

153.645 

• 

1 

0 

2 

-2 

1 

O.OQQOOnOA 

0.00000004 

-0. 00001075 

-0.00001075 

153.655 

177.455 

l 

0 

2 

-2 

2 

o. ooo on oi s 

0.000000 ] 6 

-0.00004555 

-0.00004932 

1 54.656 

1 76.454 

1 

1 

0 

0 

0 

-0.00000002 

-0,00000000 

0.00U00A91 

0.00000019 

155.435 

175.675 

-1 

0 

2 

0 

0 

-0.00000002 

-0 .00000000 

0.00000557 

0.00000021 

155.445 

175. 665 

-1 

0 

2 

u 

1 

0.00000009 

0.00000012 

-0.00002888 

-0,00003803 

155.455 

175.655 

-I 

0 

2 

0 

2 

0.00000057 

0.0 OOO OU62 

-0.00017357 

-0.00018808 

1 55.645 

175.465 

1 

0 

0 

0 

-1 

0.00000026 

-0.00000035 

-0.00007823 

0.00010708 

155.655 

175.455 

1 

0 

0 

u 

0 

0. 00000 31 9 

o. oooo on 1 2 

-0.00097083 

-0.00003502 

155.665 

175.445 

1 

0 

0 

u 

1 

0. 0000 00 29 

0.00000038 

-0.OOQOR7O4 

-0.00011451 

155.675 

. 

1 

0 

0 

0 

2 

-0. OOOOOOO J 

-0 .00000001 

0. 00000288 

0, 0000028* 

156.555 

174.555 

0 

0 

0 

1 

0 

-0. 00000002 

-o.ooonoouo 

0.00000524 

0.0000001 R 

156.654 

174.456 

1 

-1 

0 

0 

0 

0, ooooooo? 

o.onoocujon 

-0. 000 005 8 9 

-0. 0000 0020 

157.445 

173.665 

-1 

0 

0 

2 

-1 

0.00000005 

-o. on ooo 007 

-0.00001506 

0.00002046 

157.455 

173.655 

-1 

0 

0 

2 

0 

0. 00000061 

0.00000002 

-0.000 18533 

-0. 00000579 

157.465 

173.645 

-1 

0 

0 

2 

1 

0.00000006 

0.00000008 

-0,00001*08 

-0.00002370 

158.454 

L72.656 

-1 

-1 

0 

2 

0 

0. 00000 00 3 

0.0 ooooooo 

-0 • 000007 86 

-0. 0000002? 

161.557 

. 

0 . 

2 

2 

-2 

2 

-P.OOOUD.OQ? 

-0.00000002 

0,00000694 

. 0. 0000 0 694 

162.556 

168. 554 

0 

1 

2 

-2 

2 

-0. 00000053 

-0.0 000 00 58 

0.000 16284 

0.00017721 

163., 5.35 

167.575 

0 

0 

2 

-2 

U 

-0.00000002 

-0.00000.000 

0. 00 00.04 58 

0.0000.000 3 

163.545 

167.565 

0 

0 

2 

-2 

1 

0. 00000009 

0.0 OOO 00 12 

-0.00002801 

-0.00003745 

1 63.555 

. 

0 

0 

2 

-2 

2 

-0.00000002 

-0.00000002 

0.00000493 

0.00000493 

163.555 

167.555 

0 

0 

2 

-2 

2 

-0.00000909 

-0 .00000990 

0.00276942 

0.003015R7 

L63.557 

1 67.553 

0 

2 

0 

0 

0 

0,00000001 

0.00000000 

-0.00000362 

-0.00000002 


TERRESTRIAL 
L AXIS 

LUS SIM 


-0.00000355 
-U * 0 (JU020U4 
-0 . 00001)2 65 
-0 . 00002300 
-0*000071 69 
0.0U1O95U9 
0.006400 6ft 
-0*00007945 
-0*00000705 
U • 000002 1 1 
0. 00002023 
-0 .000003 69 
-O.OOOOL225 
-0.0001 6150 
0. 0( i u 01 631 
-f) . 0000 10 69 
-0. 00000240 
-0 • 0OOU1U75 
-0.00004555 
0 . 0000049 1 
0. 000005 57 
-0.00002888 
-0.00017357 
-0.00007823 
-0. 00097(184 
-0 .00008704 
O. 000002 HK 
0.00000524 
-0 . 1)0000589 
-U.0UUU1506 
-0.0001853? 
-0.00001808 
-0.000007ft 6 
0.00000 694 
0.00016284 
0.00000458 
-0 , 0000280 1 
0.00000493 
0.00276945 
-0.00000362 


-0 .1)00110 355 
-0. 00002004 
-{’.000002 66 
-0,00002 3011 
-0.00 001152 4 
0.00141187 
n.un6H966ft 
-0.00000*22 
-0.00000705 
O . 0000021 1 

n. ooo(i2o2 ? 

-0 • U0OUU3 69 
0.0000171 6 
-0.00001071 
0,00002121 
-o , oonootjuft 
-O. 00(1011240 
-0.00001075 
-('.00004932 
0 . OOOOUO l 9 
O.OOOOOU21 
-0.0OUU3tt 113 
-0 ,00018808 
0.00010708 
-0 ,0000 35i] 2 
-0.O0OU451 
0.00000288 
0,0000 Do lft 
— O .00000112 0 
. 0.00002046 
-0.00000574 
-M. 00002 379 
-0.00000 112 2 
0.00000694 
0. 00(11 772 1 
0. GOO OGLU Q3_* 
-0 .00003745 
0 ■ 00000493 
D.U0301 540 
-O.OOOOOU02 



Table 9-3— (continued) 


COEFFICIENTS FOK FUIER ANGLE P FR TU» RaT 1 ONS 
IN SECONDS OF ARC 


COS f NFS Fl |P n e L T A ThFTA 

SINFS FOR PFLTA HSI * SIMIThFTaI 


TIOAL 

CODE 

ARGUMENT 

NUMBERS 

COEFF I C I FNTS 
L LP F t) 

OF 

DM 

ROTATION 

AXIS 

Cns Sjw 

AXIS OF 

f i f;i iw f 

COS sin 

TFkrFsTkI AL 
L axis 

CHS S t 1 

163*755 

1 67.355 

2 

n 

0 

-2 

0 

o.CHjnonno3 

O.OfiOOOOOO 

-0. 000008 5 I 

-o.onnnooo^ 

-11 . OOOOOHS 1 

— o . ijnooo<>!)4 

164.554 

• 

0 

-l 

2 

-2 

2 

0. 000 DO 00 8 

■ 0.0 MOO 00 oh 

-0 *00002414 

-0.00002414 

-n.00l)U24l4 

-r«. 00002414 

164.556 

166.554 

0 

i 

0 

O 

0 

0.0(1000046 

0.00(K)O0(JU 

-0. 00013868 

-0.90000038 

-O.onoi 3 h 69 

— 0 . \ ji'lfiOi )u3H 

165.545 

165.565 

n 

0 

0 

0 

-1 

i.l . 00000 32 9 

-0.00000^42 

-0.00100274 

0.00134 62 6 

-0.00 10027 5 

{'.(<0 134 627 

1 65.555 

* 

0 

0 

0 

O 

0 

o.ounonQ(i5 

0.00000909 

-0.00275609 

-0.00276609 

-0*00275 M 3 

-M . 0027 6 .-ll 3 

165.555 

. 

0 

0 

0 

0 

£J 

o. nonu 1946 

0 .0000 1 94 6 

-0.00592617 

-0. 00692617 

-0.00602624 

-ll.l/O 692o24 

. 

165.575 

0 

0 

(1 

0 

-2 

-o.ouooooor 

o. notion oo* 

0. 000025 18 

-0.00002618 

0 » 1|I.30L12 61 K 

-n ,000025 ) n 

, 

173.445 

l 

0 

-2 

2 

-1 

0. ll 0000 001 

-0.00000001 

-0.O0U0027O 

0. 00000270 . 

-0. 01 u.hKj27o 

(•.00000270 

. 

175.475 

i 

0 

n 

0 

-2 

-o.ouounnoi 

o.oooooooi 

o. 00000205 

- 0.00000206 

<1.001100206 

-O.i.M if 1 1 li 1^06 

, 

165.365 

2 

0 

0 

o 

-L 

0. uoooooo 2 

-n.O no uoo 0 2 

-0. onyn0732 

0.00000732 

-0.00000732 

O .00000732 

. 

185.565 

0 

O 

2 

£) 

-1 

o.ooouonoi 

-0 . uo 000 001 

-0. 000002 1 3 

n. on o no ? i 3 

- 0 . r.noo!)2 \ 3 

0. Of’MM'if' L3 

, 

193.465 

1 

0 

0 

2 

-1 

0. 000 00 001 

-U .0000 000 i 

-0.00000222 

0.00000222 

-o. on 000222 

0.00000222 

• 

1 E3.455 

L 

0 

2 

2 

2 

0.0 000 on 01 

- 0 . OM 000 00 1 

-o.onyno ] 67 

o.o no ooj. 67 

-0 • IK K’kjO 1 47 

0. OOOuo 1 67 


G> 

<D 

Explanation of symbola 

I LP, F p D and OM axe Hr own* a fundamental 
arguments; 2 ,•£ , F> D and »1 

Constants 


k= 0 

K,= -7TS62430 x 10 3 Julian century- 1 
= -3T484150 x 10 3 Julian century* 1 
<C-A)/C = 3*272930 x 10 ' 3 
Cl = 360J9856 day ' 1 
(Cg/a^) = 2.343S52 x 10* 

J 3 = 1.082645 x 10 " 3 
(nig/m^ = 3.334320 x 10® 



Table 9-4 


Coefficients for Perturbations in the Euler Angles of the Rotation Axis, the Axis 
of Figure and the Angular Momentum Vector in a Deformable Earth 

f.flEFFICIFNTS FDR FULFR ANGLF PFP TIIR 8 a T I f !n$ 

IN SECONDS np arc 


CnSINFS Ff)R n FLT A TMFTA 

SI NFS FDR OFIT A PS I * STN(THFTA) 


TIDAL 

CODE 

argument 

MUMBERS 

COEFFICIENTS 
l IP F 0 

OF 

DM 

ROTATION 

AXIS 

COS SIN 

AXIS OF 
F IGHKF 

C ns ST N 

tfkrfstri al 

l AXIS 

CHS SIM 

105*955 


4 

0 

2 

0 

2 

-0.00000056 

-0. 00000056 

0.000171 49 

0.000 17 L 49 

0.00000159 

0.000(8)1 69 

107.755 


2 

0 

2 

2 

2 

-0.00000235 

-0 .00000235 

0. 0007 3.713 

0.000717 13 

0.00000 660 

0 .00000 66 11 

109.555 


0 

0 

2 

4 

2 

-0.0 0000 143 

-0.00000143 

0.00043650 

0.00043 6S0 

0. 00000399 

0.0001 'lM^9 

115.845 


3 

0 

2 

0 

1 

-O.OOOUO107 

-0.00000107 

0. 00032731 

0. 000327 3 1 

0.00000^90 

0.00000290 

115.855 


3 

0 

2 

(J 

2 

-0.00000553 

-0.00000553 

n.nn 16R328 

0.00168328 

0.00001490 

0 .0(8*11) 1.490 

117.645 


1 

0 

2 

2 

1 

-0.00000271 

-0.00000271 

0.00082603 

0. 00082 60 3 

0.00000727 

0 .UOOOU72 7 

117.655 


1 

0 

2 

2 

2 

-0.00001423 

-0.00001423 

0 . 00433 ? 74 

0.00433274 

0.00003813 

n. 00008813 

118.654 


1 

-1 

2 

2 

2 

-0.00000 107 

-0.00000107 

0.00032729 

0.00032729 

0 . 000002 H 7 

0.00000287 

119.445 


-1 

0 

2 

4 

L 

- 0 . 0 000 on 51 

-0 .000000 51 

O.nntJl 5585 

0.0001^585 

0.00000136 

n. 000001 3 6 

119.455 


-1 

0 

2 

4 

2 

-0.00000276 

-0.00000276 

0. 00084158 

0. 00084 158 

0.00000736 

O.ljOLJUfJ 73 6 

124.756 


2 

1 

2 

0 

2 

-0.00000067 

0.00000067 

-l). 00020257 

-0.00020257 

-0.00000172 

-o.ijnnooi 7? 

125.745 

1X5.365 

2 

0 

2 

0 

1 

-0. 00000783 

-0.0000 1059 

0. 00238466 

0.00322484 

0.000021 13 

U.00UU2C45 

125.755 

1X5.355 

2 

0 

2 

0 

2 

-0.00004677 

-0.00005096 

0. fU424?H3 

n.Ol 55L864 

U.0UU12215 

0.1)0013(124 

126.556 


0 

1 

2 

2 

2 

0. 00000082 

0.000000 82 

-0. 000 24 931 

-0.00024931 

-0.00000211 

-( ‘.000002 11 

126.655 


1 

0 

2 

1 

2 

0.00000056 

0.00000056 

-0.00017 ) 40 

—0 .000 171 40 

-O.OOOUOJ45 

-O.OO0OO146 

126.754 


2 

-1 

2 

0 

2 

-0. 00000077 

-0 .0 00000 7 7 

0.0OQ23372 

0* 00023372 

0.000U019K 

0 ,00000 19 1' 

127.545 

1X3.565 

0 

0 

2 

2 

1 

-0.0000095? 

-0*00001279 

0.00289886 

0.0038Q465 

0 .00002548 

(‘.0000318? 

127.555 

1X3.555 

0 

0 

2 

2 

2 

-0. 00005644 

-0 *0000 615 5 

0.017 1P741 

0. 01*74333 

0.00014 65.8 

0.00015 o4 6 

128.544 


0 

-1 

2 

2 

1 

-0.000000 7? 

-0.0000007? 

0.00021814 

0.00021814 

0 • 0000 0 1 8 3 

n . 00000183 

128.550 

♦ 

0 

-5 

6 

-2 

6 

-0. 0000040 4 

-0 .0 0000404 

0. 00123091 

0.00123091 

0.00001035 

o . U Ooo 1 o3 5 

129.355 


-2 

0 

2 

4 

2 

-0.00000179 

-0.00000179 

O.OOU54^33 

0.00054533 

0.001)00457 

n . 00000467 

133.855 


3 

0 

2 

-2 

2 

0.00000118 

0.0000011 H 

-0.001)38831 

-0.00033831 

— 0. (>00002 93 

-0.00000293 

1 34.656 


1 

1 

2 

(J 

2 

0.0000031? 

0.00000312 

-0. 0OU9S029 

-0. 00096029 

-0.00000776 

-0.00000776 

135.435 


-1 

0 

4 

0 

2 

0. 00000143 

0.0 0000143 

-0. 00043619 

-0.00043619 

-0.00000355 

— 0 .00000 355 

135.635 

195.475 

L 

0 

2 

0 

0 

0.00000429 

0. 00 000 00 0 

-0.00130785 

-0. 0000007? 

-0.00000960 

-0,00000105 

135.645 

195.465 

1 

0 

2 

0 

1 

-0.00005940 

-0 .0000 7974 

0.01808976 

0.02428307 

0.00015219 

0. 0001 9273 

135.655 

195.455 

1 

0 

2 

0 

2 

-0.00035324 

-0.0003850? 

0. 10757303 

0.11725199 

0.00088321 

0.00094658 

135.855 

195.255 

3 

0 

0 

0 

0 

n. no ooo 1.9 4 

n.onooonoo 

-0.0006Q165 

-0.00000033 

-0.00000434 

-0.00000047 

136.456 


-1 

1 

2 

2 

2 

0.00000066 

0.00000066 

-0.00020251 

-0.000 20251 

-0.0 0000 1 64 

-0. Of JO 001 64 

136*555 


0 

0 

2 

l 

2 

0. 00000 199 

0.00000199 

-0.00060754 

-0.00060754 

-0.00000493 

-0,00000493 

136.644 


1 

-1 

2 

0 

1 

-0.00000056 

-0. 000000 56 

0.000171 36 

0.00017136 

0.00000139 

0.00000139 

136.654 


1 

-1 

2 

0 

2 

-0. 00000348 

-0 .00000348 

0*00105930 

0.00105930 

0.00000860 

0.1)0000 6 60 

137.445 

193.665 

-1 

0 

2 

2 

1 

-0.0000 1126 

-0.00001514 

0.00342774 

0.00461040 

0.00002865 

0.00003643 

137.455 

193.655 

-1 

0 

2 

2 

2 

-0. 0000671 2 

-0 .0000 7314 

0.02043900 

0.02227 5? 4 

0.00016686 

0.00017892 

137.655 

193.455 

1 

0 

0 

2 

0 

0.00000798 

0.00000000 

-0.OO242885 

-0.00000127 

-O.U00017H1 

-0.00000184 

137.663 


1 

0 

0 

2 

1 

-0.00000 12 3 

-0.00000123 

0.00037386 

0.00C137 38 6 

U.0OUU03Q2 

0.00 00030 2 

1 38.444 


-1 

-1 

2 

2 

1 

-0. 00000056 

-0.00000056 

0.000171 35 

0.00017135 

0.00000138 

0.00000138 

138.454 


-1 

-1 

2 

2 

2 

-0. 00000 327 

-0 .0 0000 32 7 

0.00099696 

0.00099696 

0.00000805 

0.00000605 

139.455 

191.655 

-1 

0 

0 

4 

0 

0.000001 48 

-0.00000005 

-0.00045151 

0.00001 534 

-0.00000330 

-O.OOOOOU21 

143.535 


0 

0 

4 

-2 

2 

0. 00000087 

0.00000087 

-0.00026478 

-0.00026478 

-0.00000208 

-0,00000208 



Table 9-4— (continued) 


COEFFICIENTS FOR FlILEft ANGIF PERTURPATI'Ji'iS 
[M SFCflMOS np ABC 

COS IMPS FOR n^LTA THFTA 

SINFS FDR nFLTA RS] * StWTHFTA) 


ROTATION AXIS nF TFKRFSTRtAL 


TIDAL ARGUMENT 
CODE NUMBERS 

CDEFF I Cl FNTS UF 
L LP F L) UM 

AXIS 

CCS 

SIM 

F I Gllk F 

cns 

STM 

L AXIS 

cns 

S IN 

143*745 


2 

0 

2 

-2 

1 

0.00000 102 

0.00000102 

-0.00031151 

-0.00031 151 

-0.00000245 

-0 .0(100(1245 

143-755 


2 

0 

2 

-2 

2 

0.00000578 

0.00000578 

-0.00 176001 

-0.00176001 

-0.0000 1 364 

-0 . 00001384 

144.546 

. 

0 

1 

2 

0 

1 

0. 00000077 

0. 0000 00 7 7 

-0.00023363 

-0.0nn?3363 

-0.00000183 

-(1.(1(1000 i 8 3 

144-556 

. 

0 

1 

2 

0 

2 

0.00000665 

0.00000665 

-0.00202476 

-0. 00?o?47 6 

-0,00001589 

-0 . 0001)1.5 89 

145.535 

L85. 575 

0 

0 

2 

0 

U 

0. 00002229 

0.00000001 

-0 . 0067 14 

-0.00000250 

-0. (5000495? 

-n .00000 362 

145.545 

185.565 

Q 

0 

2 

U 

L 

-0. 00031027 

-0.00041647 

0 . 0944 ft R 86 

0. 126829 62 

. 0.00U75648 

0.0009/527 

145-555 

185.555 

a 

0 

2 

0 

2 

-0.001 84457 

-0.0020 10 46 

0.56174046 

0,61225933 

0.00442146 

0.00 4 7 632 7 

145-755 

185.355 

2 

0 

0 

U 

0 

0.00002^69 

0*00000016 

-0.0075 1990 

-0.00004945 

-0.001)05488 

-0 .UOOOO430 

145.765 

• 

2 

0 

0 

0 

1 

0.00000 20 5 

0.00000205 

-0.00062299 

-0. O0O6229Q 

-0.00000487 

-P .00000487 

146.544 

• 

0 

-1 

2 

0 

1 

-n. 00000061 

-0.00000061 

0,00018690 

0*000 1ft 690 

0,00000146 

0 . 00000146 

146.554 

4 

n 

-1 

2 

0 

2 

— 0 • 00000 5ft 8 

-0 .0000058 8 

0.00179108 

0.0017 91 OP 

0.00001397 

0.00001397 

14.7. 355 

. 

-2 

0 

2 

2 

2 

0.00000107 

0.00000107 

-0. 0003 2 7 06 

-0, GOO 3270 6 

-0.00000255 

-0 .000002 55 

147.545 

183.565 

0 

0 

0 

2 

-1 

0.. 00000 419 

-0.0000056? 

-0. 001 27608 

0.0017 121 6 

-0.00000846 

0.00001186 

147.555 

183.555 

0 

V 

0 

2 

U 

0.0000 5025 

-0.00000003 

-0.01530437 

0,00001038 

-0.00011155 

-0.00000740 

147.565 

183.545 

0 

0- 

0 

2 

1 

-0. 00000465 

-0.0 0000629 

0,00141743 

0,0019 1 547 

0.00001127 

n .0000 1465 

148.554 

182.556 

0 

-1 

0 

2 

0 

U. 00000332 

0 .00000 005 

-0.00101199 

-0.00001590 

-O.U0O0O738 

-O.UOOU0U59 

152.656 


1 

1 

2 

-2 

2 

0. 0000007? 

0.0000 0072 

-0*00021801 

-0.000? ) ftOl 

-0.00000165 

-0 ,00000 l 65 

153.645 

, _ 

1 

0 

2 

-2 

L 

0.00000322 

0.00000322 

-0.00098 104 

-0,00098104 

-0,00000743 

-0.00UU0743 

153.655 

177.455 

1 

0 

2 

-2 

2 

0. 00001360 

0 .0000 1483 

-0.00414225 

-0,0045}. 58 3 

-0, 0(101)3147 

-0,00003407 

1 54.656 

176.454 

1 

1 

0 

0 

0 

“0.000001 5? 

-0.00 0000 00 

0.00046706 

0.00(100009 

O.OOOO0339 

o* uonooui 3 

155.435 

175.675 

-1 

0 

2 

0 

0 

-0.00000174 

-0 .oooo oo 00 

0.00052934 

0,(10000010 

0.00000385 

0.00000014 

155,445 

1 75.665 

-1 

0 

2 ‘ 

0 

1 

0 . 0 00008 59 

0.00001156 

-0.00261621 

-0.00351904 

-0.00001995 

-O.UO002627 

155.455 

175.655 

-1 

0 

2 

u 

2 

0. 00005210 

0. 00005681 

-0.01586778 

-0.01729985 

-0*00011989 

-0.000] ? 94? 

155.645 

175.465 

1 

0 

0 

0 

-1 

0.00002566 

-0 • UOU03435 

-0.00781363 

0.01046081 

-6.00005403 

O.UU0O7396 

155.655 

175.455 

1 

0 

0 

0 

0 

0, 000 30 306 

0.0 000000 5 

-0.09229195 

-0,00001667 

-0.00067058 

-0. 1)000242 L 

155.665 

175.445 

1 

0 

0 

0 

1 

0.00002593 

0.00003482 

-0. 00789530 

-0.01060378 

-0,0000 6012 

-0.00007910 

155.675 


1 

0 

0 

0 

2 

-0.00000087 

-0.00000087 

0.00026472 

0.00026472 

0.00000199 

0,00000199 

156.555 

174.555 

0 

0 

0 

1 

0 

-0.00000164 

-0.00000000 

0.00049820 

0.00000008 

O.OOU0O362 

O.OOUOOU12 

156.654 

174.456 

1 

-1 

0 

u 

0' 

0. OUOUO 184 

0.00000000 

-0.00056048 

-0.00000009 

-0.00 0004 07 

-fl .00000014 

157.445 

173.665 

-1 

0 

0 

2 

-1 

0.00000491 

-0.00000654 

-0.00149429 

.0.00199257 

-0*00001040 

. 0*00001 41X . 

157.455 

173.655 

-1 

0 

0 

2 

0 

0.00005787 

0*00000001 

-0.01762385 

-0,00000275 

-0.00012801 

-0 ,00000400 

157.465 

1.73.6.45 

-1 

6 

0 

2 

1 

0.00000542 

0.00000 726 

-rO* 00 165 064 

-0,00221102 

-Q. 0000X249 

-O.OOOQl 643 

158.454 

172.656 

-1 

-1 

0 

2 

0 

0. 00000245 

0.00000000 

-0.00074730 

-o.ooooooii 

-0.00000543 

-0,00000015 

161.557 


.0. 

2 

2 

-2 

2 

-Q. 0.00002 15 ■ 

-0.00000215 

0 ,.00065362*. 

0...QO06.53 62 . 

-.0. 0000 U.4 7.9 

0 ,0000-04.79 

162.556 

16 ft. 554 

0 

1 

2 

-2 

2 

-0. 00005045 

-0.0000 5498 

0.01536357 

0.01674252 

0.0001124ft 

0.00012241 

163.535 

16.7^5 75 

0 

0 

2 

-2 

0 

-0.00000143 

-0.00000000 

. 0.00043592 

0.00000001 . 

0. 0000 03 16 

o. 00000.0 u2._. 

163.. 545 

167.565 

0 

0 

2 

-2 

1 

0. 00000869 

0.0000 1166 

-0.00264678 

-0.00354974 

-O. 0000 1935 

— 0 .00002587 

163.555 

m 

0 

0 

2 

-2 

2 

-0. 0 0000 1 53 

-0.00000153 

0.00046707 

0.00046707 

0 . Q00U0341 

0 « U0OOO34l_ 

163.555 

167*555 

0 

0 

2 

-2 

2 

-0. 0008605? 

-0 .00093797 

0.26206034 

0,28564650 

0.00191293 

0.UO20831 6 

163.557 

167.553 

0 

2 

0 

0 

0 

0.00000113 

o.onoooooo 

-0.00034475 

-0.00000001 

-0.00000250 

-o.uoooomu 



TIDAL ARGUMENT 
CODE NUMBERS 


coefficients of 
LLP F D DM 


163.755 

16 7.355 

2 

0 

0 

-2 

0 

164. 554 

. 

0 

-1 

2 

-2 

2 

164.556 

166.554 

0 

1 

0 

0 

0 

165.545 

165.565 

0 

0 

0 

0 

-1 

l 65.555 

. 

0 

0 

0 

0 

0 

165.555 

, 

0 

0 

0 

0 

0 

. 

165.575 

0 

0 

0 

0 

-2 

. 

173.445 

1 

0 

-2 

2 

-1 

. 

175.475 

1 

0 

0 

6 

-2 

• 

1 35.365 

2 

0 

0 

0 

-1 

. 

135.585 

0 

0 

2 

0 

-1 

. 

193.465 

1 

0 

0 

2 

-1 

• 

1E3.455 

1 

n 

2 

2 

2 


Explanation of sym bols 

T<« LP, F, D and OM arc Brown's fundamental 
arguments; F, D and 0 

Constants 

k = 0 

K # = -7!'J>52430 X 10 3 Julian centurv -1 
K e> = -3"484150 x 10 3 Julian century' 1 
<C-A)/C = 3.272930 x 10 3 
il = 36019856 day' 1 
{Cq/\) - 2 *343862 x 10< 

J 2 = 1.082645 >: 10 3 
(n^/ni^) = 3,334320 X 10 6 


Table 9- 4- (continued) 


COEFFICIENTS FOR FULER ANGLE P ER Tl JR KAT I DnS 
IN SECONDS OF ARC 

COSINES fur DELTA THFTA 

SINES FOR DELTA P5I * S I N ( TH FT A ) 


R nr AT TON 
AXIS 

COS SIN 


o. 000002 66 o.onooonuo 

0.00000732 . 0.00000752 

0.00004335 0 . 00 0 on 0 00 

0.000 31348 -0 .00042084 

0 . 000841 Aft 0.00086146 

0 .00185233 0,00185233 

-0* 00000 7 87 0.000007 87 

0.00000087 -0.00000087 

-0.0000006* 0 . 00000066 

n, oon on ?45 -0.00000245 

0,00000072 -0.00000072 

0.O0000077 -0.00000077 

0.00000061 . -0.00000061 


AXIS OF 
FIGURE 

Cns SIN 


-0. OOOROQ57 -0. 00000002 
-0.00229118 -0 . 007291 1 P 

-0.01320104 -0.000000 1« 

- 0.09546715 0 . 12 R 1614 Q 

-0, 26234698 — 0.26234A9R 

-0.56410116 -0.56410116 

0.00239757 -0,00239757 
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Figure 9.1. The Effect of Polar Motion on Latitude and Time. 


Figure 9.2 shows the diurnal part of the latitude variation for Goddard Space 
Flight Center during the summer of 1970. The polar motion components m x and 
m 2 in Equation (9.14) were computed using all of the coefficients from Table 9.2. 
Of the 135 tidal constituents, those with argument numbers 


135.655 

145.545 

145.555 

163.555 

165.555 
165.565 


(9.16) 
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LATITUDE (SECONDS OF ARC) 




have amplitudes large enough to dominate the motion. In Figure 9.3 the diurnal 
latitude variation based upon all 135 tidal constituents is compared to that ob- 
tained by using only the terms with the argument numbers (9.16). The error due 
to the neglected terms is at most OV0009 or 2.8 cm at the Earth's surface. 



Figure 9.3. Effect of Neglecting all but the 6 Largest Tidal Components in Computing the Diurnal 
Variation of Latitude for Goddard Space Flight Center. Both curves are for a Deformable Earth 
with k = 0.29. 


The effect of rotational and tidal deformation on the diurnal polar motion 
components is shown in Figure 9.4. The diurnal motion of the rotation axis 
within the Earth is affected only slightly by deformation, being decreased in 
amplitude by OV00024 which corresponds to 0.7 cm. In contrast to this, defor- 
mation has a substantial effect on the motion of the angular momentum vector 
within the Earth. As shown in Figure 9.5, the diurnal motion of the angular 
momentum vector is reduced in amplitude by 0 "00645 or 20 cm. The shift in 
the angular momentum vector is due primarily to the factor 



(9.17) 
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which appears in the diurnal part of Equation (7.31) for H/Cfl but is absent from 
the diurnal part of Equation (7.18) for m. Rotational and tidal deformation cause 
the direction cosines of the Earth's axis of figure to oscillate with an amplitude 
of 179 or 58.7 m, as shown in Figure 9.6. The motion of the axis of figure is due 
mainly to the mass redistribution associated with the diurnal tidal bulge. The 
direction cosines of the angular momentum vector are given by Equation (2.24) 
with h - 0 as 


_H_ 

cn 



(9.18) 


and the direction cosines of the axis of figure are, from Equation (E.9), 


= 


c 

(T-~a 


(9.19) 


Therefore 


H 

cn 



+ 



(9.20) 


so that the 58.7 m diurnal motion of «/< f is reduced by the factor (C-A)/C to the 
20 cm departure of H/C Q from the position that it would occupy in a rigid Earth. 

The theory developed in Section 7 is directly applicable to the problem of 
transforming from an inertial coordinate system to a terrestrial system that is 
fixed to a set of observatories in some prescribed manner. The geometry in- 
volved in transforming from inertial to terrestrial coordinates is shown in 
Figure 9.7. The precessional transformation is 


me an of date 


= v- 


z ) 
p ; 




mean of epoch 


(9.21) 


where z , 6 , and £ are the precessional elements written as z ,8, and £ q in 
the Explanatory Supplement to the AENA [1961, p. 29] . The transformation 
from the mean sidereal system of date to the true sidereal system of date is 


X true of date “ ®l( C TD^ X mean of date 

where the true obliquity of date is given by 


(9.22) 


'TD 


= e 


Ae 


rD 


(9.23) 
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The nutations in longitude and obliquity referred to the mean ecliptic of date are 
denoted by and Ae rD respectively. The nutation in longitude is reckoned 

positive westward. The transformation from the true sidereal system of date to 
the terrestrial system is 


terrestrial 


R 2 (-mi) RjOj) R 3 (GASTI) 


rue of date 


(9.24) 


The transformations (9.21), (9.22) and (9,24) are combined to form the complete 
transformation from inertial to terrestrial coordinates. 


In Woolard's development of the theory of precession and nutation, the nuta- 
tions A<A, and Ae for the angular momentum vector referred to the fixed mean 
ecliptic of epoch are obtained as a result of integrating Poisson's equations and 
are given in Woolard’s Table 24 [1953, p. 138] . The Euler angle perturbations 
A0 r and A# r computed from Woolard's Equations (55) are used to convert Ai p H 
and Ae H into the nutations A^ and Ae^ corresponding to the rotation axis. 


A0 r = Ai p H - Si p t 
Ae r = Ae H + id t 


(9.25) 

(9.26) 


In (9.25), Si/» r is subtracted because it is reckoned positive eastward whereas 
Ai// and A0 are reckoned positive westward. The values of Ai/» and Ae 
are listed in Woolard's Table 24 as well as his values for the Euler angle per- 
turbations and id t . The precession and nutation results represented in 
Table 24 undergo a reduction to the mean ecliptic of date. This reduction con- 
sists of adding the corrections in Woolard's Table 25 [1953, p. 1521 to the 
entries in Table 24 so as to produce Table 26 [1953, p. 153] . 


In order to apply the deformable-Earth theory of Section 7 in connection 
with Woolard's solution for the angular momentum vector, his rigid-Earth 
values for i^p T and S0 T must first be removed in Table 24 so as to give the 
precession and nutation of the angular momentum vector referred to the mean 
ecliptic of epoch. Reduction to the mean ecliptic of date is accomplished by 
adding the terms listed in Woolard's Table 25 to the modified entries from 
Table 24. 


The rigid-Earth values of i4 l T and id t from the diurnal terms in Woolard's 
Equations (55) are of order OVOGOOSor 0.15 cm. They are therefore small enough to 
be lost in the numerical round-off error involved in truncating the 5th decimal 
place in Table 24 prior to presenting the nutation series in seconds of arc to 4 
decimal places as Table 26 and as Table 2.5 in the Explanatory Supplement to 
the AENA [1961, p. 44] . The Eulerian terms in Woolard's (55) are of order 
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tf.'0005 or 1.5 cm and are therefore large enough to affect the 4th decimal place 
in the nutation series. The Eulerian terms are not included in Woolard's nuta- 
tion tables, however, because they involve the arbitrary constant of integration 
m Q . The nutation series as presented in Woolard's Table 26 and in Table 25 of 
the Explanatory Supplement to the AENA therefore represent the direction in 
space of the Earth's angular momentum vector to within the 4 decimal places 
given in the tables, and subject to the understanding that the Eulerian terms in 
Equations (7.39) and (7.40) are to be added to the tabular values of S(9 rD and 

S ^rD* 


The Euler angle perturbations 88 and St p of Equations (7.39) and (7.40) 
must be reduced to the mean ecliptic of date before they can be applied to the 
nutation series in Table 2.5 of the Explanatory Supplement to the AENA. The 
Euler angle perturbations S6 l rD and S0 rD referred to the mean ecliptic of date 
are given by 


S *rD + iS ^r D Sin SD = " ie^ 0 * 8 ”) (m - ^ < 9 - 27 ) 

which is analogous to Equation (6.10). It follows from (9.27) and (6.10) that 

S(9 rD + iSi/« rD sin = (Si9 r + iS 0 r sin#) e i( ^ GASTI ”^> (9.28) 

The Greenwich apparent sidereal time is related to the Greenwich mean sidereal 
time 0 M by 


GAST1 r 0 + A0 cos e (9.29) 

M rD TD 

and 0 M is related to the Euler angle <p by Equation (C.4). Equations (9.29) and 
(C.4) are combined to give 


GAST1 - 0 = - a + y cos(r + 4 > » ) cot e (9.30) 

M M 

The planetary precession a is, from Woolard’s Equation (67), 

a = 12M73T - 2."3804T 2 

-0!'00133T 3 (9-31) 

where T is measured in Julian centuries since 1900 Jan 0.5 ET, and the polar 
motion amplitude y is of order O’.'IS. The reduction (9.28) of the Euler angle 
perturbations to the mean ecliptic of date therefore involves adding terms of 
order OVOOOOOOJ and can be neglected in numerical computations. 
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The nutations Ai/> D and Ae rD in longitude and obliquity referred to the mean 
ecliptic of date are given by 


Ax P iD = ~ St A, 


HD 


rD 


A€ r D = A€ HD + S ^rD 


(9.32) 

(9.33) 


Neglecting second order terms, 


80 = ( 9 - 3 4) 

80 _ = he ( 9 * 35 ) 

rD r 

The nutations Ai// hd and Ae HD are taken from Table 2.5 of the Explanatory Sup- 
plement to the AENA. 


It is now possible to write down a transformation from inertial to terrestrial 
coordinates which takes into account the diurnal motion of the rotation axis and 
the angular momentum vector within a deformable Earth. The form of the trans- 
formation is 


ter re s 


trial = v-v R i( m 2 ) RjC-Hj/cn + »!> wen - m 2 > 


R 3 (GAST1 + S^ rD cos e^) 

e HD^ ^3 A '^HD^ 

V-VVW-!,) , (9,36) 

and the associated geometry is shown in Figure 9.7. The transformation from 
the mean system of epoch to the mean system of date is the same as in Equation 
(9.21). Instead of a transformation from the mean equinox of date T M directly 
to the true equinox of dateT T as in Equation (9.22), (9.36) involves a transfor- 
mation from T m to the ascending node of the mean ecliptic of date on the equator 
normal to the angular momentum vector. The obliquity e is given by 

HD 


e HD € M + At HD 


(9.37) 


The Greenwich apparent sidereal time is modified so as to account for the dis- 
placement Si// cose • The transformation from the angular momentum vector 
to the rotation axis is D made using the components of m - H/CH . The diurnal 
terms in m - H/Cfi are computed by differencing the respective components 
given in Table 9.2. The Eulerian part of m - H/C fi is from (9.6) and (9.8), 
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(m - H/Cfi), 


(9.38) 


Eui eri an 





% 


t 


The vector m 0 e icrc,t must be determined from observations. As shown in Fig- 
ure 9.8, it represents the displacement of the rotation pole, with diurnal motion 
filtered out, relative to the mean pole position. Because of the small factor 


00023 


(9.39) 


in (9.38), an uncertainty of, say 1 m, in e^ 0 * will produce an uncertainty of 
only 0.23 cm in (m - H/Ci2L , 

J ' 'Eulcnan 



Figure 9.8. Determination of m^e' ° from Observations. 



The transformation from the rotation axis to the z axis of the terrestrial 
system is based upon the components of the polar motion m. There is as yet no 
reliable model for the "non-diurnal" or "Eulerian" part of the polar motion and 
it is certainly true that the idealized circular motion given by the terms, 


% 





<9.40) 


in (9.6), does not adequately represent non-diumal polar motion in the real Earth. 
In order to represent the polar motion a "semi -analytical'' pole path is deter- 
mined by adding the diurnal terms of (9.6) to an observed pole position with 
diurnal motion filtered out. 


m = m 


'diurnal "non- diu rna I 


(9.41) 


The polar motion transformation is then broken into an empirical part and an 
analytical part as follows : 


R^C-m.) R.(m ) = R (-m . ) R,(m. ,) 

2^ l 7 l v 2 7 2 V l,non- diurnal ; 1^ 2 , non- diu rn al / 

RJ- . ) R.flTL j. . ) 

2 V 1 , diurnal 7 2, diurnal / 

An alternative form of the transformation (9.36) is 


(9.42) 


terrestrial = R 2 (“ m i) W R 3 ( GAST1 ) 

R iC- € td) r 3 (- A ^d> R i( e M> 

R 3 C- Z p) R 2 ( 0 p) R 3(-^>^ean epoch < 9 ‘ 43 > 

where the true obliquity of date e TD is given in terms of the nutation Ae fD in 
obliquity by Equation (9.23). The nutations Ai p tD and Ae rD are given by Equa- 
tions (9.32) and (9.33). The Euler angle perturbations S0. D and S0 rD are com- 
puted from Equations (7.37) and (7.38) with the Eulerian terms determined 
empirically. The diurnal terms in Si9 rD and S^ rD are given in Table 9.4. The 
polar motion transformation can be broken into diurnal and non-diurnal parts 
as in Equation (9.42), 
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A second alternative form of Equation (9.36) is 


\errestrial = R 3 fGAST1 “ “ S ’/' rD > COS € Tl) 

VA) V' A ^ d ) R i( £ m) 

^3^“ Z p^ R 2^p^ R 3^“^0^ X roean of epoch 

where the obliquity of the terrestrial equator is given by 

9 U = e « + Ae tD 

The nutations Ai p and Ae are given by 


A ^D = A ^nn " Sx P 


HD 


*D 


A€ ,D = Ae HD + S ^D 


(9.44) 

(9.45) 

(9.46) 

(9.47) 


Except for second order terms, the Euler angle perturbations S 8 zV) and S</» zD 
are the same as the corresponding perturbations referred to the fixed mean 
ecliptic of epoch. 


H ZJ> = H z < 9 * 48 > 

S# r S 6 (9.49) 

iD x 

Equations (7.43) and (7.44) are used to compute S0 ZJj and S<// zD . The Eulerian 
terms are determined empirically and the diurnal terms are given in Table 9.4. 
No polar motion rotations are needed because the transformation goes directly 
from the mean sidereal system of date to the terrestrial system, bypassing the 
rotation axis altogether. 

The choice of the best alternative form of the transformation from inertial 
to terrestrial coordinates depends upon the particular application at hand. Equa- 
tion (9.36) has the advantage that the nutational part is made using A« P HD and 
Ae HD which come directly from Table 2.5 of the Explanatory Supplement to the 
AENA. Modification of an existing transformation that does not include diurnal 
polar motion and deformable Earth effects is accomplished by inserting the 
transformation 
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(9.50) 


m l , diu rnal ) diurnal ) 

R^-Hj/CO + nij) RjCHj/Cfl - m 2 ) 
R 3 (GAST1 + St/r rD cos e HD ) 


in place of 


R 3 (GAST1) (9.51) 

Equation (9.43) involves fewer individual rotations than (9.36) but the nutations 
Ai// and Ae HD must be modified using (9.32) and (9.33) to form A </> rD and Ae rD . 
Equation (9.44) involves still fewer individual rotations but the absence of a 
polar motion transformation makes it impossible to incorporate the semi- 
analytical pole path defined by Equation (9.41). 
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APPENDIX A 


THE FUNCTIONAL FORM OF DOODSON'S EXPANSION 
OF THE TIDAL POTENTIAL 


Doodson [1922] developed a formal expansion of the tide generating potential 
based upon Brown’s [1905] lunar theory and Newcomb's [1898] theory of the sun. 
Doodson’s tabulated results are represented here in functional form as 

U = J2 [G 205 COS Hoja 1 + ^20j ,) 
j 

+ G 2 03 A 2 0j© COS(o, 2Oj0 t +^ 2Oj0 )] 

+ J2 [G 2U A 21js s in(“ 21j) t^ 21j) + *) 

j 

+ G 2li A 21jo sin ( w 2 i j© t + Aj lj0 + X)] 

+ 2-J ^ G 2 23 A 22js COS (^2213 t + ^22p + 2 A) 


+ G 223 A 22 . 0 cos (<^ 22 j ffl t + fi 2 2je + 2^)) 


+ 2-1 G3031 A 30j3 Sln ^ al 30j3 t +/3 30j3i) 

i 

+ L G313 Salp COS(W 31i 3 t +/3 31j 5 + M 


+ 7 ^323) ^32 


jo ain C fl, 32p t +/? 32P + 2*0 


+ 7 °333 A 33 j3 cos (Wj 3 j 3 t +/ 3 33i3 + 3 \) 
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Solar terms of degree 3 and all terms of degree 4 and higher are neglected be- 
cause they are small. The coefficients denoted by A n m ■ d are the numbers that 
appear in Doodson's tables [1922, pp. 322-325]. The spherical harmonic degree 
and order are denoted respectively by n and m. Individual tabular entries are 
denoted by the index j. The subscript d stands for "disturbing body" which is 
either the moon, ’ , or the sun, o . 

The geodetic coefficients G J are functions of the latitude 0 and Doodson's 
tidal parameter G Dd „ 


where 


G 20d=iG Dd (l-3 S i„2,fc) 

(A.2) 

G 21 d = ^d Sin 2< ^ 

(A.3) 

G 22d = COS ^ 

(A.4) 

G 30d = 1.11803 ^ sin 0(3 _ 5 sin 2 0) 

(A.5) 

G 3 id = 0.72618 G Dd cos 0(1-5 sin 2 0) 

(A.6) 

G 32 d ~ 2.59808 Gpj sin 0 cos 2 0 

(A.7) 

G 33d = God cos3 4 > 

(A. 8) 

3 Gr 2 

^=4 C 3 

(A. 9) 


m d = mass of disturbing body 
c d = mean distance of disturbing body 
G = universal gravitational constant 
r = geocentric radius 

In order to consolidate his results, Doodson adjusts each solar coefficient 
A nmj© so tliat Si yes the correct result when multiplied by the corresponding 
lunar geodetic coefficient G nmS . This adjustment involves the choice of specific 
values for th§_ ratio of lunar to solar mass and the ratio of lunar to solar mean 
distance. If A nnjo denotes the solar coefficient before adjustment, then 
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(A.10) 


vj i _ — vj ri 

nmj nmjO rnn 0 n mj© 

By first dividing out Doodson’s numerical factor so as to recover the A nm . 0 , it 
is possible to incorporate revised values of the lunar and solar masses and mean 
distances into subsequent calculations. The only solar terms in Doodson's table 
are those of degree 2, and for these 


G 2ms _ _ / m e\ / c aV 

G 2m:S \ m 3/ \ C ®/ 


Doodson [1922, p. 318] used 


2m© 

G 2m3 


0.46040 


(A.ll) 


(A.12) 


The solar coefficients X 2mjo for use with solar geodetic_coefficients are there- 
fore given in terms of Doodson’s tabulated coefficients, A 2mj0 , by 


A - 2n, j® 

ft 2mj© 


0.46040 


(A.13) 


The terms of degree 2 and order 1 in the tidal potential are the only ones 
that enter into the torque components (4.19) and (4.20). The corresponding geo- 
detic coefficient (A.3) is written in terms the Legendre function P*(sin$) as 




P* (sin <j>) 


(A. 14) 


The arguments of the trigonometric functions appearing in (A.l) are given 
by linear combinations of the following six standard variables chosen by Doodson 
[1922, p. 310] 



local mean lunar hour angle referred to lower transit 

lunar mean longitude 

solar mean longitude 

longitude of lunar perigee 

longitude of the lunar ascending node 

longitude of solar perigee 
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The standard variables are discussed in greater detail in Appendix B. Doodson's 
arguments are of the form 

dj r + (d a - 5) s + (d 3 - 5) h 

+ (d 4 - 5) p + (d 5 - 5) N' + (d 6 - 5) p s ( A - 15 > 

where d x through d fi are the integers in a code number written as 

d i d 2 d 3' d 4 d s d 6 

The integer d x is always equal to the order m of the harmonic in which the argu- 
ment arises, or equivalently, to Doodson's "schedule number" 0, 1, 2, or 3. The 
standard variables are given in terms of Brown's fundamental arguments by 
Equations (B.l) through <B.5). When terms involving second and higher powers 
of the time are neglected in Brown's arguments, the standard variables are 
expressible as 


W.t + P. (A. 16) 

where t is the Greenwich mean solar time, W. is a frequency and P. is a phase 
angle. In the expression for the local mean lunar time r , the longitude A. is 
separated from the rest of the phase angle so that 

t = W T t + P T + X. (A. 17) 

Doodson's arguments may be written as 


™(V + P T ) + m\ + L (d. - 5) (W.t + P.) 

i = 2 

The frequency and phase of an argument are defined by 

6 

L («, - s) w, 


Cl) . = m W_ + 

nmj d T 


i-2 


o 

4,.^ = + ]2 <d ‘ ■ 5) p ‘ 

i a 2 


(A. 18) 


(A. 19) 


(A.20) 
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and the arguments are written in the form 

"nmjdt +4 1 mjd (A.21) 

which appears in (A.l). 

For the analysis of polar motion and tidal deformation, the tide generating 
potential (A.l) is rewritten as 


U = 




P™(sin <f>) * 


n = 2 m— 0 


•Z> 


COS 

m j a 


“Uid* +/ 3 „ mjd + ^ + (n-m)~ 


(A.22) 


The tesseral diurnal coefficients A 2J jd in (A.21) are related to those in (A.l) by 


A --1a 

21 jd j 21 jd 


(A.23) 
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APPENDIX B 

DEFINITIONS OF THE STANDARD TIDAL VARIABLES 


Explicit relationships are presented here between Doodson's standard var- 
iables, the similar variables defined by Melchior [1966, p. 26] and Brown's 
fundamental arguments as given in the Explanatory Supplement to the AENA 
[1961, p. 44]. The astronomical variables s, h, p, N T , and p s are defined in 
Table B.l and the time variables T, 4 > and t are defined in Table B.2. 

M 


Table B.l 

Astronomical Variables 


Definition 

Symbols Used by 
Doodson [1922, p. 310] 
and Melchior 

[1966, p. 26] 

Symbols 
Used in 
the Expl. 
Supp. 

[1961, p. 107] 

lunar mean longitude, measured 
in the ecliptic from the mean 
equinox of date to the mean as- 
cending node of the lunar orbit, 
and then along the orbit. 

S 


solar mean longitude, measured 
in the ecliptic from the mean 
equinox of date 

h 

L 

mean longitude of lunar perigee, 
measured in the ecliptic from 
the mean equinox of date to the 
mean ascending node of the 
lunar orbit, and then along the 
orbit 

P 

r 

longitude of the mean ascending 
node of the lunar orbit on the 
ecliptic, measured from the 
mean equinox of date 

n 

i 

2 

n 

mean longitude of solar perigee, 
measured in the ecliptic from 
the mean equinox of date 

Ps 

r 
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Table B2 
Time Variables 


Definition 

This 

Paper 

Doodson 
[1922, p. 310] 

Melchior 
[1966, p. 26] 

Greenwich mean sidereal 
hour angle 

<b 


e 

local mean lunar hour 
angle, measured from 
lower transit of the mean 
moon past the local 
meridian 

r 

T 


Greenwich mean lunar 
hour angle, measured 
from upper transit of 
the mean moon past the 
Greenwich meridian 



r 

Greenwich mean solar 
hour angle, measured 
from lower transit of the 
mean sun past the Green- 
wich meridian 


t 


Greenwich mean solar 
hour angle, measured 
from upper transit of 
the mean sun past the 
Greenwich meridian 



t 

Greenwich mean solar 
time 

t 



Greenwich mean solar 
hour angle 

t 

s 
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Doodson's standard variables are given in terms of Brown's fundamental 
arguments, ■i , V , F, D, and Q by 

s = F + n (B.l) 

h = F - D + ft (B.2) 

p = -'e+F + a (B.3) 

N' = - a (B.4) 

/ 

) 

p s = - V + F - D + H (B.5) / 

and the fundamental arguments are given as polynomials in time in the Explana- 
tory Supplement to the AENA [1961, p. 44] . 

Doodson's mean lunar time is reckoned from lower transit of the moon just 
as conventional mean solar time is reckoned from lower transit of the mean sun. 

The relationship between , r, and t is shown in Figure B.l. 

M S 

<£ M = T + S- 7T-k (B.6) 

^ M = t s+ h-7, (B.7) 

From equations (B.6) and (B.7), 

T = t s +h-s+\ (B.8) 

Melchior's mean lunar and solar hour angles r and t are each referred to 
the Greenwich meridian and measured from upper transit. Thus Melchior [1966, 
p. 26] writes 


9 = r + s (B.9) 

9 = t + h (B.10) 

which have the same meaning as (B.6) and (B.7). 

Doodson's standard variables are each expressible in terms of the Green- 
wich mean sidereal hour angle d> u . Equating the forms (A.15) and (A .21) gives 

"nnyd 11 + ^nmjd + mk = d i T + ( d 2 - 5) S + (d 3 - 5) h 

+ ( d 4 - 5) P + (d s - 5) N' + (d 6 - 5) p s (B.ll) 
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G 


Figure B.l. Time Variables 


Substituting for r from (B.6) into (B.ll), with n = 2andm = d 1 =1 gives 

w.t + 0. = (4> u + 7t) + a. (B.12) 

where the subscripts 2,1, and d are omitted and 

a. = (d 2 - 6) s + (d 3 - S) h + (d 4 - 5) p 

+ (d 5 - 5) N' + (d 6 - 5) p s (B.13) 

The frequency of a particular term is written in terms of the Earth's constant 
nominal rotation rate Q [not the Q in Table B.l and Equations (B.l) through 
( 13 . b)] as 

a>. = Q - n. 

i j 

where 

P j = Q ~ " “j (B.14) 
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APPENDIX C 

EULER ANGLES RELATED TO THE MEAN SIDEREAL TIME 


In order to apply Doodson's expansion of the tidal potential to problems in 
polar motion dynamics it is necessary to relate the Greenwich mean sidereal 
hour angle 0 M to the Euler angle 0 defined in Figure 3.1. The difference (0 M - 0) 
is small and is due to polar motion, lunisolar precession and nutation, and 
planetary precession. 

As shown in Figure C.l, the Euler angle 0 is measured in the terrestrial 
equator from its descending node < T > on the fixed ecliptic. Polar motion causes 
0 to depart from the angle 0 r measured in the true equator of date. The differ- 
ence (0 - 0) is expressed in terms of the phase f and amplitude y of the polar 
motion, shown in Figure C.2. Solution of the spherical triangle shown in Figure 
C.3 gives 


0 r - 0 = y COS (F + 0) cot e 1T (C.l) 

in which second order terms in y are neglected. 

The true equator of date moves relative to the mean equator of date because 
of the lunisolar nutation. The nutation A 0 fD in longitude is responsible for the 
difference between mean and apparent sidereal time. From Figure C.l, 

GAST1 = 0 M + A 0 fD cos e TD (C.2) 

where second order terms in A0 rD are neglected. 

Planetary precession of the mean ecliptic causes the true equinoxT T to 
move relative to the equinox V on the fixed ecliptic. From Figure C.l, the 
angle 0 f is related to the Greenwich apparent sidereal time by 

0 r = GASH + a (C.3) 

where second order terms in A0 rD are neglected. The second order terms in 
A 0 rD appear because a is measured in the mean equator of date rather than in 
the true equator of date. 
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V 


90° EAST 



ON THE TRUE EQUATOR OF 
DATE 

Figure C.2. Amplitude and Phase of the Polar Motion 


Equations (C.l), (C.2) and (C.3) are combined to give the difference (0 - 

M 

as 


Al -*= -^rD COS £ TD - 3 

+ 7 cos (P + (^) cot ej-j. (C*4) 

On the right hand side of (C.3) the replacement of e TD and e 1T by e M and the 
replacement of <£ by ^ introduces only additional second order terms. There- 
fore, with second order terms neglected, 

0 M - 0 = - cos e M - a 

+ y cos (r + 0 M ) cot e M (C.5) 
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ECLIPTIC OF 
EPOCH 

Figure C.3. Departure of </> From 0 r Due to Polar Motion 
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APPENDIX D 


GEOPOTENTIAL COEFFICIENTS IN TERMS 
OF THE EARTH'S INERTIA TENSOR 


The Earth’s gravitational potential is expanded in the form 




<t>) 


» n f \n 

P :C Si ^> ^ C nm COS + S nm Sin ^) 

n= 2 m= 1 \ / 

The spherical harmonic coefficients are given by 


(D.l) 


J n = fr n P n ( s in 

m E a E J 


0) dm 


(D.2) 


where 




m K | 
mA. I 


(D.3) 


w — — (m= 1, 2,..., n) 
nm (n + m) ! ' 


(D.4) 


The origin of the x, y, z system is at the Earth’s center of mass so that 

Ji=C„=S,, = 0 (D.5) 


J 2 is given by 


J 2 



(D.6) 


101 



Substituting for r and cf> in terms of x, y and z gives 


J 2 = - f(2z 2 - x 2 - y 2 ) dm (D.7) 

2*e^J 

Equation (D.7) is written in terms of the principal moments of inertia defined in 
Equation (2.3) as 


J* = - -V 1 " + - 2i 33> (D ' 8> 

Similar developments give the other second degree geopotential coefficients in 
terms of inertial integrals as 


C - - 1 T 

21 - r 1 i3 


S - _ 1 T 

b 21 “ “ 1 23 


C 22 = — ^22 - I,l> 


4 V 


(D.9) 


(D.10) 


(D.11) 


S 


22 


1 

2 m E a E 


12 


(D.12) 


Equations (D.8) through (D.12) are rewritten in terms of the inertia tensor per- 
turbations defined in Equation (2.5) as 


_ £ - A + 2c 3 3 C 1I " C 22 

\ a l 


(D.13) 


C 


21 



(D.14) 
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(D.15) 


S 21 = " 


'23 


^*8 


C 22 = 


C 22 " C 11 


4m E^ 


(D.16) 


S 


22 



(D.17) 


The trace of the inertia tensor (2.5) is invariant under a small deformation. In 
order to show this, the trace is written as 


2^I jj =2|r 2 dm (D.18) 

The mass redistribution resulting from a small deformation is treated as a sur- 
face layer as shown in Figure D.l. The mass element dm is written in terms of 
a surface density p{6, k) as 


dm - p(<p, k) ds (D.19) 

where ds denotes a surface element. The density is expanded in spherical sur- 
face harmonics 


where 


p(i>. \) =2^ s „ 

n=0 


(D.20) 


S o = 0 (D.21) 

in order that mass be conserved. If V Q and S denote the volume and surface 
corresponding to the sphere of radius a , then 

E 


J*r 2 dm= J* r 2 dm + ^ 


a |p(0> k) ds 


(D.22) 
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APPENDIX E 


DIRECTION COSINES OF THE AXIS OF FIGURE RELATED 
TO PERTURBATIONS IN THE INERTIA TENSOR 


The Earth’s inertia tensor is written in the form 


(~ 



A + c n 

C 1 2 

C 13 

I = 

C 12 

A + C 22 

C 23 


C 13 

C 23 

C + C 3_3 


(E.l) 


given by Equation (2.5). The perturbations c^ are small in relation to the 
principal moments of inertia so that the xyz system of Figure 2.1 is almost a 
principal axis system. The xyz system must be rotated through the small 
angles l and p shown in Figure E.l in order to make the z' axis coincide with 
the principal axis of inertia. 


* z’ 



Figure E.l. Rotation of the terrestrial System to Make the x axis 
Coincide With the Principal Axis of Inertia 


105 



Neglecting second order terms in -t and p, the coordinate transformation 
from xyz to x'y'z* is 


R = 


1 0 


0 1 


-V 

- P 
0 


The inertia tensor transforms as 


I ' = RIRT 


The transformed inertia tensor is 


I' = 


A + c 


ii 


'12 


t^(A - C) + c l3 ] 
The direction cosines 


tf(A - C) + c 13 ] 
'22 [ PCA - C) + c 23 ] 

[p(A - C) + c 2J ) C + c 33 


'-12 

A + c. 


(E.2) 


(E.3) 


(E.4) 


i = 

C- A 


(E.5) 


'23 


c - 


(E.6) 


will diagonalize I’ except for the c 12 terms. If c 12 / 0 an additional rotation 
about z r is necessary to produce a principal axis system. 

The direction cosines of the principal axis of inertia are combined into a 
complex number 0 f called the axis of figure 

0 f = -6 + ip (E.7) 

In terms of the complex representation, 

C=C 13 + iC 23 (E ’ 8) 
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of the products of inertia, the axis of figure is given by 



APPENDIX F 

SYMMETRIC TIDAL ARGUMENTS 


Summations of the form 


(F.l) 

j 

arise in deriving the formulas for the Euler angle perturbations in Sections 6 
and 7. The A. denote general coefficients,^-. t +/3, . is a tidal argument and 4> 
is the Euler angle defined in Figure 3.1 which describes the Earth’s diurnal 
rotation. Two tidal arguments with distinct indices j and j_ are called sym- 
metric when 


co. t + /6. . = (d>, + tt) + a. 
j+ j+ vy m * i 

Wj.t + /3._ = (0 M + tt) - a. 

where a. denotes a linear combination of Doodson’s arguments and ^is the 
Greenwich mean sidereal hour angle. 

The arguments of symmetric terms in (F.l) are 

"j ± t + /3. ± - 0 = (0 M - 0) + 77 ± a. (F.4) 

Two symmetric terms are combined to form 

X. + e~ i<w »+ t+ *J ± ”* ) + A. e _i <‘V t+/ 3 j-'* ) 

= [-cos (0 M - 0) cos a. (A j+ + A^> 

- sin (^ M ~ 0) sin a. (-A j+ + A._)] 

+ + i [-cos(0 - 0) sin a (-A.. + A._) 

+ s in(0 M - 0) cos a. (A. + + A._)] (F.5) 


(F.2) 

(F.3) 
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